CN109829975A - 一种变孔隙度多孔介质构建方法 - Google Patents
一种变孔隙度多孔介质构建方法 Download PDFInfo
- Publication number
- CN109829975A CN109829975A CN201910193765.4A CN201910193765A CN109829975A CN 109829975 A CN109829975 A CN 109829975A CN 201910193765 A CN201910193765 A CN 201910193765A CN 109829975 A CN109829975 A CN 109829975A
- Authority
- CN
- China
- Prior art keywords
- porosity
- particle
- model
- radius
- target
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 20
- 239000002245 particle Substances 0.000 claims abstract description 91
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000012798 spherical particle Substances 0.000 claims abstract description 5
- 230000012010 growth Effects 0.000 claims description 5
- 230000034655 secondary growth Effects 0.000 claims 1
- 239000012530 fluid Substances 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 239000003344 environmental pollutant Substances 0.000 description 3
- 231100000719 pollutant Toxicity 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 230000029052 metamorphosis Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种变孔隙度多孔介质构建方法，包括在区域内随机分布多个球状颗粒，确定多孔介质目标模型的目标孔隙度，确定颗粒半径每次的增加长度，保持颗粒球心坐标固定，游历区域内的所有颗粒，在原颗粒半径的基础上，使每个颗粒的半径增加相同长度，直至达到目标孔隙度，输出多孔介质的目标模型。该方法还包括对多个目标孔隙度从打到小排序，在前一目标模型的基础上构建后一目标模型。该方法的构建过程可控，不同孔隙度多孔介质的通道特征之间的具有关联，并且可通过具体参数反映，构建过程与实际情况更接近，有助于多孔介质问题的研究。
Description
技术领域
本发明涉及多孔介质的构建，具体涉及一种变孔隙度多孔介质构建方法。
背景技术
多孔介质是一种典型的复杂三维结构，具有连接孔隙，流体可以很容易地穿透其中。多孔介质问题常见于诸多工业应用领域，例如低渗透和特地渗透油气田开发、地下水利用、煤气层的开采以及金属材料的制备等。在地下水污染问题研究中，含水层中多孔介质通道形态特征对于地下水流流动和污染物的迁移过程有重要影响，如优先流现象、不同渗透性含水层之间的污染物扩散。多孔介质的微观结构导致了宏观性质的差异，如孔隙度大小、孔隙通道形态特征、颗粒间距等不同导致污染物在其中的移动、扩散及反应过程的差异。对于这类问题的研究，不同孔隙通道特征的多孔介质模型的构建和对其特征的描述十分关键。
多孔介质内的微小空隙的总体积与该多孔介质的外表体积的比值被定义为多孔介质的孔隙度。孔隙度是影响多孔介质内流体容量和流体渗流状况的重要参量，孔隙度改变会导致孔隙通道形态变化，使得多孔介质中流体的流速分布以及溶质运移下浓度分布显著不同。变孔隙度多孔介质模型是开展此类研究的基础，准确构建符合真实情况的变孔隙度多孔介质模型十分关键。目前，对于不同孔隙度多孔介质模型的构建，往往直接采用两个独立模型。两个模型颗粒的数量、位置以及孔隙通道特征完全不同，导致了模型的结构差异较大影响了计算结果。
发明内容
发明目的：本发明目的是提供一种变孔隙多孔介质构建方法，通过控制初始多孔介质中颗粒半径的增长来改变多孔介质模型的孔隙度，进而生成新的多孔介质的模型。另外，在预设多个目标孔隙度的情况下，可以构建以下列孔隙度不同的多孔介质模型，且其孔隙通道呈现不断缩减的变化过程。
技术方案：一方面，本发明公开了一种变孔隙度多孔介质构建方法，该方法包括：(1)生成初始模型：在区域内随机分布多个球状颗粒，生成多孔介质初始模型；(2)设定目标参数：确定多孔介质目标模型的目标孔隙度；(3)根据颗粒分布情况和目标孔隙度确定颗粒半径每次的增加长度；(4)改变孔隙度：保持初始模型中的颗粒球心坐标固定，游历区域内的所有颗粒，在原颗粒半径的基础上，逐一使得颗粒的半径增加相同长度，改变模型孔隙度；(5)生成目标模型：若在游历过程中的模型孔隙度达到目标孔隙度的要求，则结束游历，生成并输出目标模型；若一轮游历结束后的模型孔隙度未达到目标孔隙度的要求，则重复步骤(3)(4)继续构建。
进一步地，步骤(4)还包括，在颗粒半径增长后对其与相邻颗粒的接触状态进行判断，若该颗粒与相邻颗粒接触，则该颗粒的半径恢复至此次增长前大小，且在之后的游历中，保持该颗粒半径不变。
进一步地，若相邻颗粒的球心距离小于颗粒半径之和，则判定此次增长会导致颗粒接触。
进一步地，目标参数包括K个目标孔隙度，K为大于1的正整数，步骤(2)还包括对K个目标孔隙度从大到小进行排序，依次作为第一至第K个目标模型的目标孔隙度。
进一步地，本发明的构建方法还包括，以前一目标模型作为初始模型，根据后一目标模型的目标孔隙度，重复步骤(3)(4)(5)继续构造后一目标模型。
进一步地，后一目标模型构建过程中颗粒半径每次的增加长度不超过前一目标模型构建过程中颗粒半径每次的增加长度。
另一方面，一种包括计算机可执行指令的计算机可读存储介质，计算机可执行指令在执行时，使计算机系统执行变孔隙度多孔介质构建方法，该方法包括：(1)生成初始模型：在区域内随机分布多个球状颗粒，生成多孔介质初始模型；(2)设定目标参数：确定多孔介质目标模型的目标孔隙度；(3)根据颗粒分布情况和目标孔隙度确定颗粒半径每次的增加长度；(4)改变孔隙度：保持初始模型中的颗粒球心坐标固定，游历区域内的所有颗粒，在原颗粒半径的基础上，逐一使得颗粒的半径增加相同长度，改变模型孔隙度；(5)生成目标模型：若在游历过程中的模型孔隙度达到目标孔隙度的要求，则结束游历，生成并输出目标模型；若一轮游历结束后的模型孔隙度未达到目标孔隙度的要求，则重复步骤(3)(4)继续构建。
作为对上述计算机可读存储介质的进一步限制，步骤(4)还包括，在颗粒半径增长后对其与相邻颗粒的接触状态进行判断，若该颗粒与相邻颗粒接触，则该颗粒的半径恢复至此次增长前大小，且在之后的游历中，保持该颗粒半径不变；其中，若相邻颗粒的球心距离小于颗粒半径之和，则判定此次增长会导致颗粒接触。
作为对上述计算机可读存储介质的进一步限制，目标参数包括K个目标孔隙度，K为大于1的正整数，步骤(2)还包括对K个目标孔隙度从大到小进行排序，依次作为第一至第K个目标模型的目标孔隙度。
作为对上述计算机可读存储介质的进一步限制，构建方法还包括，以前一目标模型作为初始模型，根据后一目标模型的目标孔隙度，重复步骤(3)(4)(5)继续构造后一目标模型。
有益效果：与现有技术相比，根据本发明的多孔介质构建方法，在保持颗粒位置不变的情况下增长颗粒半径达到孔隙度减小的目的。此外，本发明的构建方法的构建过程可控，不同孔隙度多孔介质的通道特征之间具有关联，并且可通过具体的参数反映。
附图说明
图1为变孔隙度多孔介质构建方法的流程图；
图2为实施例中的初始模型以及构建出的三个变孔隙度模型二维平面图；
图3为实施例中构建出的变孔隙多孔介质模型的局部孔隙通道。
具体实施方式
下面结合附图和实施例对本发明做进一步描述：
如图1所示，利用本发明的变孔隙多孔介质的构建方法构建目标孔隙度分别为0.50、0.45、0.40的多孔介质模型，步骤如下：
(1)生成初始模型：在区域内随机分布多个球形颗粒，生成一个多孔介质初始模型PM0。经测算，该初始模型PM0的初始孔隙度为n0＝0.55，初始平均半径为r0＝0.001m，本实施例中初始模型的颗粒数量为343个。
(2)设定目标参数：根据构建需求，确定目标模型的数量以及对应的孔隙度大小nk。本实施例中，目标模型的数量设置为3个，对应的孔隙度大小分别为n1＝0.50，n2＝0.45，n3＝0.40，将上述孔隙度按照从大到小的进行排序，分别对应第一目标模型、第二目标模型和第三目标模型。
(3)根据颗粒分布情况和目标孔隙度nk确定颗粒半径每次的增加长度dk。在本实施例中设定第一目标模型的构建过程中，颗粒半径每次的增加长度d1设置为d1＝r0*1/500；第二目标模型和第三目标模型构建过程中颗粒半径每次的增加长度d2、d3依次减小，本实施例中分别取d2＝r0*1/500、d3＝r0*1/2500。
(4)改变孔隙度：保持初始模型PM0中的343个颗粒球心坐标固定，游历区域内的每个颗粒，在原颗粒半径的基础上，使每个颗粒的半径增加相同长度d1，改变模型孔隙度；每个颗粒半径改变后，计算该颗粒与相邻颗粒之间的球心距离，若相邻颗粒的球心距离小于颗粒半径之和，则判定此次增长会导致相邻颗粒接触，则该颗粒半径恢复至此次增长之前大小，且在之后的游历中，保持该颗粒半径不变，以此实现模型孔隙度的改变。
(5)生成目标模型：每个颗粒后，判定当前模型孔隙度是否小于第一目标孔隙度，若达到要求，则生成并输出第一目标模型；若改变后的模型孔隙度未达到目标孔隙度的要求，则重复步骤(3)(4)继续构建。
(6)在第一目标模型的基础上，即以第一目标模型作为初始模型，重复重复步骤(3)(4)(5)继续构造第二目标模型；同理，在第二目标模型的基础上，通过同样的方法可构建出第三目标模型。
根据本发明的构建方法构建的初始模型PM0、第一目标模型PM1、第二目标模型PM2、第三目标模型PM3如图2、3所示，各目标模型的颗粒球心均固定为初始模型的颗粒球心，通过改变颗粒半径实现多孔介质模型孔隙度的改变，使本发明的孔隙度变化过程可控。
本发明中权利要求中方法步骤的序号不作为对构建方法执行顺序的限制。本发明公开的实施例可采用计算机可用或计算机可读介质上的计算机程序产品的形式，计算机可用或计算机可读存储介质具有具体实现在该介质中的供指令执行系统使用或结合指令执行系统使用的计算机可用或可读程序代码。
Claims (10)
1.一种变孔隙度多孔介质构建方法，其特征在于，包括以下步骤：
(1)生成初始模型：在区域内随机分布多个球状颗粒，生成多孔介质初始模型；
(2)设定目标参数：确定多孔介质目标模型的目标孔隙度；
(3)根据颗粒分布情况和目标孔隙度确定颗粒半径每次的增加长度；
(4)改变孔隙度：保持初始模型中的颗粒球心坐标固定，游历区域内的所有颗粒，在原颗粒半径的基础上，逐一使得颗粒的半径增加相同长度，改变模型孔隙度；
(5)生成目标模型：若在游历过程中的模型孔隙度达到目标孔隙度的要求，则结束游历，生成并输出目标模型；若一轮游历结束后的模型孔隙度未达到目标孔隙度的要求，则重复步骤(3)(4)继续构建。
2.根据权利要求1所述的方法，其特征在于，所述步骤(4)还包括，在颗粒半径增长后对其与相邻颗粒的接触状态进行判断，若该颗粒与相邻颗粒接触，则该颗粒的半径恢复至此次增长前大小，且在之后的游历中，保持该颗粒半径不变。
3.根据权利要求2所述的方法，其特征在于，若相邻颗粒的球心距离小于颗粒半径之和，则判定此次增长会导致颗粒接触。
4.根据权利要求1所述的方法，其特征在于，所述目标参数包括K个目标孔隙度，K为大于1的正整数，所述步骤(2)还包括对K个目标孔隙度从大到小进行排序，依次作为第一至第K个目标模型的目标孔隙度。
5.根据权利要求4所述的方法，其特征在于，所述方法还包括，以前一目标模型作为初始模型，根据后一目标模型的目标孔隙度，重复步骤(3)(4)(5)继续构造后一目标模型。
6.根据权利要求5所述的方法，其特征在于，后一目标模型构建过程中颗粒半径每次的增加长度不超过前一目标模型构建过程中颗粒半径每次的增加长度。
7.一种包括计算机可执行指令的计算机可读存储介质，所述计算机可执行指令在执行时，使计算机系统执行变孔隙度多孔介质构建方法，所述方法包括：
(1)生成初始模型：在区域内随机分布多个球状颗粒，生成多孔介质初始模型；
(2)设定目标参数：确定多孔介质目标模型的目标孔隙度；
(3)根据颗粒分布情况和目标孔隙度确定颗粒半径每次的增加长度；
(4)改变孔隙度：保持初始模型中的颗粒球心坐标固定，游历区域内的所有颗粒，在原颗粒半径的基础上，逐一使得颗粒的半径增加相同长度，改变模型孔隙度；
(5)生成目标模型：若在游历过程中的模型孔隙度达到目标孔隙度的要求，则结束游历，生成并输出目标模型；若一轮游历结束后的模型孔隙度未达到目标孔隙度的要求，则重复步骤(3)(4)继续构建。
8.根据权利要求7所述的计算机可读存储介质，其特征在于，所述步骤(4)还包括，在颗粒半径增长后对其与相邻颗粒的接触状态进行判断，若该颗粒与相邻颗粒接触，则该颗粒的半径恢复至此次增长前大小，且在之后的游历中，保持该颗粒半径不变；其中，若相邻颗粒的球心距离小于颗粒半径之和，则判定此次增长会导致颗粒接触。
9.根据权利要求8所述的计算机可读存储介质，其特征在于，所述目标参数包括K个目标孔隙度，K为大于1的正整数，所述步骤(2)还包括对K个目标孔隙度从大到小进行排序，依次作为第一至第K个目标模型的目标孔隙度。
10.根据权利要求9所述的计算机可读存储介质，其特征在于，所述方法还包括，以前一目标模型作为初始模型，根据后一目标模型的目标孔隙度，重复步骤(3)(4)(5)继续构造后一目标模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910193765.4A CN109829975A (zh) | 2019-03-14 | 2019-03-14 | 一种变孔隙度多孔介质构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910193765.4A CN109829975A (zh) | 2019-03-14 | 2019-03-14 | 一种变孔隙度多孔介质构建方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109829975A true CN109829975A (zh) | 2019-05-31 |
Family
ID=66869017
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910193765.4A Pending CN109829975A (zh) | 2019-03-14 | 2019-03-14 | 一种变孔隙度多孔介质构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109829975A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427732A (zh) * | 2019-09-11 | 2019-11-08 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
| CN112668228A (zh) * | 2021-01-05 | 2021-04-16 | 中国石油天然气股份有限公司 | 二维多孔介质的模型构建方法、装置、设备和存储介质 |
| CN112668254A (zh) * | 2020-12-28 | 2021-04-16 | 大连理工大学 | 一种考虑颗粒大小分布和孔隙度的致密数字岩心生成方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105677963A (zh) * | 2016-01-05 | 2016-06-15 | 中国石油天然气股份有限公司 | 一种构建多孔介质模型的方法、服务器及系统 |
| CN106227943A (zh) * | 2016-07-26 | 2016-12-14 | 中国石油大学(华东) | 一种二维不同表面粗糙度的随机孔隙生成方法 |
| CN107423498A (zh) * | 2017-07-13 | 2017-12-01 | 山东大学 | 一种高致密度离散颗粒多相体系的建模方法 |
| CN107784169A (zh) * | 2017-10-11 | 2018-03-09 | 上海电力学院 | 基于变差函数及孔隙度的多孔介质表征单元体确定方法 |
| CN108734751A (zh) * | 2018-05-08 | 2018-11-02 | 中山大学 | 一种基于随机分形的多孔介质细观结构模型生成方法 |
-
2019
- 2019-03-14 CN CN201910193765.4A patent/CN109829975A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105677963A (zh) * | 2016-01-05 | 2016-06-15 | 中国石油天然气股份有限公司 | 一种构建多孔介质模型的方法、服务器及系统 |
| CN106227943A (zh) * | 2016-07-26 | 2016-12-14 | 中国石油大学(华东) | 一种二维不同表面粗糙度的随机孔隙生成方法 |
| CN107423498A (zh) * | 2017-07-13 | 2017-12-01 | 山东大学 | 一种高致密度离散颗粒多相体系的建模方法 |
| CN107784169A (zh) * | 2017-10-11 | 2018-03-09 | 上海电力学院 | 基于变差函数及孔隙度的多孔介质表征单元体确定方法 |
| CN108734751A (zh) * | 2018-05-08 | 2018-11-02 | 中山大学 | 一种基于随机分形的多孔介质细观结构模型生成方法 |
Non-Patent Citations (2)
| Title |
|---|
| PEIPEI YANG 等: "Permeability in multi-sized structures of random packed porous media using three-dimensional lattice Boltzmann method", 《INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER》 * |
| 乔太斌: "多孔介质中固体体积分数与颗粒尺度对流体绝对渗透率的影响", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427732A (zh) * | 2019-09-11 | 2019-11-08 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
| CN110427732B (zh) * | 2019-09-11 | 2023-04-18 | 河海大学 | 一种双层非均质土层结构模型构建方法 |
| CN112668254A (zh) * | 2020-12-28 | 2021-04-16 | 大连理工大学 | 一种考虑颗粒大小分布和孔隙度的致密数字岩心生成方法 |
| CN112668228A (zh) * | 2021-01-05 | 2021-04-16 | 中国石油天然气股份有限公司 | 二维多孔介质的模型构建方法、装置、设备和存储介质 |
| CN112668228B (zh) * | 2021-01-05 | 2022-11-04 | 中国石油天然气股份有限公司 | 二维多孔介质的模型构建方法、装置、设备和存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109829975A (zh) | 一种变孔隙度多孔介质构建方法 | |
| CN110555523B (zh) | 一种基于脉冲神经网络的短程跟踪方法及系统 | |
| CN111308896B (zh) | 基于可变误差的非线性系统自适应最优控制方法 | |
| CN107703751A (zh) | 基于蜻蜓算法的pid控制器优化方法 | |
| CN105260607A (zh) | 一种串并联耦合的多模型水文预报方法 | |
| Lin et al. | Traffic signal optimization based on fuzzy control and differential evolution algorithm | |
| CN113744527B (zh) | 一种面向高速公路合流区的智能靶向疏堵方法 | |
| CN105045091A (zh) | 基于模糊神经控制系统的疏浚工艺智能决策分析方法 | |
| CN105808885B (zh) | 一种基于遗传算法的日照约束下建筑容积率计算方法 | |
| CN113569438A (zh) | 基于多源降雨融合和实时校正的城市洪涝模型构建方法 | |
| Zhao et al. | Traffic signal control with deep reinforcement learning | |
| CN113468466A (zh) | 基于神经网络的多工况一维波动方程求解方法 | |
| CN109359320A (zh) | 基于多采样率自回归分布滞后模型的高炉指标预测方法 | |
| CN109492768A (zh) | 空气质量预测模型训练方法及装置 | |
| Chen et al. | An optimal algorithm design of rssi indoor location based on neural network | |
| CN110555225B (zh) | 一种基于分层粒子群优化算法的rbpf-slam计算方法 | |
| CN117172116A (zh) | 基于pso-bp组合模型的混凝土抗压强度预测方法 | |
| CN111445076A (zh) | 一种基于lstm的空气污染分析方法 | |
| CN105631169A (zh) | 一种索杆结构初始应变确定方法 | |
| CN115993843A (zh) | 一种应用于群体智能系统的最优编队控制方法 | |
| CN104616507A (zh) | 一种交通子区信号周期的协调控制方法及系统 | |
| CN116256468A (zh) | 一种基于动态自适应粒子群的空气质量预测优化方法 | |
| CN114896909A (zh) | 一种基于水位高度的明渠流量计算方法 | |
| CN109063306B (zh) | 一种网格化河北模型的土壤下渗能力空间离散方法 | |
| Yu et al. | Minimize pressure difference traffic signal control based on deep reinforcement learning |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190531 |