CN111208052A - 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 - Google Patents
基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 Download PDFInfo
- Publication number
- CN111208052A CN111208052A CN202010130359.6A CN202010130359A CN111208052A CN 111208052 A CN111208052 A CN 111208052A CN 202010130359 A CN202010130359 A CN 202010130359A CN 111208052 A CN111208052 A CN 111208052A
- Authority
- CN
- China
- Prior art keywords
- pore
- pores
- shale
- kozeny
- shale reservoir
- 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.)
- Granted
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011148 porous material Substances 0.000 claims abstract description 233
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 238000012512 characterization method Methods 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims description 13
- 238000009792 diffusion process Methods 0.000 claims description 13
- 238000001179 sorption measurement Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 11
- 230000035945 sensitivity Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 25
- 238000004364 calculation method Methods 0.000 description 9
- 230000000704 physical effect Effects 0.000 description 6
- 230000002195 synergetic effect Effects 0.000 description 4
- 238000003795 desorption Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Porous Artificial Stone Or Porous Ceramic Products (AREA)
Abstract
本发明公开了基于改进Kozeny‑Carman模型的页岩储层渗透率预测方法,将页岩储层孔隙截面为非圆形孔隙校正为圆形孔隙;利用孔隙等效直径将孔隙尺寸演化进行定量表征;基于分形粗化理论确定孔隙尺寸分形因子;根据孔隙尺寸分形因子,计算毛管长度以及预设范围的孔隙数,利用实时孔隙等效半径和毛管长度计算孔隙体积,利用毛管长度以及预设范围的孔隙数计算基质体积,利用孔隙体积和基质体积计算动态孔隙度,利用实时孔隙等效半径计算不同截面的单根毛管表观渗透率,根据不同截面孔隙占比以及不同截面的单根毛管表观渗透率计算页岩表观渗透率。本发明综合考虑页岩储层特性,基于分形理论动态表征页岩储层孔隙度,并将此动态孔隙度应用于页岩储层渗透率定量表征。
Description
技术领域
本发明所属技术领域为非常规油气开发中页岩储层渗透率的理论计算领域,具体涉及基于改进Kozeny-Carman模型的页岩储层渗透率预测方法。
背景技术
页岩气作为一种新型的非常规天然气资源,以其分布范围广、资源量大等特点正日益受到关注和重视。如何科学高效地开发这类非常规能源,需要从页岩储层特性出发,针对其截然不同的孔隙结构、储集方式等特征,建立相应的定量表征储层物性参数的计算方法或模型,孔隙度和渗透率作为储层物性的关键基础参数,是储量计算及气藏动态预测的基础,准确定量表征相关物性参数,有助于合理配产及制定开发方案,对页岩气高效大规模开发起着至关重要的意义。
经典的Kozeny-Carman方程,被广泛应用于常规储层的渗透率预测,该方式基于孔隙度将单根毛管Posenille流动进行粗化,从而预测岩芯尺度渗透率,通过毛管束模型计算的孔隙度过于简化孔隙截面几何形态及孔隙尺寸等微观非均质性,其假设所有毛管为等径圆形毛管,因此将单根毛管粗化至毛管束过程中,因过于简化孔隙结构的复杂性及非连续性,在一定程度上造成粗化后的渗透率理论值与实测数据的偏差。同时页岩基质内复杂的多尺度孔隙结构、有机孔隙及应力敏感的狭缝形孔隙大量发育,这一系列的储层特征使得页岩储层微观非均质性通常强于常规储层。基于等径圆形毛管建立的经典Kozeny-Carman方程虽然被被广泛应用于常规储层的渗透率预测,但由于页岩储层孔隙结构复杂,有机孔、无机孔及应力敏感狭缝形孔隙共同发育,微观非均质性较强,经典Kozeny-Carman模型用于表征页岩储层相关物性存在如下问题:1、未考虑页岩储层孔隙截面几何形态的多样性;2、未考虑页岩储层跨尺度孔隙结构特征;3、未考虑页岩储层泄压动态开发过程中,由于吸附及应力敏感引起的孔隙尺寸演化;4、未考虑滑脱边界、表面及Knudsen扩散等传质行为。因此应用传统Kozeny-Carman方程理论计算得到的页岩储层物性参数与实测数据偏差较大,从而导致产能预测结果出现较大偏差。因此需要对一种新的方法对页岩渗透率进行有效预测。
发明内容
为解决现有技术中存在的问题,本发明的目的在于提供基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,本发明综合考虑页岩储层特性,基于分形理论动态表征页岩储层孔隙度,并将此动态孔隙度应用于页岩储层渗透率定量表征。
本发明采用的技术方案如下:
基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,包括如下步骤:
S1,针对页岩储层孔隙截面几何形态多样性进行归一化校正,将页岩储层孔隙截面为非圆形孔隙校正为圆形孔隙,得到孔隙等效直径;
S2,针对页岩气藏泄压动态开发过程中的孔隙尺寸演化,利用孔隙等效直径将孔隙尺寸演化进行定量表征,得到实时孔隙等效半径;
S3,基于分形粗化理论确定孔隙尺寸分形因子;根据孔隙尺寸分形因子,计算毛管长度以及预设范围的孔隙数,利用实时孔隙等效半径和毛管长度计算孔隙体积,利用毛管长度以及预设范围的孔隙数计算基质体积,利用孔隙体积和基质体积计算动态孔隙度,利用实时孔隙等效半径计算不同截面的单根毛管表观渗透率,根据不同截面孔隙占比以及不同截面的单根毛管表观渗透率计算单根毛管平均表观渗透率,根据单根毛管平均表观渗透率和动态孔隙度计算得到改进后的Kozeny-Carman方程,从而计算得到页岩表观渗透率对页岩储层渗透率进行预测。
优选的,S1中,通过校正因子γ将截面为非圆形的孔隙校正为圆形孔隙,校正后的孔隙等效直径Dequ为:
Dequ=γD
上式中:
γ为校正因子,孔隙截面为圆形时γ=1,孔隙截面为正方形时γ=1.094,孔隙截面为等边三角形时γ=1.186;
D为孔隙特征长度,单位为m。
优选的,S2中,实时孔隙等效半径rp为:
rp=rin+drdis+drpm-ra=f(γD,pp)
其中,drdis为基质收缩引起的孔隙半径变化,单位为m;drpm为应力敏感引起的孔隙半径变化,单位为m;ra为吸附层厚度,单位为m。
优选的,基质收缩引起的孔隙半径变化drdis为:
其中,rin为初始孔隙等效半径(Dequ/2),单位为m;Φint为初始孔隙度;εL为Langmuir应变;PL为Langmuir压力,单位为Pa;Pin为初始孔隙压力,单位为Pa;Pp为实时孔隙压力,单位为Pa。
优选的,应力敏感引起的孔隙半径变化drpm为:
其中,αf为页岩孔隙的Biot系数;αm为页岩基质的Biot系数;Kn为页岩孔隙平均方向刚度,单位为Pa/m;Em为页岩基质杨氏模量,单位为Pa;s为页岩基质特征长度,单位为m;Pin为初始孔隙压力,单位为Pa;Pp为实时孔隙压力,单位为Pa。
优选的,吸附层厚度ra为:
ra=ppdm/(pL+pp)
其中,dm为气体分子直径,单位为m;PL为Langmuir压力,单位为Pa;Pp为实时孔隙压力,单位为Pa。
优选的,S3中,基于分形粗化理论,孔隙尺寸大于或等于D的累积孔隙数N为:
其中:Dmin为孔隙特征长度下限,单位为m;Dmax为孔隙特征长度上限,单位为m;Df为孔隙尺寸分形因子,
Df为:
基于分形理论,计算在[D,D+dD]之间的孔隙数-dN:
毛管长度L0为:
优选的,S3中,孔隙体积Vpi为:
其中,ωi为不同截面孔隙占比;Vpi为不同截面孔隙对应的孔隙体积,单位为m3;Dmin为孔隙特征长度下限;Dmax为孔隙特征长度上限;ra为吸附层厚度,rp为实时孔隙等效半径,单位为m;L0为毛管长度,单位为m;
S3中,基质体积Vm为:
其中,ωi为不同截面孔隙占比;Vmi为不同截面孔隙对应的基质体积,单位为m3;Dmin为孔隙特征长度下限,单位为m;Dmax为孔隙特征长度上限,单位为m;L0为毛管长度,单位为m;rin为初始孔隙等效半径,大小为孔隙等效直径的一半,单位为m;Φint为初始孔隙度;εL为Langmuir应变;PL为Langmuir压力,单位为Pa;Pin为初始孔隙压力,单位为Pa;Pp为实时孔隙压力,单位为Pa;i代表不同截面孔隙的种类,即当孔隙截面为圆形、正方形或等边三角形时各自应一个i;λ为不同截面孔隙的种类数,当只有圆形截面的孔隙、正方形截面的孔隙或等边三角形截面的孔隙时λ为1,当有圆形截面的孔隙、正方形截面的孔隙和等边三角形截面的孔隙中的任意两种时λ为2,当同时具有圆形截面的孔隙、正方形截面的孔隙和等边三角形截面的空隙时λ为3。
优选的,S3中,动态孔隙度Φp为:
其中,Vm为基质体积,单位为m3;Vpi为孔隙体积,单位为m3;
考虑不同截面单根毛管的平均表观渗透率Ka为:
其中,ωi为不同截面孔隙占比;μ为气体粘度,单位为pa.s;ρavg为气体平均密度,单位为kg/m;R为气体常数;α为切向动量调节系数;Pavg为毛管中的平均压力;T为温度,单位为T;Ds为气体表面扩散系数,单位为m2/s;Cs为吸附气体浓度,单位为kg/m3;M为气体分子量,单位为kg/mol;rp为实时孔隙等效半径,单位为m;ζms为表面扩散修正因子;i代表不同截面孔隙的种类;λ为不同截面孔隙的种类数。
优选的,S3中,基于改进的Kozeny-Carman方程计算页岩表观渗透率K为:
K=ΦpKa。
本发明具有如下有益效果:
本发明基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,基于经典的Kozeny-Carman方程,就其对页岩储层进行定量表征的局限性进行修正,本发明的特征主要包括如下三方面:1、考虑页岩储层孔隙截面几何形态的多样性,将页岩储层孔隙截面为非圆形孔隙校正为圆形孔隙;2、考虑页岩储层跨尺度孔隙结构特征,基于分形理论将微尺度非均质性融入模型粗化过程;3、综合考虑页岩储层泄压动态开发过程中,孔隙尺寸演化造成的渗流通道的实时改变引起的储层物性基础参数孔隙度及渗透率的实时改变。本发明的方法客观还原页岩储层特征,建立的相关模型不仅是孔隙几何形态、尺寸的函数,同时也是压力的函数,有效弥补现有理论计算方法的不足,提高了页岩储层孔隙度及渗透率预测的准确性,利于掌握页岩气藏动态开发规律,为页岩气藏高效开发提供理论依据。
具体实施方式
下面结合实施例来对本发明做进一步的说明。
针对经典的Kozeny-Carman方程表征页岩储层渗透率的局限性,基于页岩储层特征,对经典的Kozeny-Carman方程进行改进,改进后的模型考虑页岩储层泄压动态开发过程中,因吸附气解吸造成的孔隙壁面吸附层厚度变化及页岩基质收缩协同影响下的渗流通道的改变,同时耦合储层应力敏感造成的孔隙渗流通道的改变,综合考虑以上三重机理协同作用下的孔隙尺寸动态演化并结合页岩岩芯静态参数(孔隙尺寸分布及孔隙截面几何形态),基于分形理论求出多机理协同影响下的动态孔隙度,然后基于滑脱边界、表面扩散及Knudsen扩散建立的单根毛管传质模型,最终求得页岩渗透率。
本发明基于页岩储层特征,从如下四个方面对经典的Kozeny-Carman方程进行修正:1、考虑页岩储层孔隙截面几何形态的多样性;2、考虑页岩储层跨尺度孔隙结构特征;3、考虑页岩储层泄压动态开发过程中,由于吸附及应力敏感引起的孔隙尺寸演化;4、考虑基于滑脱边界、表面扩散及Knudsen扩散建立的传质模型替换传统的Posenille流动模型。
本发明所述基于改进Kozeny-Carman模型的页岩储层渗透率预测方法主要包括三部分:1、针对页岩储层孔隙截面几何形态多样性进行归一化校正,提出校正因子,将截面为非圆形孔隙校正为圆形孔隙;2、针对页岩气藏泄压动态开发过程中,孔隙尺寸演化进行定量表征;3、针对页岩储层多尺度孔隙特征,结合分形理论考虑页岩储层孔隙尺度的非均质性,将PSD数据融入孔隙体积及基质体积的求解过程中。最后将不同几何截面孔隙占比的权重分别乘以该类孔隙计算出的孔隙体积及基质体积然后分别进行求和得到总的孔隙体积及总的基质体积,最终得到页岩样品的动态孔隙度理论计算值。基于该动态孔隙度并结合考虑滑脱边界、表面扩散及Knudsen扩散建立单根毛管传质模型,最终求得页岩渗透率。本发明基于改进Kozeny-Carman模型的页岩储层渗透率预测方法的具体步骤如下:
(1)通过校正因子γ将截面为非圆形的孔隙校正为圆形孔隙Dequ:
Dequ=γD
上式中:Dequ:孔隙等效直径,单位为m;
γ:校正因子,孔隙截面为圆形时γ=1,孔隙截面为正方形时γ=1.094,孔隙截面为等边三角形时γ=1.186;
D:孔隙特征长度,单位为m;
(2)针对页岩储层泄压动态开发过程中,孔隙尺寸演化,提出如下计算方法:
由于气体解吸导致基质收缩,孔隙半径增大:
上式中:drdis:基质收缩引起的孔隙半径变化,单位为m;
rin:初始孔隙等效半径(Dequ/2),单位为m;
Φint:初始孔隙度;
εL:Langmuir应变;
PL:Langmuir压力,单位为Pa;
Pin:初始孔隙压力,单位为Pa;
Pp:实时孔隙压力,单位为Pa;
由于应力敏感导致孔隙半径变小:
上式中:drpm:应力敏感引起的孔隙半径变化,单位为m;
αf:页岩孔隙的Biot系数;
αm:页岩基质的Biot系数;
Kn:页岩孔隙平均方向刚度,单位为Pa/m;
Em:页岩基质杨氏模量,单位为Pa;
s:页岩基质特征长度,单位为m;
Pin:初始孔隙压力,单位为Pa;
Pp:实时孔隙压力,单位为Pa;
由于吸附层气体解吸导致吸附层厚度变小,有效孔隙半径变大:
ra=ppdm/(pL+pp)
上式中:ra:吸附层厚度,单位为m;
dm:气体分子直径,单位为m;
PL:Langmuir压力,单位为Pa;
Pp:实时孔隙压力,单位为Pa;
三重机理协同作用下实时孔隙等效半径为:
rp=rin+drdis+drpm-ra=f(γD,pp)
上式中:rp:实时孔隙等效半径,单位为m;
基于以上计算,实时孔隙等效半径为孔隙截面几何形态、孔隙尺寸及孔隙压力的函数。(3)针对其多尺度孔隙结构,提出如下粗化计算流程:
基于分形粗化理论,孔隙尺寸大于或等于D的累积孔隙数N:
上式中:N:累积孔隙数目;
Dmin:孔隙特征长度下限,单位为m;
Dmax:孔隙特征长度上限,单位为m;
Df为孔隙尺寸分形因子,通过下式进行计算:
基于分形理论,计算在[D,D+dD]之间的孔隙数-dN:
毛管长度,单位为m:
孔隙体积,单位为m3:
上式中:ωi:不同截面孔隙占比,可通过页岩薄片扫描电镜图片进行统计得到;
Vpi:不同截面孔隙对应的孔隙体积,单位为m3;
基质体积(考虑基质收缩):
上式中:ωi:不同截面孔隙占比;
Vmi:不同截面孔隙对应的基质体积,单位为m3;
动态孔隙度:
单根毛管表观渗透率:
μ:气体粘度,单位为pa.s;
ρavg:气体密度,单位为kg/m;
R:气体常数;
α:切向动量调节系数;
T:温度,单位为K;
ωi:不同截面孔隙占比;
Ds:气体表面扩散系数,单位为m2/s;
ζms:表面扩散修正因子;
Cs:吸附气体浓度,单位为kg/m3;
M:气体分子量,单位为kg/mol;
基于改进的Kozeny-Carman方程计算页岩表观渗透率为:
K=ΦpKa
根据基于改进的Kozeny-Carman方程计算页岩表观渗透率对页岩储层渗透率进行预测。
Claims (10)
1.基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,包括如下步骤:
S1,针对页岩储层孔隙截面几何形态多样性进行归一化校正,将页岩储层孔隙截面为非圆形孔隙校正为圆形孔隙,得到孔隙等效直径;
S2,针对页岩气藏泄压动态开发过程中的孔隙尺寸演化,利用孔隙等效直径将孔隙尺寸演化进行定量表征,得到实时孔隙等效半径;
S3,基于分形粗化理论确定孔隙尺寸分形因子;根据孔隙尺寸分形因子,计算毛管长度以及预设范围的孔隙数,利用实时孔隙等效半径和毛管长度计算孔隙体积,利用毛管长度以及预设范围的孔隙数计算基质体积,利用孔隙体积和基质体积计算动态孔隙度,利用实时孔隙等效半径计算不同截面的单根毛管表观渗透率,根据不同截面孔隙占比以及不同截面的单根毛管表观渗透率计算单根毛管平均表观渗透率,根据单根毛管平均表观渗透率和动态孔隙度计算得到改进后的Kozeny-Carman方程,从而计算得到页岩表观渗透率对页岩储层渗透率进行预测。
2.根据权利要求1所述的基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,S1中,通过校正因子γ将截面为非圆形的孔隙校正为圆形孔隙,校正后的孔隙等效直径Dequ为:
Dequ=γD
上式中:
γ为校正因子,孔隙截面为圆形时γ=1,孔隙截面为正方形时γ=1.094,孔隙截面为等边三角形时γ=1.186;
D为孔隙特征长度。
3.根据权利要求1所述的基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,S2中,实时孔隙等效半径rp为:
rp=rin+drdis+drpm-ra=f(γD,pp)
其中,drdis为基质收缩引起的孔隙半径变化;drpm为应力敏感引起的孔隙半径变化;ra为吸附层厚度。
6.根据权利要求3所述的基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,吸附层厚度ra为:
ra=ppdm/(pL+pp)
其中,dm为气体分子直径;PL为Langmuir压力;Pp为实时孔隙压力。
8.根据权利要求1所述的基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,S3中,孔隙体积Vpi为:
S3中,基质体积Vm为:
10.根据权利要求1所述的基于改进Kozeny-Carman模型的页岩储层渗透率预测方法,其特征在于,S3中,基于改进的Kozeny-Carman方程计算页岩表观渗透率K为:
K=ΦpKa
Φp为动态孔隙度,Ka为考虑不同截面单根毛管的平均表观渗透率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010130359.6A CN111208052B (zh) | 2020-02-28 | 2020-02-28 | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010130359.6A CN111208052B (zh) | 2020-02-28 | 2020-02-28 | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111208052A true CN111208052A (zh) | 2020-05-29 |
CN111208052B CN111208052B (zh) | 2021-08-13 |
Family
ID=70788619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010130359.6A Active CN111208052B (zh) | 2020-02-28 | 2020-02-28 | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111208052B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220519A (zh) * | 2018-11-23 | 2020-06-02 | 中国石油天然气股份有限公司 | 标准岩心模型及其制造方法 |
CN113029864A (zh) * | 2021-03-03 | 2021-06-25 | 重庆科技学院 | 一种循环测试聚合物溶液动态吸附量的监测装置及方法 |
CN113128040A (zh) * | 2021-04-12 | 2021-07-16 | 陕西延长石油(集团)有限责任公司 | 一种基于孔隙结构特征的页岩表观渗透率预测方法 |
CN113188976A (zh) * | 2021-05-11 | 2021-07-30 | 中国石油大学(华东) | 一种夹层状页岩各向异性渗透率确定方法及系统 |
CN114136862A (zh) * | 2021-11-29 | 2022-03-04 | 西南石油大学 | 一种双润湿性页岩的液体表观渗透率计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181219A (zh) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | 一种基于碳酸盐岩孔隙结构预测渗透率的方法及其装置 |
CN108229048A (zh) * | 2018-01-17 | 2018-06-29 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑页岩基质差异的多场耦合渗透率计算方法 |
CN108710723A (zh) * | 2018-04-09 | 2018-10-26 | 西南石油大学 | 一种页岩储层多孔介质表观渗透率的计算方法 |
CN108763711A (zh) * | 2018-05-22 | 2018-11-06 | 中国石油大学(华东) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 |
CN109375283A (zh) * | 2018-11-07 | 2019-02-22 | 中国石油大学(华东) | 一种砂岩储层3d渗透率演化史的分析方法 |
-
2020
- 2020-02-28 CN CN202010130359.6A patent/CN111208052B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181219A (zh) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | 一种基于碳酸盐岩孔隙结构预测渗透率的方法及其装置 |
CN108229048A (zh) * | 2018-01-17 | 2018-06-29 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑页岩基质差异的多场耦合渗透率计算方法 |
CN108710723A (zh) * | 2018-04-09 | 2018-10-26 | 西南石油大学 | 一种页岩储层多孔介质表观渗透率的计算方法 |
CN108763711A (zh) * | 2018-05-22 | 2018-11-06 | 中国石油大学(华东) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 |
CN109375283A (zh) * | 2018-11-07 | 2019-02-22 | 中国石油大学(华东) | 一种砂岩储层3d渗透率演化史的分析方法 |
Non-Patent Citations (4)
Title |
---|
李玉丹等: "基于分形理论的页岩基质表观渗透率研究", 《油气地质与采收率》 * |
段永刚等: "页岩储层纳米孔隙流动模拟研究", 《西南石油大学学报(自然科学版)》 * |
白瑞婷等: "考虑启动压力梯度的致密砂岩储层渗透率分形模型", 《天然气地球科学》 * |
耿黎东: "基于微纳尺度流动的页岩储层产能预测模型与完井参数优化", 《中国博士学位论文全文数据库工程科技I辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220519A (zh) * | 2018-11-23 | 2020-06-02 | 中国石油天然气股份有限公司 | 标准岩心模型及其制造方法 |
CN113029864A (zh) * | 2021-03-03 | 2021-06-25 | 重庆科技学院 | 一种循环测试聚合物溶液动态吸附量的监测装置及方法 |
CN113029864B (zh) * | 2021-03-03 | 2021-12-14 | 重庆科技学院 | 一种循环测试聚合物溶液动态吸附量的监测装置及方法 |
CN113128040A (zh) * | 2021-04-12 | 2021-07-16 | 陕西延长石油(集团)有限责任公司 | 一种基于孔隙结构特征的页岩表观渗透率预测方法 |
CN113128040B (zh) * | 2021-04-12 | 2024-02-20 | 陕西延长石油(集团)有限责任公司 | 一种基于孔隙结构特征的页岩表观渗透率预测方法 |
CN113188976A (zh) * | 2021-05-11 | 2021-07-30 | 中国石油大学(华东) | 一种夹层状页岩各向异性渗透率确定方法及系统 |
CN114136862A (zh) * | 2021-11-29 | 2022-03-04 | 西南石油大学 | 一种双润湿性页岩的液体表观渗透率计算方法 |
CN114136862B (zh) * | 2021-11-29 | 2024-04-09 | 西南石油大学 | 一种双润湿性页岩的液体表观渗透率计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111208052B (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111208052B (zh) | 基于改进Kozeny-Carman模型的页岩储层渗透率预测方法 | |
CN111208051B (zh) | 复杂机理协同影响下页岩储层渗透率动态预测方法 | |
CN108959789B (zh) | 一种考虑应力敏感效应的页岩表观渗透率计算方法 | |
CN109902918B (zh) | 一种页岩压裂自支撑裂缝缝网渗透率计算方法 | |
CN107461192B (zh) | 储层条件下页岩动态表观渗透率的计算方法 | |
CN108229048B (zh) | 一种考虑页岩基质差异的多场耦合渗透率计算方法 | |
CN110685651B (zh) | 一种多层合采气井产量劈分方法及系统 | |
CN111340298B (zh) | 页岩储层诱导非均质性动态演化预测方法 | |
CN110472372B (zh) | 基于双重介质的渗透率预测方法及系统 | |
CN109033541B (zh) | 一种基于eur的压后页岩气藏非均质性评价方法 | |
Anderberg et al. | Method for simultaneous determination of sorption isotherms and diffusivity of cement-based materials | |
CN108982320B (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
CN108629471A (zh) | 储气库气井交互注采的方法 | |
CN113128040B (zh) | 一种基于孔隙结构特征的页岩表观渗透率预测方法 | |
CN112102487A (zh) | 一种基于多重混合分形的非常规储层三维渗透率确定方法 | |
CN108280288B (zh) | 基于应力敏感性分析的储层渗透率模型建立方法 | |
CN108416116B (zh) | 拱脚增大截面加固抛物线拱时拱顶压重量的确定方法 | |
CN108071392B (zh) | 一种海上异常高压气藏动态储量计算方法 | |
CN115017691A (zh) | 基于页岩孔隙性质及结构特征的表观渗透率预测方法及系统 | |
CN115713163A (zh) | 深层页岩气压裂水平井产量及可采储量预测方法 | |
CN114169204A (zh) | 一种用于海上油气田开发生产的防砂时机确定方法 | |
CN112967147B (zh) | 一种考虑多尺度裂缝的基岩产量贡献率计算方法 | |
CN114462690A (zh) | 一种基于突变理论的深层页岩气压裂效果预测方法 | |
CN108487087B (zh) | 拱脚加固时拱顶压重水箱水量实时控制方法 | |
CN112112639B (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 |