CN107291993A - 一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 - Google Patents
一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 Download PDFInfo
- Publication number
- CN107291993A CN107291993A CN201710387669.4A CN201710387669A CN107291993A CN 107291993 A CN107291993 A CN 107291993A CN 201710387669 A CN201710387669 A CN 201710387669A CN 107291993 A CN107291993 A CN 107291993A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- crosslinked gel
- virtual boundary
- mark point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000000725 suspension Substances 0.000 title claims abstract description 21
- 239000007863 gel particle Substances 0.000 claims abstract description 29
- 239000002245 particle Substances 0.000 claims abstract description 17
- 239000012530 fluid Substances 0.000 claims abstract description 15
- 239000008187 granular material Substances 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000005452 bending Methods 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims abstract description 6
- 230000009471 action Effects 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000011438 discrete method Methods 0.000 claims description 2
- 239000011148 porous material Substances 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 10
- 230000007246 mechanism Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005325 percolation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G01V20/00—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/645—Fluid contacts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
本发明公开了一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,包括以下步骤:基于二维岩心CT切片图像,建立多孔介质模拟区域,进行数值网格离散;设定预交联凝胶颗粒的初始位置和原始形状,生成虚拟边界标记点;标记数值网格节点类型;根据预交联凝胶颗粒边界上数值网格节点的动量变化计算流体对各虚拟边界标记点的作用力;采用颗粒接触作用模型计算各虚拟边界标记点受到的接触作用力;根据预交联凝胶颗粒的当前形状和原始形状,计算各虚拟边界标记点受到的弯曲恢复力和拉伸恢复力;基于牛顿第二定律,计算各虚拟边界标记点的速度和位移,将各虚拟边界标记点相连即得到预交联凝胶颗粒的整体形状和位置。该方法能够处理复杂多孔介质、且流固耦合模拟精度高,可用于预交联凝胶等弹性颗粒悬浮液的微观渗流模拟研究。
Description
技术领域
本发明涉及一种多孔介质中预交联凝胶悬浮液微观流动的数值模拟方法,属于油气田开发提高原油采收率、油藏数值模拟和计算流体动力学的交叉领域。
背景技术
预交联凝胶颗粒驱是高含水油田进一步提高原油采收率的新兴深部调驱技术之一。然而,预交联凝胶悬浮液在多孔介质中的流动十分复杂,不仅具备刚性颗粒悬浮液的基本渗流特征,同时还具有弹性变形等特性。目前,考虑这些特性的油藏数值模拟技术多基于经典渗滤理论或尺寸排斥理论,属于宏观数值模拟的范畴,受限于连续介质假设,无法模拟预交联凝胶颗粒在多孔介质中“运移、封堵、变形、再运移”的真实过程。
为此,本发明提出一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,兼顾了颗粒变形的精细化处理、颗粒之间的接触作用模拟以及颗粒流体间的高效耦合,为深入认识预交联凝胶悬浮液的复杂渗流机理提供了一种研究手段。
发明内容
为弥补现有技术的不足,本发明将燃烧算法和转折点算法结合,提出一种可以识别孔喉末端的孔隙喉道识别方法,其识别过程能准确识别末端并且显著降低识别后的喉道长度与实际长度的相对误差。
本发明技术方案具体步骤如下:
(1)对图像进行二值化处理,区分岩石颗粒与孔喉空间;
(2)以岩石颗粒为前景像素对孔喉空间进行欧式长度变换,长度变换计算公式为:
其中p为孔喉空间内一点,q为p点周围任意的岩石颗粒点。
(3)由基于长度变换的转折点算法提取图像中线,此时的中线未考虑末端;
(4)基于转折点中线,进行不考虑末端时的孔喉识别;
(5)利用燃烧算法提取中线,查找燃烧中线尾端端点,确定末端所在区域;
(6)借助燃烧中线尾端端点,利用转折点算法提取末端中线;
(7)将步骤(4)中的转折点中线与步骤(6)中的末端转折点中线合并,对末端区域进行孔喉识别,然后将末端识别结果与步骤(4)中未考虑末端区域的孔喉识别结果合并,得到考虑末端的孔喉识别结果。
所述步骤(4)中,孔喉识别的步骤如下:
①删除位于孔隙内部的转折点中线路径,将剩余中线路径所对应像素按顺序编号,并以每个像素的长度变换结果作为此像素对应的孔隙半径;
②以像素编号为横坐标,对应的孔隙半径为纵坐标,做出每条路径的沿程曲线;
③将沿程曲线看作信号,选择离散Dmey波为小波基函数,对信号进行二尺度分解,选择合理阈值截取小波系数,然后重构信号完成小波消噪(参考文献:孙延奎.小波分析及其应用[M].DynoMedia Inc.,2005.),阈值选择公式为:
其中σ为噪声标准方差,N是信号长度;
④在去噪后的沿程曲线上查找最靠近两端的两个局部极小值,将其对应的像素作为孔隙和喉道的分界点即喉道的瓶颈,在瓶颈对应的位置画出孔喉识别线;
⑤孔喉识别线和岩石边界将非基质空间分成了若干封闭区域,将包含中线结点的区域作为孔隙,其余部分作为喉道。
所述步骤(6)中,提取末端中线步骤如下:
①从步骤(4)识别后的孔隙或喉道中,提取包含燃烧中线端点像素的喉道或孔隙区域即为末端区域;
②以燃烧中线的末端像素为原点,沿正交方向将提取的末端区域切分为四个子图,应用转折点算法分别提取四个子图的中线;
③将四个子图的中线合并在一起并清理分支,所得结果即为转折点中线在末端处的分支。
本发明具有以下有益效果及优点:
(1)以转折点算法为主要算法提取孔隙空间的中线,使中线的整体居中效果得到改进;
(2)利用小波消噪后的沿程曲线查找瓶颈,识别后的喉道长度与实际喉道长度相对误差显著降低;
(3)由燃烧算法获取的端点引导转折点算法提取末端部分的中线,可以将末端部分单独识别出来。
附图说明
图1为多孔介质中预交联凝胶悬浮液微观流动的模拟方法流程图。
图2为多孔介质模拟区域示意图。
图3为多孔介质中数值网格节点的类型标记示意图。
图4为多孔介质中预交联凝胶悬浮液微观流动的模拟结果图。
具体实施方式
图1给出了本发明的详细实施流程,下面将结合其他附图和具体实施例作进一步的说明,以便更好的理解本发明,但不限定本发明的范围。
(1)基于二维岩心CT切片图像,建立多孔介质模拟区域,如图2所示。该模拟区域多孔介质部分长度为260μm,为了模拟预交联凝胶颗粒进入和离开多孔介质的过程,分别在多孔介质入口端和出口端设置长度为40μm的管道,因此模拟区域总长度为340μm,宽度为100μm。设定网格步长为0.2μm,将该模拟区域进行数值网格离散,则共有1700×500=850000个计算节点。
(2)设定预交联凝胶颗粒的初始位置在多孔介质入口端,原始形状为圆形,直径为16μm。将颗粒的真实曲线边界离散生成18个虚拟边界标记点,并由其连线组成的18边形近似表征预交联凝胶颗粒。
(3)标记数值网格节点类型,若某数值网格节点及其相邻8个节点均被预交联凝胶颗粒覆盖则为颗粒内部节点(图3中实心圆点),若某数值网格节点被预交联凝胶颗粒覆盖但存在未被覆盖的相邻节点则为颗粒边界节点(图3中实心三角形),若某数值网格节点及其相邻8个节点均未被预交联凝胶颗粒覆盖则为流体内部节点(图3中空心圆点),若某数值网格节点未被预交联凝胶颗粒覆盖但存在被覆盖的相邻节点则为流体边界节点(图3中空心三角形)。图3中数值网格节点A为流体内部节点,B为颗粒内部节点,S为流体边界节点。
(4)输入网格步长Δx=0.2μm,时间步长Δt=0.0025s,松弛时间τ=0.65,根据预交联凝胶颗粒边界上数值网格节点的动量变化计算流体对各虚拟边界标记点的作用力,其计算公式为:
其中,Δx表示网格步长,Δt表示时间步长,Vt表示颗粒边界上以某数值网格节点为中心边长为Δx的控制单元的总体积,Vs表示上述控制单元中被颗粒覆盖部分的体积;τ表示松弛时间;Cj为附加碰撞项;lj表示j方向上的离散速度;p表示各虚拟边界标记点控制范围内属于颗粒边界类型的数值网格节点个数。
(5)计算虚拟边界标记点之间的接触重叠量δ,输入接触刚度k=104kg/s2,阻尼系数λ=0.3,计算各虚拟边界标记点受到的接触作用力,其计算公式为:
其中,k表示接触刚度,λ表示阻尼系数,δ表示各虚拟边界标记点之间的接触重叠量。
(6)输入预交联凝胶颗粒的弯曲恢复刚度γb=102kg/s2,拉伸恢复刚度γs=10kg/s2,根据预交联凝胶颗粒的当前形状和原始形状,计算各虚拟边界标记点受到的弯曲恢复力和拉伸恢复力,其计算公式分别为:
Fb=γb(α-α0)
其中,α0和α分别表示预交联凝胶颗粒变形前后相邻三个虚拟边界标记点之间的向量夹角,和分别表示预交联凝胶颗粒变形前后相邻两个虚拟边界标记点之间的距离,γb表示弯曲恢复刚度,γs表示拉伸恢复刚度。
(7)基于牛顿第二定律,计算各虚拟边界标记点的速度u和位移x,将各虚拟边界标记点相连即得到预交联凝胶颗粒的整体形状和位置;
x=xold+uΔt
其中,uold和u分别表示迭代前后虚拟边界标记点的速度大小,m表示预交联凝胶颗粒的质量,n表示虚拟边界标记点的个数,xold和x分别表示迭代前后虚拟边界标记点的位移大小。
(8)重复步骤(4)至步骤(7),当前后两次模拟结果相对误差小于0.01或者迭代次数大于100000时,结束模拟过程,输出模拟结果,如图4所示。
由上述实施例可以看出,本发明公开的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法能够直观反映预交联凝胶颗粒在多孔介质中的真实变形和运移过程,并能够定量表征预交联凝胶悬浮液在多孔介质中的微观渗流机制。因此,本发明为预交联凝胶悬浮液等含弹性颗粒悬浮液的微观渗流机制研究提供了一种有效的数值模拟手段。
Claims (7)
1.一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,包括以下步骤:
(1)基于二维岩心CT切片图像,建立多孔介质模拟区域,进行数值网格离散;
(2)设定预交联凝胶颗粒的初始位置和原始形状,生成虚拟边界标记点;
(3)标记数值网格节点类型,包括颗粒内部、颗粒边界、流体内部和流体边界四种类型;
(4)根据预交联凝胶颗粒边界上数值网格节点的动量变化计算流体对各虚拟边界标记点的作用力Ff;
(5)采用颗粒接触作用模型计算各虚拟边界标记点受到的接触作用力Fp;
(6)根据预交联凝胶颗粒的当前形状和原始形状,计算各虚拟边界标记点受到的弯曲恢复力Fb和拉伸恢复力Fs;
(7)基于牛顿第二定律,计算各虚拟边界标记点的速度u和位移x,将各虚拟边界标记点相连即得到预交联凝胶颗粒的整体形状和位置;
(8)重复步骤(4)至步骤(7),当前后两次模拟结果相对误差小于0.01或者迭代次数大于100000时,结束模拟过程,输出模拟结果。
2.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(2)中所述虚拟边界标记点为将颗粒的真实曲线边界离散后生成的虚拟边界标记点,可以离散为12、18或24个。
3.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(3)中所述标记数值网格节点类型的判断标准为:若某数值网格节点及其相邻8个节点均被预交联凝胶颗粒覆盖则为颗粒内部节点,若某数值网格节点被预交联凝胶颗粒覆盖但存在未被覆盖的相邻节点则为颗粒边界节点,若某数值网格节点及其相邻8个节点均未被预交联凝胶颗粒覆盖则为流体内部节点,若某数值网格节点未被预交联凝胶颗粒覆盖但存在被覆盖的相邻节点则为流体边界节点。
4.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(4)中所述流体对各虚拟边界标记点的作用力计算公式为:
<mrow>
<msub>
<mi>F</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<msub>
<mi>V</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<msub>
<mi>V</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>8</mn>
</munderover>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<msub>
<mi>l</mi>
<mi>j</mi>
</msub>
</mrow>
其中,Δx表示网格步长,Δt表示时间步长,Vt表示颗粒边界上以某数值网格节点为中心边长为Δx的控制单元的总体积,Vs表示上述控制单元中被预交联凝胶颗粒覆盖部分的体积;τ表示松弛时间;Cj为附加碰撞项;lj表示j方向上的离散速度;p表示各虚拟边界标记点控制范围内属于颗粒边界类型的数值网格节点个数。
5.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(5)中所述采用颗粒接触作用模型计算虚拟边界标记点接触作用力的计算公式为:
<mrow>
<msub>
<mi>F</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mi>k</mi>
<mi>&delta;</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&delta;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,k表示接触刚度,λ表示阻尼系数,δ表示各虚拟边界标记点之间的接触重叠量。
6.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(6)中所述各虚拟边界标记点受到的弯曲恢复力和拉伸恢复力计算公式分别为:
Fb=γb(α-α0)
其中,α0和α分别表示预交联凝胶颗粒变形前后相邻三个虚拟边界标记点之间的向量夹角,和分别表示预交联凝胶颗粒变形前后相邻两个虚拟边界标记点之间的距离,γb表示弯曲恢复刚度,γs表示拉伸恢复刚度。
7.根据权利要求1所述的一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法,其特征在于,步骤(7)中所述各虚拟边界标记点的速度和位移计算公式为:
<mrow>
<mi>u</mi>
<mo>=</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>o</mi>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>F</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mi>s</mi>
</msub>
</mrow>
<mrow>
<mi>m</mi>
<mo>/</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
x=xold+u Δt
其中,uold和u分别表示迭代前后虚拟边界标记点的速度大小,m表示预交联凝胶颗粒的质量,n表示虚拟边界标记点的个数,xold和x分别表示迭代前后虚拟边界标记点的位移大小。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710387669.4A CN107291993A (zh) | 2017-05-27 | 2017-05-27 | 一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 |
US15/989,861 US10371859B2 (en) | 2017-05-27 | 2018-05-25 | Method for simulation of microscopic flow of pre-crosslinked gel suspension liquid in porous medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710387669.4A CN107291993A (zh) | 2017-05-27 | 2017-05-27 | 一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107291993A true CN107291993A (zh) | 2017-10-24 |
Family
ID=60094859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710387669.4A Pending CN107291993A (zh) | 2017-05-27 | 2017-05-27 | 一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10371859B2 (zh) |
CN (1) | CN107291993A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109509220A (zh) * | 2018-11-06 | 2019-03-22 | 北京理工大学 | 一种模拟多孔介质固相转换器内部流体流动的方法 |
CN112749468A (zh) * | 2019-10-29 | 2021-05-04 | 中国石油化工股份有限公司 | 固相悬浮物在油气储集层孔隙中通过能力的数值模拟方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107085850B (zh) * | 2017-03-24 | 2019-10-01 | 大连成者科技有限公司 | 遮蔽采集图像中异物的标记体、识别图像中异物标记体的方法以及书籍扫描方法 |
CN109885894B (zh) * | 2019-01-24 | 2023-05-05 | 江西理工大学 | 一种基于Trefftz多源点配点法的稀土矿层地下水渗流模型 |
CN110263362B (zh) * | 2019-04-25 | 2022-11-29 | 南京大学(苏州)高新技术研究院 | 基于孔隙密度流的岩土体离散元流固耦合数值模拟方法 |
CN111737894B (zh) * | 2020-06-01 | 2023-04-07 | 电子科技大学 | 一种三角面网格划分有限元粒子模拟的粒子快速定位方法 |
CN111695281B (zh) * | 2020-06-01 | 2023-04-25 | 电子科技大学 | 一种四面体网格划分有限元粒子模拟的粒子快速定位方法 |
CN114757119B (zh) * | 2022-04-01 | 2023-08-08 | 河海大学 | 一种利用外包多边形构造管网的二维水力计算方法 |
CN116872499B (zh) * | 2023-08-03 | 2023-12-19 | 武汉必盈生物科技有限公司 | 一种可变层高的3d打印方法及系统 |
CN117390934B (zh) * | 2023-12-11 | 2024-03-29 | 武汉理工大学三亚科教创新园 | 多孔介质礁灰岩溶蚀过程的有限元模型构建方法及终端 |
CN117688823A (zh) * | 2024-02-04 | 2024-03-12 | 北京航空航天大学 | 一种岩土颗粒轨迹预测方法、电子设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129548A (zh) * | 2011-03-15 | 2011-07-20 | 中国石油大学(华东) | 多孔介质中流动泡沫结构图像实时采集与识别方法 |
CN103745082A (zh) * | 2013-12-11 | 2014-04-23 | 中国石油化工股份有限公司 | 非均相复合驱油体系数值模拟方法 |
CN104573215A (zh) * | 2014-12-31 | 2015-04-29 | 大港油田集团有限责任公司 | 一种调剖剂颗粒浓度及直径的计算方法 |
CN105043943A (zh) * | 2015-06-29 | 2015-11-11 | 中国石油大学(华东) | 一种动态测量预交联凝胶颗粒在多孔介质中粒径分布的实验装置及其测量方法 |
CN105158287A (zh) * | 2015-07-08 | 2015-12-16 | 上海大学 | 一种凝胶在岩石微通道中的运移性的可视化评价方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3140757B1 (en) * | 2014-05-05 | 2020-06-24 | Siemens Healthcare GmbH | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in models based on diseased patients and hypothetical normal anatomical models |
US10482215B2 (en) * | 2014-10-22 | 2019-11-19 | Indiana University Research And Technology Corporation | Unified computational method and system for patient-specific hemodynamics |
CN106598912B (zh) * | 2016-10-20 | 2023-09-01 | 浙江工业大学 | 一种基于cfd-dem耦合模型的磨粒流场分析方法 |
-
2017
- 2017-05-27 CN CN201710387669.4A patent/CN107291993A/zh active Pending
-
2018
- 2018-05-25 US US15/989,861 patent/US10371859B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129548A (zh) * | 2011-03-15 | 2011-07-20 | 中国石油大学(华东) | 多孔介质中流动泡沫结构图像实时采集与识别方法 |
CN103745082A (zh) * | 2013-12-11 | 2014-04-23 | 中国石油化工股份有限公司 | 非均相复合驱油体系数值模拟方法 |
CN104573215A (zh) * | 2014-12-31 | 2015-04-29 | 大港油田集团有限责任公司 | 一种调剖剂颗粒浓度及直径的计算方法 |
CN105043943A (zh) * | 2015-06-29 | 2015-11-11 | 中国石油大学(华东) | 一种动态测量预交联凝胶颗粒在多孔介质中粒径分布的实验装置及其测量方法 |
CN105158287A (zh) * | 2015-07-08 | 2015-12-16 | 上海大学 | 一种凝胶在岩石微通道中的运移性的可视化评价方法 |
Non-Patent Citations (1)
Title |
---|
KANG ZHOU等: "An efficient LBM-DEM simulation method for suspensions of deformable preformed particle gels", 《CHEMICAL ENGINEERING SCIENCE》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109509220A (zh) * | 2018-11-06 | 2019-03-22 | 北京理工大学 | 一种模拟多孔介质固相转换器内部流体流动的方法 |
CN109509220B (zh) * | 2018-11-06 | 2021-12-14 | 北京理工大学 | 一种模拟多孔介质固相转换器内部流体流动的方法 |
CN112749468A (zh) * | 2019-10-29 | 2021-05-04 | 中国石油化工股份有限公司 | 固相悬浮物在油气储集层孔隙中通过能力的数值模拟方法 |
CN112749468B (zh) * | 2019-10-29 | 2022-09-02 | 中国石油化工股份有限公司 | 固相悬浮物在油气储集层孔隙中通过能力的数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
US10371859B2 (en) | 2019-08-06 |
US20180341727A1 (en) | 2018-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107291993A (zh) | 一种多孔介质中预交联凝胶悬浮液微观流动的模拟方法 | |
Wang et al. | Training deeper convolutional networks with deep supervision | |
CN109543502B (zh) | 一种基于深度多尺度神经网络的语义分割方法 | |
Froyland et al. | Robust FEM-based extraction of finite-time coherent sets using scattered, sparse, and incomplete trajectories | |
CN109145983A (zh) | 一种基于轻量化网络的实时场景图像语义分割方法 | |
CN103400372B (zh) | 一种基于Reeb图描述的三维拓扑信息提取方法 | |
CN107679462A (zh) | 一种基于小波的深度多特征融合分类方法 | |
CN105096268A (zh) | 一种点云去噪平滑方法 | |
CN104834922A (zh) | 基于混合神经网络的手势识别方法 | |
CN105069415A (zh) | 车道线检测方法和装置 | |
CN109559315B (zh) | 一种基于多路径深度神经网络的水面分割方法 | |
Qu et al. | An improved YOLOv5 method for large objects detection with multi-scale feature cross-layer fusion network | |
CN105069434A (zh) | 一种视频中人体动作行为识别方法 | |
CN102867094A (zh) | 一种移动粒子半隐式算法中自由表面流动模型的构建方法 | |
CN111767810A (zh) | 一种基于D-LinkNet的遥感图像道路提取方法 | |
CN106022355A (zh) | 基于3dcnn的高光谱图像空谱联合分类方法 | |
CN103678483A (zh) | 基于自适应概率超图和半监督学习的视频语义分析方法 | |
CN104794501A (zh) | 模式识别方法及装置 | |
CN106295604A (zh) | 基于联合滤波的遥感影像路网提取技术 | |
CN104794455A (zh) | 一种东巴象形文字识别方法 | |
CN107730536A (zh) | 一种基于深度特征的高速相关滤波物体追踪方法 | |
CN106155540A (zh) | 电子毛笔笔形处理方法和装置 | |
Jiang et al. | Dfnet: Semantic segmentation on panoramic images with dynamic loss weights and residual fusion block | |
CN109544561A (zh) | 细胞标注方法、系统和装置 | |
CN109255181A (zh) | 一种基于多模型的障碍物分布仿真方法、装置以及终端 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171024 |