CN114136862B - 一种双润湿性页岩的液体表观渗透率计算方法 - Google Patents
一种双润湿性页岩的液体表观渗透率计算方法 Download PDFInfo
- Publication number
- CN114136862B CN114136862B CN202111436249.3A CN202111436249A CN114136862B CN 114136862 B CN114136862 B CN 114136862B CN 202111436249 A CN202111436249 A CN 202111436249A CN 114136862 B CN114136862 B CN 114136862B
- Authority
- CN
- China
- Prior art keywords
- shale
- wettability
- liquid
- contact angle
- double
- 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.)
- Active
Links
- 239000007788 liquid Substances 0.000 title claims abstract description 60
- 230000035699 permeability Effects 0.000 title claims abstract description 40
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 239000011148 porous material Substances 0.000 claims abstract description 46
- 239000005416 organic matter Substances 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 9
- 239000002071 nanotube Substances 0.000 claims description 8
- 230000009977 dual effect Effects 0.000 claims description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 6
- 229910052799 carbon Inorganic materials 0.000 claims description 6
- 238000001179 sorption measurement Methods 0.000 claims description 6
- 238000005481 NMR spectroscopy Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 5
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 5
- 229910052753 mercury Inorganic materials 0.000 claims description 5
- 238000000691 measurement method Methods 0.000 claims description 4
- 229910001873 dinitrogen Inorganic materials 0.000 claims 1
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000007723 transport mechanism Effects 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 abstract description 2
- 238000012512 characterization method Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 239000012530 fluid Substances 0.000 description 12
- 239000007789 gas Substances 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 229910052757 nitrogen Inorganic materials 0.000 description 4
- 230000002209 hydrophobic effect Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 229910019093 NaOCl Inorganic materials 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- SUKJFIGYRHOWBL-UHFFFAOYSA-N sodium hypochlorite Chemical compound [Na+].Cl[O-] SUKJFIGYRHOWBL-UHFFFAOYSA-N 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 108010046685 Rho Factor Proteins 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000013043 chemical agent Substances 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 229910010272 inorganic material Inorganic materials 0.000 description 1
- 239000011147 inorganic material Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012900 molecular simulation Methods 0.000 description 1
- 239000002090 nanochannel Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Biochemistry (AREA)
- Dispersion Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Immunology (AREA)
Abstract
本发明公开一种双润湿性页岩的液体表观渗透率计算方法,包括以下步骤:获取双润湿性页岩基质的和液体的相关参数;确定双润湿性页岩的有机质孔径分布、无机质孔径分布以及孔隙度;测量双润湿性页岩基质的表观接触角,并根据表观接触角计算有机质接触角和无机质接触角;最后计算双润湿性页岩的液体表观渗透率。本发明考虑了有机质与无机质的传输机理,并将速度分布、密度分布、表观粘度、滑移长度、页岩孔隙的迂曲度、孔隙度等参数引入页岩基质的宏观渗透率表征中;然后,将多孔介质的结构参数与单管模型相结合,耦合多种输运机制,建立了统一的液体表观渗透率模型,为页岩气的精确数值模拟提供了关键的输入参数,具有重要的理论与现实意义。
Description
技术领域
本发明涉及一种双润湿性页岩的液体表观渗透率计算方法,属于非常规油气开发技术领域。
背景技术
21世纪以来,随着经济的快速发展,世界能源消费的不断攀升,煤、石油、常规天然气等传统能源的消耗逐步增加,能源压力日益增大。在对能源巨大需求的新形势下,作为常规油气能源的有效补充,包括页岩气在内的非常规能源越来越受到重视。
页岩气储集层一般呈低孔、低渗透率的物性特征,开采难度较大。因此页岩气井需要实施储层压裂改造才能开采出来。在水力压裂作业过程中,数千立方米的压裂液被强行注入地下储层,而现场数据表明,在压裂液返排期间,仅部分注入流体可采出,大部分则不可采出。学者们提出了造成流体损失的原因,包括自吸、压裂液圈闭、粘土水化和蒸发扩散等。除上述因素外,页岩储层纳米级孔隙中的液体滑移使流体容易渗透到超致密的页岩基质中,可能是流体损失大的一个重要原因。因此,为了能准确预测出页岩气藏压裂施工过程压裂液的滤失量,计算页岩基质中液体的流动能力(表观渗透率)非常重要。
页岩孔隙可分为无机孔(亲水)和有机孔(疏水),因此,页岩是一种典型的双润湿性/混合湿性多孔介质。由于页岩的纳米孔隙直径与液体分子直径相当,这些孔隙中水的输运行为与常规储层有很大的不同。大量的实验和分子模拟表明,无论是疏水还是亲水的纳米通道,水在其中的运移都不遵守泊肃叶定律中经典的无滑移假设。液体结构易受到液体和壁面间的相互作用的强烈影响,从而改变了壁面边界条件。目前,针对气体的表观渗透率开展了大量的研究,但是针对液体的表观渗透率研究较少。为此,本发明基于单个纳米管中水的流动特性,结合无机和有机纳米孔不同的输运机制,引入弯曲度、孔隙度和总有机碳等结构参数,建立了非均质双润湿性页岩基质中液体流动的表观渗透率模型。
发明内容
为了克服现有技术中的问题,本发明提供一种双润湿性页岩的液体表观渗透率计算方法。
本发明解决上述技术问题所提供的技术方案是:一种双润湿性页岩的液体表观渗透率计算方法,包括以下步骤:
步骤S1、获取双润湿性页岩基质的和液体的相关参数;
步骤S2、确定双润湿性页岩的有机质孔径分布、无机质孔径分布以及孔隙度;
步骤S3、测量双润湿性页岩基质的表观接触角,并根据表观接触角计算有机质接触角和无机质接触角;
步骤S4、最后根据下式计算双润湿性页岩的液体表观渗透率;
ξin=(1-α)φ/τ
ξor=αφ/τ
式中:kapp为双润湿性页岩的液体表观渗透率;ρin为有机孔隙中的液体密度,ρor为无机孔隙中的液体密度;ρw为液体密度;ξin为有机质修正因子;ξor为无机质修正因子;α为总有机碳含量;τ为迂曲度;φ为页岩孔隙度;rin为有机质孔径;ror为无机质孔径;ls(in)为液体/有机质边界的滑移长度;ls(or)为液体/无机质边界的滑移长度。
进一步的技术方案是,所述步骤S2中根据压汞/氮气吸附实验确定双润湿性页岩的有机质孔径分布、无机质孔径分布。
进一步的技术方案是,所述步骤S2中根据核磁共振实验确定双润湿性页岩的孔隙度。
进一步的技术方案是,所述步骤S3中采用液滴角度测量法测量双润湿性页岩基质的表观接触角。
进一步的技术方案是,所述步骤S3中无机质接触角直接测量得到,再根据下式计算有机质接触角;
cos(θapp)=f1cos(θin)+f2cos(θor)
f1+f2=1
式中:f1、f2分别为有机质和无机质的表面积分数;θapp为表观接触角;θor为无机质接触角;θin为机质接触角。
进一步的技术方案是,所述步骤S4中的公式通过以下步骤建立:
对纳米孔隙中的流动行为进行表征,并计算纳米管中水的体积通量;
采用修正因子将纳米管的流动扩展到多孔介质中得到不同类型孔隙中的质量通量;
最后,将多孔介质的结构参数与单管模型相结合,得到双润湿性页岩的液体表观渗透率的计算公式。
本发明具有以下有益效果:本发明在计算页岩液体的表观渗透率时创新性的考虑了有机质与无机质的传输机理,针对亲水无机纳米孔模型考虑了液固吸附(密度与粘度空间分布),针对疏水有机质还考虑了速度滑移(近壁流动)等机理,并将速度分布、密度分布、表观粘度、滑移长度、页岩孔隙的迂曲度、孔隙度等参数引入页岩基质的宏观渗透率表征中。然后,将多孔介质的结构参数与单管模型相结合,耦合多种输运机制,建立了统一的液体表观渗透率模型,为页岩气的精确数值模拟提供了关键的输入参数,具有重要的理论与现实意义。
附图说明
图1为双润湿性页岩的液体表观渗透率计算流程图;
图2为有机质与无机质的孔径分布输入参数图;
图3为不同孔隙度、不同TOC含量对双润湿性页岩基质液体表观渗透率(ALP)的影响图;
图4为不同TOC含量下无机质接触角对双润湿性页岩基质ALP的影响图;
图5为不同有机质孔径下有机质接触角对双润湿性页岩基质ALP的影响图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的一种双润湿性页岩的液体表观渗透率计算方法,包括以下步骤:
步骤S1、采用压汞法获取目标岩心稍大尺度的孔径分布,采用低温氮气吸附法获取目标岩心稍小尺度的孔径分布,进而获取岩心的全尺度孔径分布A;
选取相邻位置所取得的岩心采用NaOCl等化学剂腐蚀岩心中的有机质,再次采用压汞法和氮气吸附法获取页岩的全尺度孔径分布B,该孔径分布校正有机质团簇的大小后即为无机质孔隙的孔径分布;通过孔径分布A减去无机质的孔径分布,则为有机质的孔径分布;
步骤S2、采用低场核磁共振仪对页岩进行孔隙度测试;通过核磁共振试验,获得横向弛豫时间与信号幅度的变化曲线;为了获得孔隙度,需要采用流体标样进行刻度,获得核磁信号量与孔隙度之间的线性拟合关系,从而将核磁信号量转化为孔隙度。
步骤S3、采用液滴角度测量法测量页岩的表观接触角;
其中将被测页岩表面磨成光面,将页岩放在烘箱中进行干燥,冷却一段时间后在页岩光面上滴一滴水,直径约1mm~2mm,然后通过光学***,将液滴放大、投影到屏幕上,拍照后上传到电脑中,经过放大处理后在图片上测出表观接触角;
再用测出的接触角计算有机质接触角和无机质接触角,一般无机质接触角2易于获取(选取无有机质的地方测量),有机质接触可由下式求得:
cos(θapp)=f1cos(θin)+f2cos(θor)
f1+f2=1
式中:f1、f2分别为有机质和无机质的表面积分数;θapp为表观接触角;θor为无机质接触角;θin为机质接触角;
步骤S4、对纳米孔隙中的流动行为进行表征,不可压缩流体在轴向压力驱动下通过稳定缓慢流过圆管(雷诺数远小于1)时的流量可表示为:
式中,μ∞为液体粘度;r为参考半径;是压力梯度;ls是液体/固体边界的滑移长度。其中,滑移长度ls描述了液体和固体之间的速度不连续,通常定义为切向速度分量消失时的外推长度:
式中:v为轴向速度;r0是管的半径。
对于特定的纳米管,滑移长度ls可以用原子显微镜或MD测量。对于光滑且无溶解气体的表面,滑移长度是接触角的函数,表示为:
式中:θ为接触角;c=0.41。
所以有:
步骤S5、计算流体的有效粘度;
除了液体滑移外,流体粘度的直径相关变化也可能影响流动。界面区粘度μi与体积水μ∞之间存在线性关系,可表示为:
空间粘度变化对流动行为的影响采用有效粘度表征,可以通过界面和体积区域粘度的加权平均得到:
式中:Ai和At分别为界面区域和整个管的横截面积。
步骤S6、界面区域的临界厚度可确定为0.7nm,约为两层分子厚度,其中水分子受到壁面分子的强烈影响,综上所述,可以得到纳米管中水的体积通量是:
步骤S7、引入页岩的两个基本结构参数(孔隙度和迂曲度),采用修正因子将无机和有机纳米管的流动模型扩展到多孔介质中,不同类型孔隙中的质量通量可表示为:
对于无机多孔介质,质量通量可表示为:
对于有机多孔介质,总质量通量为:
式中:A为孔的有效截面积,ξ为修正因子,ρin为有机孔隙中的液体密度,ρor为无机孔中的液体密度;
界面区域密度ρif与体相水密度ρb的比值与粘度比近似相等
修正因子ξ的表达式为:
ξ=φ/τ
无机孔隙和有机孔隙的孔隙度与储层固有孔隙度之间的关系如下:
φin=(1-α)φ
φor=αφ
式中:α为总有机碳含量;
因此,无机质修正因子和有机质修正因子为:
ξin=(1-α)φ/τ
ξor=αφ/τ
步骤S8、结合有机质(无机质)孔隙中的质量通量和修正因子,得到无机多孔介质和有机多孔介质的渗透率分别为:
最后得到页岩储层双润湿性纳米级多孔介质中液体表观渗透率:
步骤S9、最后根据上式计算双润湿性页岩的液体表观渗透率。
实施例
(1)收集页岩基质和水的相关参数以及其它基本参数,如表1所示;
(2)采用压汞法、低温氮气吸附法以及用NaOCl化学处理相结合的方法,获得了页岩的有机质和无机质孔径分布,无机质孔径取30纳米,无机质孔径取5纳米,如图2所示;
(3)通过核磁共振试验,获得核磁信号量与孔隙度之间的线性拟合关系,从而得到页岩孔隙度;
(4)采用液滴角度测量法直接测量接触角大小,再计算出有机质和无机质的接触角大小,无机质的接触角大小在0~90°范围内,有机质的接触角大小在90~180°范围内;
(5)将这些参数代入液体表观渗透率的计算公式,得到不同参数对液体表观渗透率的影响情况。
首先,保持除TOC以外的其他参数不变,以0.02为增幅,代入不同的TOC值,得到液体表观渗透率与TOC的变化关系,然后,改变孔隙度值,可以得到不同孔隙度页岩多孔介质的表观液体渗透率随TOC含量的变化情况,如图3所示;
同样,改变无机质接触角的值,从0~90°变化,以15°为间隔,得到液体表观渗透率与无机质接触角的变化关系,然后,改变TOC值,可以得到不同TOC含量页岩多孔介质的表观液体渗透率随无机质接触角的变化情况,如图4所示;
改变有机质接触角的大小,从90~160°变化,以10°为间隔,计算不同有机孔径下有机质接触角对双润湿性页岩基质ALP,如图5所示。
表1算例的模型输入参数
参数 | 符号 | 取值 |
孔隙度 | φ | 0.048 |
总有机碳含量 | α | 0.1 |
无机多孔介质接触角 | θin | 30° |
有机多孔介质接触角 | θor | 150° |
迂曲度 | τ | 2.540 |
近壁区厚度 | δ | 0.700 |
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (6)
1.一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,包括以下步骤:
步骤S1、获取双润湿性页岩基质的和液体的相关参数;
相关参数包括孔隙度、总有机碳含量、无机多孔介质接触角、有机多孔介质接触角、迂曲度、近壁区厚度;
步骤S2、确定双润湿性页岩的有机质孔径分布、无机质孔径分布以及孔隙度;
步骤S3、测量双润湿性页岩基质的表观接触角,并根据表观接触角计算有机质接触角和无机质接触角;
步骤S4、最后根据下式计算双润湿性页岩的液体表观渗透率;
ξin=(1-α)φ/τ
ξor=αφ/τ
式中:kapp为双润湿性页岩的液体表观渗透率;ρin为有机孔隙中的液体密度,ρor为无机孔隙中的液体密度;ρw为液体密度;ξin为有机质修正因子;ξor为无机质修正因子;α为总有机碳含量;τ为迂曲度;φ为页岩孔隙度;rin为有机质孔径;ror为无机质孔径;ls(in)为液体/有机质边界的滑移长度;ls(or)为液体/无机质边界的滑移长度。
2.根据权利要求1所述的一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,所述步骤S2中根据压汞/氮气吸附实验确定双润湿性页岩的有机质孔径分布、无机质孔径分布。
3.根据权利要求1所述的一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,所述步骤S2中根据核磁共振实验确定双润湿性页岩的孔隙度。
4.根据权利要求1所述的一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,所述步骤S3中采用液滴角度测量法测量双润湿性页岩基质的表观接触角。
5.根据权利要求4所述的一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,所述步骤S3中无机质接触角直接测量得到,再根据下式计算有机质接触角;
cos(θapp)=f1cos(θin)+f2cos(θor)
f1+f2=1
式中:f1、f2分别为有机质和无机质的表面积分数;θapp为表观接触角;θor为无机质接触角;θin为机质接触角。
6.根据权利要求1所述的一种双润湿性页岩的液体表观渗透率计算方法,其特征在于,所述步骤S4中的公式通过以下步骤建立:
对纳米孔隙中的流动行为进行表征,并计算纳米管中水的体积通量;
采用修正因子将纳米管的流动扩展到多孔介质中得到不同类型孔隙中的质量通量;
最后,将多孔介质的结构参数与单管模型相结合,得到双润湿性页岩的液体表观渗透率的计算公式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111436249.3A CN114136862B (zh) | 2021-11-29 | 2021-11-29 | 一种双润湿性页岩的液体表观渗透率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111436249.3A CN114136862B (zh) | 2021-11-29 | 2021-11-29 | 一种双润湿性页岩的液体表观渗透率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114136862A CN114136862A (zh) | 2022-03-04 |
CN114136862B true CN114136862B (zh) | 2024-04-09 |
Family
ID=80389989
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111436249.3A Active CN114136862B (zh) | 2021-11-29 | 2021-11-29 | 一种双润湿性页岩的液体表观渗透率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114136862B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108229048A (zh) * | 2018-01-17 | 2018-06-29 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑页岩基质差异的多场耦合渗透率计算方法 |
CN108710723A (zh) * | 2018-04-09 | 2018-10-26 | 西南石油大学 | 一种页岩储层多孔介质表观渗透率的计算方法 |
CN108959789A (zh) * | 2018-07-12 | 2018-12-07 | 西南石油大学 | 一种考虑应力敏感效应的页岩表观渗透率计算方法 |
CN109100278A (zh) * | 2018-07-18 | 2018-12-28 | 西南石油大学 | 一种考虑页岩孔径分布特征的表观渗透率计算方法 |
CN109632578A (zh) * | 2018-12-29 | 2019-04-16 | 西南石油大学 | 一种页岩储层强制自吸量预测方法 |
CN110210460A (zh) * | 2019-06-26 | 2019-09-06 | 中国石油大学(华东) | 一种考虑多重因素影响的页岩气表观渗透率计算方法 |
CA3042922A1 (en) * | 2018-05-15 | 2019-11-15 | Uti Limited Partnership | Method and system for measuring a parameter of a core sample |
CN111208052A (zh) * | 2020-02-28 | 2020-05-29 | 西安石油大学 | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 |
CN111291521A (zh) * | 2020-03-10 | 2020-06-16 | 西南石油大学 | 一种致密储层液相表观渗透率计算方法及装置 |
CN111401669A (zh) * | 2020-06-05 | 2020-07-10 | 西南石油大学 | 一种基于小波神经网络的页岩油压后返排率预测方法 |
CN111929219A (zh) * | 2020-08-12 | 2020-11-13 | 西南石油大学 | 一种页岩油藏油水两相相对渗透率计算方法 |
CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
CN112102487A (zh) * | 2020-09-17 | 2020-12-18 | 西南石油大学 | 一种基于多重混合分形的非常规储层三维渗透率确定方法 |
CN113051841A (zh) * | 2021-02-10 | 2021-06-29 | 中国石油大学(华东) | 一种强非均质页岩油储层表观渗透率确定方法和装置 |
CN113128040A (zh) * | 2021-04-12 | 2021-07-16 | 陕西延长石油(集团)有限责任公司 | 一种基于孔隙结构特征的页岩表观渗透率预测方法 |
-
2021
- 2021-11-29 CN CN202111436249.3A patent/CN114136862B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108229048A (zh) * | 2018-01-17 | 2018-06-29 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑页岩基质差异的多场耦合渗透率计算方法 |
CN108710723A (zh) * | 2018-04-09 | 2018-10-26 | 西南石油大学 | 一种页岩储层多孔介质表观渗透率的计算方法 |
CA3042922A1 (en) * | 2018-05-15 | 2019-11-15 | Uti Limited Partnership | Method and system for measuring a parameter of a core sample |
CN108959789A (zh) * | 2018-07-12 | 2018-12-07 | 西南石油大学 | 一种考虑应力敏感效应的页岩表观渗透率计算方法 |
CN109100278A (zh) * | 2018-07-18 | 2018-12-28 | 西南石油大学 | 一种考虑页岩孔径分布特征的表观渗透率计算方法 |
CN109632578A (zh) * | 2018-12-29 | 2019-04-16 | 西南石油大学 | 一种页岩储层强制自吸量预测方法 |
CN110210460A (zh) * | 2019-06-26 | 2019-09-06 | 中国石油大学(华东) | 一种考虑多重因素影响的页岩气表观渗透率计算方法 |
CN111208052A (zh) * | 2020-02-28 | 2020-05-29 | 西安石油大学 | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 |
CN111291521A (zh) * | 2020-03-10 | 2020-06-16 | 西南石油大学 | 一种致密储层液相表观渗透率计算方法及装置 |
CN111401669A (zh) * | 2020-06-05 | 2020-07-10 | 西南石油大学 | 一种基于小波神经网络的页岩油压后返排率预测方法 |
CN111929219A (zh) * | 2020-08-12 | 2020-11-13 | 西南石油大学 | 一种页岩油藏油水两相相对渗透率计算方法 |
CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
CN112102487A (zh) * | 2020-09-17 | 2020-12-18 | 西南石油大学 | 一种基于多重混合分形的非常规储层三维渗透率确定方法 |
CN113051841A (zh) * | 2021-02-10 | 2021-06-29 | 中国石油大学(华东) | 一种强非均质页岩油储层表观渗透率确定方法和装置 |
CN113128040A (zh) * | 2021-04-12 | 2021-07-16 | 陕西延长石油(集团)有限责任公司 | 一种基于孔隙结构特征的页岩表观渗透率预测方法 |
Non-Patent Citations (4)
Title |
---|
Apparent permeability model for shale oil with multiple mechanisms;Feng, Q., Xu, S., Wang, S., Li, Y., Gao, F., & Xu, Y;《Journal of Petroleum Science and Engineering》;20191231(第175期);814-827 * |
Numerical evaluation of shear and tensile stimulation volumes based on natural fracture failure mechanism in tight and shale reservoirs;Luo, Z., Zhang, N., Zhao, L., Ran, L., & Zhang, Y;《Environmental Earth Sciences》;20190228(第78期);1-15 * |
微纳米孔隙页岩气藏表观渗透率计算模型研究进展;刘化普;;《中外能源》;20170215(第02期);33-36 * |
考虑天然微裂缝的酸蚀蚓孔扩展数值模拟研究;赵立强;王润宇;刘平礼;梁冲;邹宏岚;罗志锋;;《油气藏评价与开发》;20200426(第02期);80-86 * |
Also Published As
Publication number | Publication date |
---|---|
CN114136862A (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100278B (zh) | 一种考虑页岩孔径分布特征的表观渗透率计算方法 | |
Cai et al. | Shale gas transport model in 3D fractal porous media with variable pore sizes | |
CN107461192B (zh) | 储层条件下页岩动态表观渗透率的计算方法 | |
CN105910971B (zh) | 富有机质致密岩心气体渗透率和扩散系数的联测方法 | |
CN108710723B (zh) | 一种页岩储层多孔介质表观渗透率的计算方法 | |
CN108959789B (zh) | 一种考虑应力敏感效应的页岩表观渗透率计算方法 | |
CN110598167B (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
CN101968423B (zh) | 低渗透储层启动压力测试方法 | |
Li et al. | Effect of water saturation on gas slippage in circular and angular pores | |
Lu et al. | Spontaneous imbibition in coal with in-situ dynamic micro-CT imaging | |
Gao et al. | Apparent permeability and gas flow behavior in carboniferous shale from the Qaidam Basin, China: an experimental study | |
CN108801879B (zh) | 一种页岩基质颗粒孔隙率和渗透率一体化测量***及方法 | |
Cao et al. | A new approach for measuring the permeability of shale featuring adsorption and ultra-low permeability | |
CN114136861B (zh) | 一种储气库近井地带干化盐析效应实验***及评价方法 | |
Wang et al. | Experimental study on permeability in tight porous media considering gas adsorption and slippage effect | |
CN106706492A (zh) | 一种研究多孔介质全区域渗流机制的渗透装置 | |
CN107525746A (zh) | 一种表征致密油藏低速非达西渗流特征的方法及装置 | |
CN109085104B (zh) | 一种致密气藏的气体相对渗透率的确定方法 | |
Dai et al. | Fractal model for permeability estimation in low-permeable porous media with variable pore sizes and unevenly adsorbed water lay | |
CN113128040A (zh) | 一种基于孔隙结构特征的页岩表观渗透率预测方法 | |
Chen et al. | Shale gas transport behavior considering dynamic changes in effective flow channels | |
CN114370269B (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN114136862B (zh) | 一种双润湿性页岩的液体表观渗透率计算方法 | |
Clerke et al. | Spontaneous imbibition of water into oil saturated M_1 bimodal limestone | |
CN113297779B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |