一种确定微米毛细管通道中压力与流量关系的方法与流程
本发明涉及石油强化开采领域，特别是涉及一种利用具备微米毛细管通道的动态测量装置获取气液动态驱替相界面的位移、接触角和速度后确定微米毛细管通道中压力与流量关系的方法。
背景技术
单个微米尺度的毛细管中的两相驱替机理是自然环境多孔介质内的两相驱替的基础，多孔介质中的两相驱替发生在众多的工业或自然过程中，比如致密油藏的水驱气、溶解气驱、水驱油等过程，上述的工业过程的特点是粘性力微弱(毛管数的范围10-10<ca<10-5)，重力微弱(邦德数bo小于10-4)。影响驱替过程的关键因素为毛管力，而计算毛管力的关键参数是两相界面与管道壁面的接触角。毛管力的计算按照young-laplace公式计算，如下式所示，
γ是界面张力，接触角θ的实际范围为0～180°，当θ<90°时pc为正值，当θ>90°时，pc为负值，因此毛管力有可能是阻碍驱替的阻力也可能是推进驱替的动力。不仅如此，接触角评估的错误将会使毛管力计算存在严重的误差，从而导致驱替过程的评估失误。因此，获取上述工业条件下的动态接触角信息对于评估整个多孔介质内两相驱替过程至关重要。
现已较完善的测量相界面运动状态下(动态)接触角的技术，比如平板浸没法、毛细管液升法。平板浸没法是将一块待测固体平板，以一定速度、以平板垂直于液面的姿态进入或抽离液面，这个过程采用高速摄像机拍摄液面与平板接触点处相界面的接触角，在通过图像测量接触角的大小；或者采用测力计测量平板在这个运动过程中的受力，测算出界面张力反推出接触角。
毛细管液升法，是将一根已知半径的毛细管入口垂直伸入湿润相液体的液面下，通过高速摄像机拍摄湿润相液面上升的高度及对应的时间，通过washburn公式推算出接触角。
但是采用上述两种方案都不能够测量毛管数10-10<ca<10-5及邦德数bo<10-4工况下的动态接触角。
此外，为了观察到接触角，需要采用高倍数的显微镜获得放大的图像，也就意味着观察的视野小。然而动态的驱替过程意味着相界面会移动一定的距离，因此高倍数的显微镜无法观察到相界面的运动过程。如果采用低倍数的显微镜进行观察，可以观察到相界面位置的变化，却无法观察到相界面，也就无法直接测量接触角。
综上所述，现需要一种可清晰观察到相界面的动态测试平台，以求出相界面的接触角、位移和速度，进而得到压力与流量之间的关系，从而简化相界面的控制过程的方法。
技术实现要素：
本发明的目的是要提供一种确定微米毛细管通道中压力与流量关系的方法。
特别地，本发明提供一种确定微米毛细管通道中压力与流量关系的方法，包括如下步骤：
步骤100，将测量选定的微流动芯片安装在动态测量装置移动平台的移动座上并连接注射泵，调节并安装摄像单元，同时使各设备与控制系统连接；先将测试液体注入微流动芯片上的微米毛细管通道中，测量出当前测试液体的压力，由控制系统控制抽取泵的注入压力满足预定的测试压力；
步骤200，在抽取泵控制压力的过程中，控制系统同时控制移动平台的电机使移动座水平移动，带动微流动芯片向与微米毛细管通道中测试液体移动方向相反的方向移动，且两者的移动保持动态平衡并使测试液体形成的相界面始终位于摄像单元的视野中心；
步骤300，在微米毛细管通道中的测试压力满足要求时，由摄像单元获取相界面的图像并传送至控制系统；控制系统对该图像进行灰度转化和滤波处理，然后利用阈值分割算法将灰度图像转化成二值图像，利用连通区域方法通过像素点的颜色提取出构成管道形边界图像的边界点；
步骤400，识别边界点中代表管道形边界图像中表示左右两个相对界面的相界面和末端面及上下两个相对壁面的坐标点，并根据坐标值的大小对各坐标点进行排序和编号，依据构成未端面的坐标点特性，将所有未端面的坐标点删除，同时确定两个壁面与未端面相交的两个末端接触点坐标；
步骤500，对剩余坐标序列利用能量函数根据坐标值的变化对各坐标点进行排序，初步确定相界面和两个壁面的坐标点分布趋势，然后设定能量阈值选取变化低于该能量阈值的坐标点并进行划分和编号，通过取编号的平均值即可确定相界面两个端点的粗略坐标，进而初步划分出分别构成相界面和两个壁面的坐标点；
步骤600，对构成相界面的坐标点进行三次多项式拟合后得到三次函数，对两个壁面的坐标点进行一次函数拟合得到两条对应的表达式，求解后即可得到相界面两个端点的精确坐标，再将两者之间的所有坐标点由三次函数求出，拟合两个端点之间的所有坐标点，即得到相界面；将精确坐标分别代入三次函数的一次导数中，即得到相界面两端点处的接触角；
步骤700，已知图像分辨率的情况下，计算出相界面两个端点到视野一侧的距离，而相界面到该侧的距离为两个端点到该距离的平均值，相界面某一时刻至微米毛细管通道入口的位移在前述基础上利用移动座的位移速度和距离的加和即可得到；而相界面的速度是相界面位移的一阶导数，对位移以中心差分格式进行求导即得到相界面的当前速度；
步骤800，将相界面速度、相界面接触角、相界面位移代入计算微米毛细管通道内压力的两相流方程求出当前微米毛细管通道内压力，再根据流量和相界面速度之间的关系，即得到当前压力和流量之间的关系，进而实现通过流量控制压力的目的。
在本发明的一个实施方式中，所述移动座的位移速度获取方式如下：设所述移动座的速度为v，根据：
n＝f/k
其中，n为电机的转速，f为控制系统发出的脉冲信号频率，k为转频比；则：
v＝lpn
lp为已知丝杠的螺距。
在本发明的一个实施方式中，所述步骤300中提取出构成管道形边界图像的边界点的具体处理步骤如下：
步骤310，将二值图像中黑色区域的像素点记为数值0，白色区域的像素点记为数值1；
步骤320，将所有像素值i(xi,yi)＝1的坐标点同其相邻的像素点进行“与”运算；若运算结果为1则将像素点标记为内部点，若运算结果为0则将该像素点标记为边界点；
步骤330，运算结束后将所有边界点提出来，即得到所有构成管道形边界图像的边界点。
在本发明的一个实施方式中，所述步骤400中，构成末端面的坐标点特性如下：
由于未端面为一条竖线，因此其横坐标值变化一致或在一定范围内，而纵坐标变化为渐增，根据所有坐标点中横、纵坐标值的大小对各坐标点进行从小到大的排序，选择横坐标值变化在预定范围内的坐标点全部删除，即删除了未端面。
在本发明的一个实施方式中，所述步骤500中，初步划分出分别构成相界面和两个壁面的坐标点的过程如下：
步骤510，从任意一个坐标点
步骤520，将区域序列中编号最小者作为上壁面的一个端点，其与对应的上壁面的未端接触点坐标之间的坐标点为形成上壁面的坐标点编号范围；
将区域序列编号中最大者作为下壁面的一个端点，其与对应的下壁面的未端接触点坐标之间的坐标点为形成下壁面的坐标点编号范围；
对区域序列中所有编号求和后取平均值，则第一个小于平均值的编号为相界面的一个端点，而第一个大于平均值的编号则为相界面的另一个端点，两个端点之间的坐标点则为形成相界面的坐标点编号范围。
在本发明的一个实施方式中，还包括对坐标点排序时的顺序修正程序：
设置边界点序列的能量函数，由两个壁面的末端接触点坐标之间的区域且根据两个未端接触点坐标的纵坐标的差值确定极大值，即emax＝|yd-yc|，再通过寻找单个最大值算法确定ei＝emax的坐标编号，设其编号为i＝imax，将坐标点序列中的第一坐标点分别与能量函数进行循环位移，将编号imax变成编号0，仍然是n+1个坐标点，从0开始编号，新的编号j和原编号i的关系为：
通过对排序的坐标点进行循环位移方式调整排序后，使得随着编号的递增，代表管道形边界图像的四个边界点依次出现。
在本发明的一个实施方式中，其特征在于，
所述步骤500中，所述初步确定相界面和两个壁面的坐标点分布趋势的过程如下：
将坐标点序列上的任意一个坐标点(xi,yi)，与其在序列上相距m位的坐标点(xi+m,yi+m)按照下式进行比较：
ei＝||xi-xi+m|-|yi-yi+m||
可得到相界面上的坐标点在x方向上的变化较弱，在y方向上变化较强，因此ei≈m；而两个壁面上的坐标在x方向上的变化较强，在y方向上变化较弱，会有ei≈m；而从两个壁面过渡到相界面的这部分边界点其在x方向和y方向上的变化基本相同ei≈0。
在本发明的一个实施方式中，所述步骤600中，得到相界面两个端点的精确坐标的步骤如下：
步骤610，设构成相界面的第一个坐标点为p，最后一个坐标点为q，用p至q范围内所有坐标点进行三次多项式拟合，以横坐标x为因变量，纵坐标y为自变量，拟合得到弯月面形状的三次函数x＝fpq(y)；
步骤620，设构成上壁面的第一个坐标点为d，最后一个坐标点为o；构成下壁面的第一个坐标点为r，最后一个坐标点为c；采用一次函数拟合，以横坐标y为因变量，纵坐标x为自变量，得到两条直线的表达式y＝kdox+cdo及y＝krcx+crc，其中kdo和krc为直线的斜率；
步骤630，根据拟合得到的三次函数及两条直线表达式，求解以下两个方程组；
即得到相界面上端点a的坐标(xa,ya)以及相界相下端点b的坐标(xb,yb)。
在本发明的一个实施方式中，其特征在于，获得接触角的过程如下：
根据相界面的拟合函数x＝fpq(y)，将上端点a的坐标(xa,ya)以及下端点b的坐标(xb,yb)分别代入fpq的一次导数中即可计算出接触角θa和θb：
进行倾斜修正：
设两个壁面的倾斜角度为上壁面和下壁面倾斜角度的平均值，即(arctankdo+arctankrc)/2；则修正后的接触角θa和θb为：
在本发明的一个实施方式中，其特征在于，
利用加和得到相界面位移x的过程如下：
其中，da为上端点至视野一侧的距离，db为下站点至视野同侧的距离，s为移动座的位移，d为相界面与视野同侧的距离，i＝1,2,3…m，δt为移动座的速度记录的时间间隔。
在本发明的一个实施方式中，相对相界面的位移x，相界面的速度v是位移x的一阶导数
由此得到相界面的速度。
在本发明的一个实施方式中，其特征在于，
所述步骤800中得到当前压力和流量之间的关系的过程如下：
两相流方程为：
pt为微米毛细管通道内的希望压力值，v为相界面速度，x为相界面位移，θ为相界面接触角，μ为测试液体粘度，s为微米毛细管通道截面边长，a为微米毛细孔通道截面面积，r为微米毛细管通道的等效半径，γ为界面张力；
设注射泵的注射流量为q，则：由于q＝av，因此上式可改写为：
由此可直接计算出理想的注射流量q。
在本发明的一个实施方式中，还包括对所述移动座的位移记录步骤：设移动座从0时刻到t时刻的位移为s，根据位移的公式：
控制从时刻0开始，每隔δt的时间间隔记录移动座的速度vi，i＝1,2,3…m，到时刻t时总共记录了m个速度数据，则该式可以离散成：
其中，δt的取值依赖于实际操作过程中所述移动座的速度是否会有快速的变化。
在本发明的一个实施方式中，使相界面始终自动维持在摄像单元的视野宽度的中间的方法如下：
其中，v为实时采集的相界面速度，d为相界面在摄像单元的视野宽度中距离左边缘的距离，v为移动座的速度，τv为时间常数，其物理意义为从将相界面偏离摄像单元的视野中央的距离d-ws/2降到0所需要的时间。
在本发明的一个实施方式中，所述动态测量装置包括：
微流动芯片，用于供测试液体通过，设置有微米毛细管通道；
压力控制管路，包括通过管路与所述微流动芯片连接的注射泵，在管路上安装有检测压力的压力传感器；
移动平台，包括作为支撑的基座，和安装在基座上的丝杆，通过螺纹套在丝杆上沿直线移动的移动座，驱动丝杆转动的电机，所述微流动芯片安装在移动座上；
摄像单元，设置在所述微流动芯片的上方，用于获取测试液体的相界面图像并输出；
控制系统，根据所述压力传感器的信息控制所述注射泵的注射过程，同时通过控制所述电机使所述微流动芯片的水平移动速度与所述微流动芯片内的测试液体移动速度相同且方向相反，以使测试液体形成的相界面始终保持在所述摄像单元的视野范围内。
在本发明的一个实施方式中，所述微流动芯片为透明材料制成的板状结构，一面上设置有内凹的凹槽，在所述凹槽上盖有贴合所述凹槽上表面的盖板，所述微米毛细管通道设置在所述凹槽内且以阵列方式设置有多个，每个所述微米毛细管通道的两端分别设置有进液区和出液区；所述微米毛细管通道的宽度为10～100μm，所述进液区和所述出液区的宽度至少为所述微米毛细管通道宽度的15～20倍。
在本发明的一个实施方式中，所述进液区和所述出液区为水滴形且水滴形的尖端与所述微米毛细管通道连通，在所述进液区的中间设置有底部开口的进液孔，在所述出液区的中间设置有底部开口的排液孔；在所述进液区和所述出液区内分别均匀地设置有多根垂直的防止微米塌陷的立柱，所述立柱以环形的方式逐圈布置，且中间留出所述进液孔和所述出液孔的设置空间。
在本发明的一个实施方式中，所述移动座的上表面安装有两个相对间隔且垂直的支撑板，在两个所述支撑板的相对面上分别安装有固定所述微流动芯片的固定台；在所述移动座上设置有多个安装孔，所述支撑板通过l形折片利用所述移动座上的安装孔与所述移动座固定，在所述支撑板上设置有垂直于放置后所述微流动芯片的滑动槽，所述固定台通过所述滑动槽可调节的安装在所述支撑板上。
在本发明的一个实施方式中，所述固定台为条状的支撑条，或为设置有夹持所述微流动芯片侧边的凹槽的夹持条；在两个所述支撑板之间安装有为所述微流动芯片照明的光源。
在本发明的一个实施方式中，在所述移动座的一侧通过l形折片安装有带三通的透明接头，其中两个通孔为水平且相对设置的直通，另一个通孔为垂直于两个直通中部的竖通，在竖通上通过螺纹安装有与所述控制系统连接的压力传感器，在两个直通上通过螺纹连接有鲁尔接头，所述注射泵通过柔性软管与所述鲁尔接头连接，另一直通的鲁尔接头通过柔性软管与所述进液区的进液孔连接，所述出液区的出液孔通过柔性软管与接水盒连接。
本发明实现了对微米尺度的毛细导管内的气液两相流动过程的压力或流速的控制，给出了获取相界面和相界面位移的方法，首次实现了在低毛管数、低邦德数状态下的相界面形状获取。第一次采用了基于图像处理方法的相界面形状计算和速度计算，极大提高了数据处理的效率及测量精度，通过对相界面的速度、接触角和位移的确定，进而可确定微米毛细管通道内压力和测试液体流量之间的关系，从而可使控制系统根据流量直接控制和调整压力变化，简化实验过程和移动座的速度控制算法，提高了实际操作时的灵活性，可满足各种工况的需求。
本发明的优点在于鲁棒性强，测量结果稳定，具有较强容错性，对于质量不好的图像也能提供准确的测量结果；采用该计算方法编程简单，程序运行效率高，采用常规处理器处理一张2300×1730像素的图片耗时小于1s。
附图说明
图1是本发明一个实施方式的方法流程示意图；
图2是本发明一个实施方式的动态测试装置的结构示意图；
图3是本发明一个实施方式的移动平台结构示意图；
图4是微流动芯片的整体结构示意图；
图5是判断任意一个坐标点是否是边界点的示意图；
图6是识别出的形边界图像示意图；
图7是相界面、上下壁面和末端面的示意图；
图8是识别并提出的形边界图像示意图；
图9是形边界序列上的坐标点排列及各个顶点的示意图；
图10是能量函数算法的示意图；
图11循环位移前的坐标点序列示意图；
图12相界面和两个壁面的拟合以及接触点的确定示意图。
图13相界面拟合示意图；
图14是微米毛细管的进液区结构示意图；
图15是图3中接头的结构示意图。
具体实施方式
如图1所示，本发明提供一种确定微米毛细管通道中压力与流量关系的方法，一般包括如下步骤：
步骤100，将测量选定的微流动芯片安装在动态测量装置移动平台的移动座上并连接注射泵，调节并安装摄像单元，同时使各设备与控制系统连接；先将测试液体注入微流动芯片上的微米毛细管通道中，测量出当前测试液体的压力，由控制系统控制抽取泵的注入压力满足预定的测试压力；
如图2、3、4所示，本实施方式涉及的动态测量装置一般包括微流动芯片20、压力控制管路40、移动平台30、摄像单元50和控制系统10。
该微流动芯片20用于供测试液体通过，设置有微米毛细管通道21；具体的微流动芯片20可以是玻璃、有机材质等制作的板形结构。微米毛细管通道21是气液动态两相驱替发生的场所，其需要具有足够的强度、加工精度和可操作性。根据本实施方式的需要，微米毛细管通道21的特征尺寸在10～100μm范围内，如本实施方式采用的为20、40和80μm分别进行测试。
具体的微流动芯片20可以采用聚二甲基硅氧烷(pdms)材料，而微米毛细管通道21则可以采用软刻蚀工艺制作，微米毛细管通道21的加工截面可以为矩形或圆形。本实施方式采用正方形，实验中的条件为：1.刻蚀深度为20微米(即高度为20微米)，长度以满足相界面的足够移动距离即可；2.测试液体为交联剂，浓度为10％；3.键和温度为100摄氏度；4.硬化时间为120分钟。工艺要求2、3和4可以保证结构的弹性模量为2000kpa，可满足小于20kpa的实验要求(弹性应变小于1％)。
压力控制管路40包括通过管路与微流动芯片20连接的注射泵41，在管路上安装有检测压力的压力传感器42。压力传感器42用于测量注射泵41施加至管路内的压力，即测试注射泵41注射流量产生的实时压力并传递给控制系统10。注射泵为高精度注射泵，而且在控制系统的控制下实现注射速率的调整。
该移动平台30包括作为支撑的基座32，和安装在基座32上的丝杆34，通过螺纹套在丝杆34上沿直线移动的移动座33，驱动丝杆34转动的电机31，微流动芯片20安装在移动座33上。整个移动平台30相当于是一个丝杆结构，通过电机31驱动丝杆34的转动来使移动座33水平移动，进而使微流动芯片20实现水平移动。
摄像单元50设置在微流动芯片20的上方，用于获取测试液体在微米毛细管通道21中的相界面图像，然后输出至控制系统10进行接触角等数据分析。摄像单元50可以采用高精度的数码相机。
控制系统10可以是一个带有控制程序的计算机，其根据压力传感器37的信息控制注射泵41的注射过程，以调节压力；同时通过控制电机31使微流动芯片20的水平移动速度与微米毛细管通道21内的测试液体移动速度相同且方向相反，以使测试液体形成的相界面始终保持在摄像单元50的视野范围内，既实现高精度观察，又具备足够的观察视野。
各设备之间的配合工作过程如下：
先将设置好的微流动芯片20水平安装在移动座33上，数码相机与计算机连接后安装在微流动芯片20的正上方，调整好拍摄距离；将相应的测试液体注入微流动芯片20的微米毛细管通道21内，通过计算机对注射泵41的工作过程进行控制，使注入压力相对实验要求的微米毛细管通道21内的预定目标压力达到一致；同时计算机需要控制电机41的转动来使移动座33带动微流动芯片20移动，由于测试液体的相界面会随着注射泵31的压力波动而在微米毛细管通道21内来回移动，影响数码相机的观察；此时，可以控制微流动芯片20相对相界面的移动而同步移动，但每次移动的方向与相界面的移动方向相反，使相界面虽然移动但被动态地限制在数码相机的观察视野内，从而保证既有足够的观察视野又可获取清晰的相界面图像。
本实施方式的动态测量装置可以通过相互独立的移动结构和气液两相驱替结构实现气液相界面当前停留位置的自动控制，使显微镜能够始终观察到最合适的相界面图像，提高了实际操作时的灵活性以便满足各种工况的需求。
当前测试液体的压力可由压力传感器42得到。
步骤200，在抽取泵控制压力的过程中，控制系统同时控制移动平台的电机使移动座水平移动，带动微流动芯片向与微米毛细管通道中测试液体移动方向相反的方向移动，且两者的移动保持动态平衡并使测试液体形成的相界面始终位于摄像单元的视野中心；
控制移动座速度v的方式如下：
n＝f/k
其中，n为电机的转速，f为控制系统发出的脉冲信号频率，k为转频比；
v＝lpn
lp为已知丝杠的螺距。由此，移动座的移动速度可由控制系统进行精确控制，电机的控制器可以接受控制方向的电压信号，设比较电压为2.5v，给控制器输入一个高电压(如5v)，则电机将会正转，给控制器输入一个低电压(如1v)，则电机将会反转。因此移动座的移动速度v和方向都是可以由计算机进行实时控制和记录的。具体移动座往某个方向移动时速度v为正数，反之v为负数则可以由人为规定。
记录移动座的位移对计算相界面的位移和速度很重要，根据上述控制步进电机的思路，设计了如下的移动座位移记录方案：设移动座从0时刻到t时刻的位移为s，根据位移的公式：
控制从时刻0开始，每隔δt的时间间隔记录移动座的速度vi，i＝1,2,3…m，到时刻t时总共记录了m个速度数据，则该式可以离散成：
其中，δt的取值依赖于实际操作过程中所述移动座的速度是否会有快速的变化。只要时间间隔δt取得足够小，如10ms～100ms，则可以获得高精度的位移数据。δt的取值依赖于实际操作过程中移动座的速度是否会有快速的变化。
使相界面始终自动维持在摄像单元的视野宽度的中间的方法如下：
其中，v为实时采集的相界面速度，d为相界面在摄像单元的视野宽度中距离左边缘的距离，v为移动座的速度，τv为时间常数，其物理意义为从将相界面偏离摄像单元的视野中央的距离d-ws/2降到0所需要的时间。
自动控制模式下，当相界面以某一个速度v偏离视野中央，移动座需要采用大于该速度的反向速度将相界面维持在视野中央，并且当相界面回到视野中央后，移动座的速度应当为v＝-v。
步骤300，在微米毛细管通道中的测试压力满足要求时，由摄像单元获取相界面的图像并传送至控制系统；控制系统对该图像进行灰度转化和滤波处理，然后利用阈值分割算法将灰度图像转化成二值图像，利用连通区域方法通过像素点的颜色提取出构成管道形边界图像的边界点；
如图6所示，摄像单元获取的无论是彩色图像还是灰度图像，都进行转化，以仅得到反映相界面的微米毛细管通道的图像。
图像中的白对应灰度值255，黑对应0，其余不同程度的灰色的灰度值在0和255之间。滤波处理是为了去除图像中像素的噪点、小面积的独立斑块。
其中提取出构成管道形边界图像的边界点的具体处理步骤如下：
步骤310，将二值图像中黑色区域的像素点记为数值0，白色区域的像素点记为数值1；
这里采用阈值分割算法将灰度图像转化成二值(黑白)图像，此时液体区域和固体区域被划分为黑色区域(数值为0)，空气所占据的区域为白色区域(数值为1)，因此图像中存在两个区域黑色区域及白色区域。
步骤320，将所有像素值i(xi,yi)＝1的坐标点同其相邻的像素点进行“与”运算；若运算结果为1则将像素点标记为内部点，若运算结果为0则将该像素点标记为边界点；
如图5所示，这里是通过识别连通区域来提取边界，先将所有像素值i(xi,yi)＝1的坐标点同其相邻的8个(若在图像边缘则为5个，若在图像四个角落中则为3个)像素点的像素值为1进行“与”运算；若运算结果为1则将像素点标记为内部点，若像素点相邻的8个坐标点中存在一个像素值i＝0(黑色点)的像素点，则运算的结果为0，该坐标点标记为边界点。
步骤330，运算结束后将所有边界点提出来，即得到所有构成管道形边界图像的边界点。如图7、8所示。
步骤400，识别边界点中代表管道形边界图像中表示左右两个相对界面的相界面和末端面及上下两个相对壁面的坐标点，并根据坐标值的大小对各坐标点进行排序和编号，依据构成未端面的坐标点特性，将所有未端面的坐标点删除，同时确定两个壁面与未端面相交的两个末端接触点坐标；
对提取后各坐标点进行排序的方式如下：从任意一个边界点
这里的编号是指从0开始按1、2、3……的顺序对排序后的坐标点依次进行编号。
构成末端面的坐标点特性如下：
如图9所示，由于未端面为一条竖线，理论上应该是一条垂直的直线，但实际获取的图像中其可能会产生一定的偏差，因此其横坐标值变化一致或在一定范围内，而纵坐标变化为渐增，根据前面获取的坐标点中横、纵坐标值的大小对各坐标点进行从小到大的排序，选择横坐标值变化在预定范围内的坐标点全部删除，即删除了未端面。
本发明中的壁面指显微图像中微毛细管道的内壁，相界面指气体和液体的交界面，接触点指相界面与上下两个壁面的交点，未端面则为与相界面相对的管道末端。具体地，设管道末端的横坐标为x＝ws，对于任意边界点i，只要满足xi＝ws，则将点i从边7界点的序列中剔除，则如图所示，线段c(末端面)会被剔除，而剩余的点中与上壁面(线段d)和下壁面(线段b)两端相邻的末端接触点坐标d和c则可以获得。
步骤500，对剩余坐标序列利用能量函数根据坐标值的变化对各坐标点进行排序，初步确定相界面和两个壁面的坐标点分布趋势，然后设定能量阈值选取变化低于该能量阈值的坐标点并进行划分和编号，通过取编号的平均值即可确定相界面两个端点的粗略坐标，进而初步划分出分别构成相界面和两个壁面的坐标点；
初步确定相界面和两个壁面的坐标点分布趋势的过程如下：
设边界点中除去线段c之后还剩余n+1个像素点，坐标组成的坐标点序列为(x0,y0)、(x1,y1)…(xn,yn)，在实际图像中，相界面并不是一条光滑的曲线，上壁面和下壁面同样也不是笔直的直线。通过能量函数来计算坐标点的变化趋势，可将相界面和上壁面、下壁面进行粗略区分。
将坐标点序列上的任意一个坐标点(xi,yi)，与其在序列上相距m位的坐标点(xi+m,yi+m)按照下式进行比较：
ei＝||xi-xi+m|-|yi-yi+m||
可得到相界面上的坐标点在x方向上的变化较弱，在y方向上变化较强，因此ei≈m；而两个壁面上的坐标在x方向上的变化较强，在y方向上变化较弱，会有ei≈m；而从两个壁面过渡到相界面的这部分边界点其在x方向和y方向上的变化基本相同ei≈0。
初步划分出分别构成相界面和两个壁面的坐标点的过程如下：
步骤510，从任意一个坐标点
比如编号10、11、12、20、21、22中，10、11和12是连续的，20、21和22是连续的，但是12到20是不连续的。
步骤520，将区域序列中编号最小者作为上壁面的一个端点，其与对应的上壁面的未端接触点坐标之间的坐标点为形成上壁面的坐标点编号范围；
将区域序列编号中最大者作为下壁面的一个端点，其与对应的下壁面的未端接触点坐标之间的坐标点为形成下壁面的坐标点编号范围；
对区域序列中所有编号求和后取平均值，则第一个小于平均值的编号为相界面的一个端点，而第一个大于平均值的编号则为相界面的另一个端点，两个端点之间的坐标点则为形成相界面的坐标点编号范围。
如图10所示，假设这两段连续的区域的首末两个坐标点的编号为总共u个坐标点(如编号分别是10、11、12、20、21、22，此时u＝6)；搜索序列中最小值(例子中的编号10)，令此点为o点；搜索序列中最大值(例子中的编号22)，令此点为r点；将序列所有编号进行求和运算并处以u(在例子中为(10+11+12+20+21+22)/6＝16)，令计算结果为
前述得到的坐标点序列不一定按照上壁面右端点(d)-相界面上端点(a)-相界面下端点(b)-下壁面右端点(c)的顺序排列，为了避免保证能按前述方法正确获得点o、p、q和r，需要在采用循环位移的方式，调整坐标点点序列的排序，使得随着编号的递增，关键点d-a-b-c依次出现。具体的对坐标点排序时的顺序修正程序如下：
如图11所示，设置边界点序列的能量函数，由两个壁面的末端接触点坐标之间的区域且根据两个未端接触点坐标的纵坐标的差值确定极大值，即emax＝|yd-yc|，再通过寻找单个最大值算法确定ei＝emax的坐标编号，设其编号为i＝imax，将坐标点序列中的第一坐标点分别与能量函数进行循环位移，将编号imax变成编号0，仍然是n+1个坐标点，从0开始编号，新的编号j和原编号i的关系为：
通过对排序的坐标点进行循环位移方式调整排序后，使得随着编号的递增，代表管道形边界图像的四个边界点依次出现。
步骤600，对构成相界面的坐标点进行三次多项式拟合后得到三次函数，对两个壁面的坐标点进行一次函数拟合得到两条对应的表达式，求解后即可得到相界面两个端点的精确坐标，再将两者之间的所有坐标点由三次函数求出，拟合两个端点之间的所有坐标点，即得到相界面；将精确坐标分别代入三次函数的一次导数中，即得到相界面两端点处的接触角；
如图12所示，得到相界面两个端点精确坐标的步骤如下：
步骤610，设构成相界面的第一个坐标点为p，最后一个坐标点为q，用p至q范围内所有坐标点进行三次多项式拟合，以横坐标x为因变量，纵坐标y为自变量，拟合得到弯月面形状的三次函数x＝fpq(y)；
步骤620，设构成上壁面的第一个坐标点为d，最后一个坐标点为o；构成下壁面的第一个坐标点为r，最后一个坐标点为c；采用一次函数拟合，以横坐标y为因变量，纵坐标x为自变量，得到两条直线的表达式y＝kdox+cdo及y＝krcx+crc，其中kdo和krc为直线的斜率；
步骤630，根据拟合得到的三次函数及两条直线表达式，求解以下两个方程组；
即得到相界面上端点a的坐标(xa,ya)以及相界相下端点b的坐标(xb,yb)。如图13所示，至此相界面两个端点的坐标已知，相界面a-b范围内各点的坐标，计算的范围是对于y坐标为ya＜y＜yb，采用拟合得到的三次函数fpq进行计算，就可以得到a到b范围内所有点的坐标，就可以得到拟合的相界面a-b的曲线坐标。
获得接触角的过程如下：
根据相界面的拟合函数x＝fpq(y)，将上端点a的坐标(xa,ya)以及下端点b的坐标(xb,yb)分别代入fpq的一次导数中即可计算出接触角θa和θb：
将上端点a的坐标(xa,ya)带入fpq的一次导数中可以得到相界面的函数在上端点a的切线斜率，通过arctan函数可将斜率转化为角度，由于测量的是位于液体那一侧的接触角，因此需要加上π/2就可以得到θa；将下端点b的坐标(xb,yb)分别代入fpq的一次导数中，可以求得相界面在下端点b的切线斜率，同样采用arctan将斜率转化成角度，由于下端点b处的斜率是负值，因此转化成角度后需要乘上一个负号，再加上π/2就可以得到θb。
还包括对两条壁面的倾斜修正程序：
设两个壁面的倾斜角度为上壁面和下壁面倾斜角度的平均值，即(arctankdo+arctankrc)/2；则修正后的接触角θa和θb为：
步骤700，已知图像分辨率的情况下，计算出相界面两个端点到视野一侧的距离，而相界面到该侧的距离为两个端点到该距离的平均值，相界面某一时刻至微米毛细管通道入口的位移在前述基础上利用移动座的位移速度和距离的加和即可得到；相界面的速度是相界面位移的一阶导数，对位移以中心差分格式进行求导即得到相界面的当前速度。
显微图像的分辨率可以通过设置参数得到，用kpixel表示，即每个像素点对应的空间长度。前述的步骤已经获得了两个上端点(xa,ya)和下端点(xb,yb)的坐标，那么接触点与视野某侧边缘的实际距离可通过下述确定：
这里的侧边可以是摄像单元的视野左侧，也可以是视野右侧，本方案中以面向图像时的左侧作为测量对象，以下涉及到的视野侧边都是指左侧；相界面与视野左侧的距离d为两个接触点到左侧距离的平均，
而相界面某一时刻的真实位移，指的是相界面到微孔道入口的距离x(t)，这是一个随着时间变化的函数。令移动座从微米毛细管道入口开始追踪相界面，则x(t)是移动座的位移s和距离d的加和，根据式
其中，da为上端点至视野一侧的距离，db为下站点至视野同侧的距离，s为移动座的位移，d为相界面与视野同侧的距离，i＝1,2,3…m，δt为移动座的速度记录的时间间隔。相界面的位移测量结合了电机的控制和图像测量的结果，该结果为计算相界面的移动速度和接触角提供了基础。
相对相界面的位移x，相界面的速度v是位移x的一阶导数
由此得到相界面的速度。确定相界面的速度后，也可以反过来作为控制电机转动的条件。
步骤800，将相界面速度、相界面接触角、相界面位移代入计算微米毛细管通道内压力的两相流方程求出当前微米毛细管通道内压力，再根据流量和相界面速度之间的关系，即得到当前压力和流量之间的关系，进而实现通过流量控制压力的目的。
由于整个系统的压力损失在于微米毛经细通道内的两相流动过程，因此可以根据描述微米毛细管通道中两相流动的基本方程washburn公式推导得到流量q与目标压力pt之间的关系。根据washburn公式有：
式中v为相界面速度，x为相界面位移(即液体在微管道中的长度)，θ为相界面的接触角，μ为液体粘度，s为微米毛细管通道截面边长，a为微米毛细管通道截面面积，r为微米毛细管通道的等效半径，γ为界面张力。上式是针对矩形截面的微孔道，对于圆形截面的微孔道，系数0.5623需要更改为0.5。由于体积流量q与相界面速度v存在关系，q＝av，因此上式可改写为：
根据该式可以直接由目标压力pt、相界面接触角θ和相界面位移x直接计算出理想的注射流量q，得到流量q之后，可从控制程序将控制信号输出至注射泵控制器就可以实时控制注射泵的流量了，进而可以简化后续实验的计算量。
得到注射流量q后，可以此与电机的转速和移动座的移动速度建立联系，进而则可控制微流动芯片上相界面的位置。
本发明的优点在于鲁棒性强,测量结果稳定，具有较强容错性，对于质量不好的图像也能提供准确的测量结果；采用该计算方法编程简单，程序运行效率高，采用常规处理器处理一张2300×1730像素的图片耗时小于1s。
如图4、14所示，在本发明的一个实施方式中，为方便观察，该微流动芯片20可以是由透明材料制成的板状结构，在其一个平面上设置内凹的凹槽22，凹槽22的形状可以是圆形或是矩形，在凹槽22上盖有贴合凹槽22内上表面的活动盖板(同样为透明)，盖板可以采用卡扣或键槽一类结构固定在凹槽22上；微米毛细管通道21设置在凹槽22内的平面上，且可以阵列方式设置多个，各微米毛细管通道21相互之间平行，如本实施方式在微流动芯片20上设置六条微米毛细管通道21；每个微米毛细管通道21的两端分别设置有进液区211和出液区212，微米毛细管通道21的宽度可根据实验要求设置在10～100μm内，而进液区211和出液区212的宽度至少为微米毛细管通道21宽度的15～20倍；这是因为微米毛细管通道21本身的宽度是微米级，测试液体如果直接进入不易控制压力，采用进液区211和出液区212可以在微米毛细管通道21的两端形成一个测试液体进入和排出的缓冲通道，以引导测试液体的流动，使其在微米毛细管通道21内形成的相界面能够呈稳定的移动状态。
进液区211和出液区212的形状以使测试液体能够平稳进入微米毛细管通道21内为准，如圆形、椭圆；在本发明的一个实施方式中，该进液区211和出液区212采用水滴形，水滴形的尖端与微米毛细管通道21连通，尖端与微米毛细管通道21连接处形成的夹角可在60度左右，以使测试液体能够在收缩状态下平稳进入微米毛细管通道21。
在进液区211和出液区212内分别均匀地设置有多根垂直的防止微米塌陷的立柱213，立柱213与立柱213之间相互间隔，同时立柱213的直径可以和微米毛细管通道21的宽度相同或为其1～2倍，立柱213之间的间隔距离可以是立柱213的直径。采用立柱213是因为pdms材料属于弹塑性材料，在宽深比大于2(需要说明这个是什么意思)的孔道结构中，材料会因为重力的作用而产生微米量级的“塌陷”变形。立柱213可避免变形对测试液体流动造成的影响。而供测试液体进入和排出的进液孔和排液孔分别设置在进液区211和出液区212的中间位置且开口在凹槽22的底部。进液管路和排液管路可直接与相应的进液孔和排液孔连接。
立柱213的布置方式可以根据进液区211和出液区212的形状确定，如本实施方式中立柱213以环形的方式逐圈布置，中间留出进液孔和出液孔的设置空间。
如图3所示，在本发明的一个实施方式中，为方便移动和调整微流动芯片20，可以在移动座33的上表面安装两个相对间隔且垂直的支撑板36，支撑板36之间的间隔距离与微流动芯片20的宽度对应；在两个支撑板36的相对面上分别安装有固定或支撑微流动芯片20的固定台361；在移动座33上设置有多个安装孔331，安装孔331可以使支撑板36相对微流动芯片20的宽度调整安装位置；支撑板33通过l形折片35利用移动座33上的安装孔331与移动座33固定，l形折片35的两折边上分别设置有固定孔，其通过螺栓分别与移动座33和支撑板36固定。
进一步地，可以在微流动芯片20的下方两个支撑板36之间安装一个光源，光源用于使微流动芯片20更明亮，以方便摄像单元50获取清楚的相界面图像。具体的光源可以是led灯。
此外，为方便调整固定台361的高度，可以在每块支撑板36上分别设置垂直于放置后微流动芯片20的两道间隔滑动槽(图中未示出)，固定台361可利用穿过滑动槽的螺栓可调节的安装在支撑板36上。这里给出的是调整微流动芯片20的方法，在其它的实施方式中也可以直接调整数码相机的高度，而微流动芯片20保持不动。
此外，为方便观察，支撑板36、固定台361和l形折片35都可以采用与微流动芯片20同样的透明材料制作。
具体的固定台361可以为条状的支撑条，以对微流动芯片20形成水平支撑；还可以为设置有夹持微流动芯片20侧边的凹槽的夹持条，微流动芯片20在安装时，其相对的两侧可以夹持地安装在夹持条的凹槽内，该结构能够提高安装后微流动芯片20的稳定性。
如图15所示，在本发明的一个实施方式中，为方便安装压力传感器42，可以在移动座33的一侧通过l形折片35安装一个带三通的透明接头37，接头37其中两个通孔为水平且相对设置的直通371，另一个通孔为垂直于两个直通371中部的竖通372，在竖通372上通过螺纹安装有与控制系统10连接的压力传感器42，在两个直通371上通过螺纹连接有鲁尔接头373，注射泵41可通过柔性软管与鲁尔接头373连接，另一直通371的鲁尔接头373则通过柔性软管与进液区211的进液孔连接，而出液区212的出液孔通过柔性软管与接水盒43连接。
在测试过程中，则注射泵41推出的压力可以由压力传感器42直接测量并传送给控制系统10，控制系统10则根据该压力大小与预定的目标压力值进行比对，控制系统10在这个过程中根据压力传感器42的压力值持续控制注射泵41的注射过程，直至微米毛细管管路21中的压力值达到预定的目标压力值。
在本实施方式中，为防止接头处泄露，可仅在微米毛细管通道21中注入测试液体，保证测试液体的量能够在微米毛细管通道21中形成相应的相界面，后面则由注射泵41注入密封液体以在管路内形成注射压力，密封液体既可以是测试液体同一种液体，也可以是其它类液体，如水；在实验前将管路中预先充满液体，由液体传递压力并占据连接点处的缝隙。由于泄漏量与流体粘度成反比，空气粘度为0.01mpa·s，液体的粘度远远大于空气(如水的粘度在标况下约为1mpa·s)，采用液体作为压力传递介质可使泄漏量可以大大减少。在该结构下密封液体与测试液体之间有一段空间间隔(空气柱)，空气柱不接触涉及到管路与接头连接的区域，可起到隔离测试液体和密封液体的作用。
至此，本领域技术人员应认识到，虽然本文已详尽示出和描述了本发明的多个示例性实施例，但是，在不脱离本发明精神和范围的情况下，仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变型或修改。因此，本发明的范围应被理解和认定为覆盖了所有这些其他变型或修改。
- 还没有人留言评论。精彩留言会获得点赞！