CN112595653B - 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法 - Google Patents

粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法 Download PDF

Info

Publication number
CN112595653B
CN112595653B CN202110245188.6A CN202110245188A CN112595653B CN 112595653 B CN112595653 B CN 112595653B CN 202110245188 A CN202110245188 A CN 202110245188A CN 112595653 B CN112595653 B CN 112595653B
Authority
CN
China
Prior art keywords
clay
porosity
granular
granular clay
pore
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.)
Expired - Fee Related
Application number
CN202110245188.6A
Other languages
English (en)
Other versions
CN112595653A (zh
Inventor
李溢龙
李晓平
周克明
谭晓华
孟展
邹梦文
卢虹林
刘婷芝
刘蜀东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Exploration And Development Research Institute Of Southwest Oil And Gas Field Petrochina Co ltd
Southwest Petroleum University
Original Assignee
Exploration And Development Research Institute Of Southwest Oil And Gas Field Petrochina Co ltd
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Exploration And Development Research Institute Of Southwest Oil And Gas Field Petrochina Co ltd, Southwest Petroleum University filed Critical Exploration And Development Research Institute Of Southwest Oil And Gas Field Petrochina Co ltd
Priority to CN202110245188.6A priority Critical patent/CN112595653B/zh
Publication of CN112595653A publication Critical patent/CN112595653A/zh
Application granted granted Critical
Publication of CN112595653B publication Critical patent/CN112595653B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry

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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，属于油气田储层保护领域，它解决了现今孔隙度应力敏感分析方法少，现有方法须通过实验测定、可操作性低等问题；其技术方案是：根据岩心实验，获取粒状黏土矿物在多孔介质中的尺寸参数；根据岩心实验所获取得尺寸参数，建立粒状黏土物理模型；推导粒状黏土孔隙度有效应力系数数学模型；推导粒状黏土孔隙压缩系数数学模型；结合有效应力关系方程，推导粒状黏土孔隙度应力敏感模型；将岩心实验所获取得尺寸参数代入粒状黏土孔隙度应力敏感模型，进行孔隙度应力敏感分析；本发明具有以下有益效果：实验结合理论，可靠性强，准确度高；计算便捷，可推广性强。

Description

粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法
技术领域
本发明涉及一种粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，属于油气田储层应力敏感评价领域。
背景技术
在我国已开发的气藏中，大部分为中、低渗砂岩气藏，且其中低渗和特低渗气藏占据了非常大的比例。近几年来，随着常规气藏资源的逐渐减少，低渗和特低渗气藏的开发已经受到越来越多的关注。由于低渗气藏渗流过程受到“应力敏感”、“黏土矿物”因素的影响，这使得其渗流机理非常复杂，低渗气藏孔隙度应力敏感表现为随着有效应力的增加，储层孔隙度降低；低渗气藏在开发工程中具有产量下降快、单井产量低、稳产能力差等特点，与常规气藏的开发特点相比有很大的差异，这种应力敏感的存在进一步加大了低渗气藏开采的复杂性。
明确低渗储层孔隙应力敏感情况，对认识储层物性及预测油气储量具有重要意义，而目前有关孔隙度应力敏感分析的方法及手段极少，仅专利号CN202020079197.3的一种高温高压降内压孔隙度应力敏感性的测试装置，从实验入手，利用实验数据进行孔隙度应力敏感的分析，但该方法成本高、好时长，操作复杂且实验精度有限。因此，急需一种高效准确、可操作性强的孔隙度应力敏感分析方法进行富含黏土矿物疏松砂岩储层的应力敏感评价，为后续油气藏储量评价，高效开发提供支撑。
发明内容
本发明目的是：为了解决现今孔隙度应力敏感分析方法少，现有方法须通过实验测定、可操作性低等问题，本发明以基础实验同理论相结合的手段，提出了一种粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，计算准确，适用性强。
为实现上述目的，本发明提供了粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，其特征在于，该方法包括下列步骤：
S100、根据已获取岩样进行岩心实验，获取粒状黏土矿物在多孔介质中的尺寸参数及岩石参数，具体分为以下步骤；
S101、对岩样进行岩心实验，所述岩心实验包括铸体薄片鉴定、电镜扫描分析；
S102、获取粒状黏土矿物在多孔介质中的尺寸参数，所述尺寸参数包括孔隙内半 径为
Figure 100002_DEST_PATH_IMAGE002
、外半径为b、黏土颗粒半径为d；
S103、获取储层的岩石参数，所述岩石参数包括黏土矿物的泊松比v c、岩石骨架的泊松比v r、黏土矿物拉梅系数μc、岩石骨架拉梅系数μr。
S200、根据铸体薄片鉴定、电镜扫描分析得到的粒状黏土矿物在多孔介质中的尺寸参数，建立粒状黏土物理模型；
S201、所述粒状黏土物理模型，由岩石、孔隙、黏土三部分构成，其黏土矿物以颗粒的形式分散充填在砂岩的粒间孔隙中。
S300、根据所建立粒状黏土物理模型，推导粒状黏土孔隙度有效应力系数数学模型，具体分为以下步骤；
S301、根据孔隙度有效应力系数定义式，结合粒状黏土物理模型孔隙度定义式，通 过链式法进行展开，建立孔隙内壁面和孔隙外壁面的位移对孔隙压力和围压的偏导数关系 式及孔隙度对半径
Figure 100002_DEST_PATH_IMAGE004
、b和d的偏导数，进一步定义岩石骨架与黏土的刚度比及粒状黏土矿 物含量关系式，得到粒状黏土孔隙度有效应力系数方程；
S302、所述粒状黏土物理模型孔隙度定义式为
Figure 100002_DEST_PATH_IMAGE006
，式中
Figure 818722DEST_PATH_IMAGE002
为孔隙内半 径，单位为μm；b为孔隙外半径，单位为μm；d为黏土颗粒半径，单位为μm；φ为孔隙度，无量纲 量；
S303、所述粒状黏土孔隙度有效应力系数展开式为
Figure DEST_PATH_IMAGE008
；式中
Figure DEST_PATH_IMAGE010
为粒状黏土 孔隙度有效应力系数，无量纲量；φ为孔隙度，无量纲量；Pp为孔隙压力，单位为MPa；Pc为围 压，单位为MPa；
Figure DEST_PATH_IMAGE012
为黏土颗粒半径d对围压Pc的偏导数，
Figure DEST_PATH_IMAGE014
；
Figure DEST_PATH_IMAGE016
为黏 土颗粒半径d对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE018
；
Figure DEST_PATH_IMAGE020
为孔隙内半径
Figure 100002_DEST_PATH_IMAGE021
对孔隙压 力Pp的偏导数，
Figure DEST_PATH_IMAGE023
；
Figure DEST_PATH_IMAGE025
为孔隙内半径
Figure 640922DEST_PATH_IMAGE021
对围压Pc的偏导 数，
Figure DEST_PATH_IMAGE027
；
Figure DEST_PATH_IMAGE029
为孔隙外半径b对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE031
；
Figure DEST_PATH_IMAGE033
为孔隙外半径b对围压Pc的偏导数，
Figure DEST_PATH_IMAGE035
；v c为黏土矿物的泊松比，无量纲量；v r为岩石骨架的泊松 比，无量纲量；μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲量；
S304、对孔隙度定义式中的孔隙内半径
Figure 482976DEST_PATH_IMAGE004
、孔隙外半径b和黏土颗粒半径d求偏导 数，所述孔隙度对孔隙内半径
Figure 143765DEST_PATH_IMAGE021
、孔隙外半径b和黏土颗粒半径d的偏导数分别为
Figure DEST_PATH_IMAGE037
、
Figure DEST_PATH_IMAGE039
、
Figure DEST_PATH_IMAGE041
；
S305、定义岩石骨架与黏土的刚度比关系式，所述岩石骨架与黏土的刚度比关系 式为
Figure DEST_PATH_IMAGE043
，式中μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲 量；
S306、定义粒状黏土矿物含量关系式，所述粒状黏土矿物含量关系式为
Figure DEST_PATH_IMAGE045
；
S307、将粒状黏土矿物含量关系式、岩石骨架与黏土的刚度比关系式代入粒状黏 土孔隙度有效应力系数展开式中，得到粒状黏土孔隙度有效应力系数方程、所述粒状黏土 孔隙度有效应力系数方程为
Figure DEST_PATH_IMAGE047
，式 中M、N为无量纲量，
Figure DEST_PATH_IMAGE049
、
Figure DEST_PATH_IMAGE051
。
S400、基于所建立粒状黏土物理模型，推导粒状黏土孔隙压缩系数数学模型；
S401、所述粒状黏土孔隙压缩系数数学模型为
Figure DEST_PATH_IMAGE053
，式中C pp为粒状黏土孔隙压缩系数，单位为MPa-1。
S500、将粒状黏土孔隙度有效应力系数方程、粒状黏土孔隙压缩系数数学模型代入孔隙度与有效应力关系方程，推导粒状黏土孔隙度应力敏感模型；
S501、所述孔隙度与有效应力关系方程为
Figure DEST_PATH_IMAGE055
，Pp,0为地面条件下的孔隙压力，单位为MPa；φ 0为 地面条件下的孔隙度，无量纲量；
S502、所述粒状黏土孔隙度应力敏感模型为
Figure DEST_PATH_IMAGE057
，φ 0为地面条 件下测得初始孔隙度，无量纲量。
S600、根据所得粒状黏土孔隙度应力敏感模型，将粒状黏土矿物在多孔介质中的尺寸参数及岩石参数代入粒状黏土孔隙度应力敏感模型，进行疏松砂岩储层孔隙度应力敏感分析。
上述的粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法：所述方法除对富含黏土矿物的疏松砂岩适用外，当模型中黏土矿物含量F c为0时，对不含黏土矿物的疏松砂岩孔隙度应力敏感分析同样适用。
与现有技术相比，本发明具有以下有益效果：（1）实验结合理论，可靠性强；（2）基于粒状黏土矿物的赋存状态提出符合实际的理论模型，准确度高；（3）所需参数较少，计算便捷，可推广性强。
附图说明
图1是本方法技术路线图。
图2是粒状黏土物理模型图。
图3是无因次孔隙度随有效应力变化曲线图。
图中：
Figure 100002_DEST_PATH_IMAGE058
为孔隙内半径，单位为μm；b为孔隙外半径，单位为μm；d为黏土颗粒半径单 位为μm；P p为孔隙压力，单位为MPa；P c为围压，单位为MPa。
具体实施方式
下面结合实施方式和附图对本发明做进一步说明。
本发明提供了粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，图1为本方法的技术路线图，该方法包括下列步骤：
S100、根据已获取岩样进行岩心实验，获取粒状黏土矿物在多孔介质中的尺寸参数及岩石参数，具体分为以下步骤；
S101、对岩样进行岩心实验，所述岩心实验包括铸体薄片鉴定、电镜扫描分析；
S102、获取粒状黏土矿物在多孔介质中的尺寸参数，所述尺寸参数包括孔隙内半 径为
Figure 879508DEST_PATH_IMAGE058
、外半径为b、黏土颗粒半径为d；
S103、获取储层的岩石参数，所述岩石参数包括黏土矿物的泊松比v c、岩石骨架的泊松比v r、黏土矿物拉梅系数μc、岩石骨架拉梅系数μr。
S200、根据铸体薄片鉴定、电镜扫描分析得到的粒状黏土矿物在多孔介质中的尺寸参数，建立粒状黏土物理模型；
S201、所述粒状黏土物理模型，由岩石、孔隙、黏土三部分构成，其黏土矿物以颗粒的形式分散充填在砂岩的粒间孔隙中，如图2所示。
S300、根据所建立粒状黏土物理模型，推导粒状黏土孔隙度有效应力系数数学模型，具体分为以下步骤；
S301、根据孔隙度有效应力系数定义式，结合粒状黏土物理模型孔隙度定义式，通 过链式法进行展开，建立孔隙内壁面和孔隙外壁面的位移对孔隙压力和围压的偏导数关系 式及孔隙度对半径
Figure 967550DEST_PATH_IMAGE004
、b和d的偏导数，进一步定义岩石骨架与黏土的刚度比及粒状黏土矿 物含量关系式，得到粒状黏土孔隙度有效应力系数方程；
S302、所述粒状黏土物理模型孔隙度定义式为
Figure 234583DEST_PATH_IMAGE006
，式中
Figure 699063DEST_PATH_IMAGE058
为孔隙内半 径，单位为μm；b为孔隙外半径，单位为μm；d为黏土颗粒半径，单位为μm；φ为孔隙度，无量纲 量；
S303、所述粒状黏土孔隙度有效应力系数展开式为
Figure 633521DEST_PATH_IMAGE008
；式中
Figure 626885DEST_PATH_IMAGE010
为粒状黏土 孔隙度有效应力系数，无量纲量；φ为孔隙度，无量纲量；Pp为孔隙压力，单位为MPa；Pc为围 压，单位为MPa；
Figure DEST_PATH_IMAGE059
为黏土颗粒半径d对围压Pc的偏导数，
Figure DEST_PATH_IMAGE060
；
Figure 177952DEST_PATH_IMAGE016
为 黏土颗粒半径d对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE061
；
Figure 180543DEST_PATH_IMAGE020
为孔隙内半径
Figure 969507DEST_PATH_IMAGE021
对 孔隙压力Pp的偏导数，
Figure 664931DEST_PATH_IMAGE023
；
Figure 172135DEST_PATH_IMAGE025
为孔隙内半径
Figure 712838DEST_PATH_IMAGE021
对围压Pc 的偏导数，
Figure 621888DEST_PATH_IMAGE027
；
Figure 222634DEST_PATH_IMAGE029
为孔隙外半径b对孔隙压力Pp的偏导 数，
Figure 951556DEST_PATH_IMAGE031
；
Figure 295949DEST_PATH_IMAGE033
为孔隙外半径b对围压Pc的偏导数，
Figure DEST_PATH_IMAGE062
；v c为黏土矿物的泊松比，无量纲量；v r为岩石骨架的 泊松比，无量纲量；μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲量；
S304、对孔隙度定义式中的孔隙内半径
Figure 325085DEST_PATH_IMAGE004
、孔隙外半径b和黏土颗粒半径d求偏导 数，所述孔隙度对孔隙内半径
Figure 96732DEST_PATH_IMAGE021
、孔隙外半径b和黏土颗粒半径d的偏导数分别为
Figure 312950DEST_PATH_IMAGE037
、
Figure 461034DEST_PATH_IMAGE039
、
Figure 79097DEST_PATH_IMAGE041
；
S305、定义岩石骨架与黏土的刚度比关系式，所述岩石骨架与黏土的刚度比关系 式为
Figure 21646DEST_PATH_IMAGE043
，式中μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲量；
S306、定义粒状黏土矿物含量关系式，所述粒状黏土矿物含量关系式为
Figure DEST_PATH_IMAGE063
；
S307、将粒状黏土矿物含量关系式、岩石骨架与黏土的刚度比关系式代入粒状黏 土孔隙度有效应力系数展开式中，得到粒状黏土孔隙度有效应力系数方程、所述粒状黏土 孔隙度有效应力系数方程为
Figure DEST_PATH_IMAGE064
，式中M、N为无量纲量，
Figure 256318DEST_PATH_IMAGE049
、
Figure DEST_PATH_IMAGE065
。
S400、基于所建立粒状黏土物理模型，推导粒状黏土孔隙压缩系数数学模型；
S401、所述粒状黏土孔隙压缩系数数学模型为
Figure DEST_PATH_IMAGE066
，式中C pp为粒状黏土孔隙压缩系数，单 位为MPa-1。
S500、将粒状黏土孔隙度有效应力系数方程、粒状黏土孔隙压缩系数数学模型代入孔隙度与有效应力关系方程，推导粒状黏土孔隙度应力敏感模型；
S501、所述孔隙度与有效应力关系方程为
Figure 208093DEST_PATH_IMAGE055
，Pp,0为地面条件下孔隙压力，单位为MPa；φ 0为 地面条件下孔隙度，无量纲量；
S502、所述粒状黏土孔隙度应力敏感模型为
Figure 946242DEST_PATH_IMAGE057
，φ 0为地面条件 下测得初始孔隙度，无量纲量。
S600、根据所得粒状黏土孔隙度应力敏感模型，将粒状黏土矿物在多孔介质中的 尺寸参数及岩石参数代入粒状黏土孔隙度应力敏感模型，进行疏松砂岩储层孔隙度应力敏 感分析，得到岩心无因次孔隙度
Figure DEST_PATH_IMAGE068
随有效应力（Pc与Pp差值）的变化曲线图，如图3所示。
进一步的，所述方法除对富含黏土矿物的疏松砂岩适用外，当模型中黏土矿物含量F c为0时，对不含黏土矿物的疏松砂岩孔隙度应力敏感分析同样适用。
下面结合具体疏松砂岩岩心X01对粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法做进一步说明。
通过岩心实验包括铸体薄片鉴定、电镜扫描分析；得到X01岩心内粒状黏土矿物在多孔介质中的尺寸参数，见表1：
表1
Figure DEST_PATH_IMAGE070
利用上述所推导的粒状黏土孔隙度应力敏感模型：
Figure DEST_PATH_IMAGE071
，进行孔隙度应 力敏感分析，利用模型计算，绘制岩心无因次孔隙度
Figure DEST_PATH_IMAGE072
随有效应力（Pc与Pp差值）的变化曲 线图，得到如图3所示曲线，得到了岩心X01不同有效应力（Pc与Pp差值）下的岩心无因次孔隙 度大小，以有效应力是30MPa为例，此时无因次孔隙度
Figure DEST_PATH_IMAGE073
为0.871，表明30MPa有效应力下的 孔隙度为地面条件下孔隙度的87.1%，对比岩心X01通过实验手段测得有效应力是30MPa时 的孔隙度为地面条件下孔隙度的85.6%，误差为1.75%，该方法同实验结果又较好一致性，且 发明提供的粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法适用便捷，可操作性 强。
与现有技术相比，本发明具有以下有益效果：（1）实验结合理论，可靠性强；（2）基于黏土矿物的赋存状态提出符合实际的理论模型，准确度高；（3）所需参数较少，计算便捷，可推广性强。
最后所应说明的是：以上实施例仅用以说明而非限制本发明的技术方案，尽管参照上述实施例对本发明进行了详细说明，本领域的普通技术人员应该理解：依然可以对本发明进行修改或者等同替换，而不脱离本发明的精神和范围的任何修改或局部替换，其均应涵盖在本发明的权利要求范围当中。

Claims (2)

1.粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，其特征在于，该方法包括下列步骤：
S100、根据已获取岩样进行岩心实验，获取粒状黏土矿物在多孔介质中的尺寸参数及岩石参数，具体分为以下步骤；
S101、对岩样进行岩心实验，所述岩心实验包括铸体薄片鉴定、电镜扫描分析；
S102、获取粒状黏土矿物在多孔介质中的尺寸参数，所述尺寸参数包括孔隙内半径为
Figure DEST_PATH_IMAGE002
、外半径为b、黏土颗粒半径为d；
S103、获取储层的岩石参数，所述岩石参数包括黏土矿物的泊松比v c、岩石骨架的泊松比v r、黏土矿物拉梅系数μc、岩石骨架拉梅系数μr；
S200、根据铸体薄片鉴定、电镜扫描分析得到的粒状黏土矿物在多孔介质中的尺寸参数，建立粒状黏土物理模型；
S201、所述粒状黏土物理模型，由岩石、孔隙、黏土三部分构成，其黏土矿物以颗粒的形式分散充填在砂岩的粒间孔隙中；
S300、根据所建立粒状黏土物理模型，推导粒状黏土孔隙度有效应力系数数学模型，具体分为以下步骤；
S301、根据孔隙度有效应力系数定义式，结合粒状黏土物理模型孔隙度定义式，通过链式法进行展开，建立孔隙内壁面和孔隙外壁面的位移对孔隙压力和围压的偏导数关系式及孔隙度对半径
Figure DEST_PATH_IMAGE004
、b和d的偏导数，进一步定义岩石骨架与黏土的刚度比及粒状黏土矿物含量关系式，得到粒状黏土孔隙度有效应力系数方程；
S302、所述粒状黏土物理模型孔隙度定义式为
Figure DEST_PATH_IMAGE006
，式中
Figure DEST_PATH_IMAGE007
为孔隙内半径，单位为μm；b为孔隙外半径，单位为μm；d为黏土颗粒半径，单位为μm；φ为孔隙度，无量纲量；
S303、所述粒状黏土孔隙度有效应力系数展开式为
Figure DEST_PATH_IMAGE009
；式中
Figure DEST_PATH_IMAGE011
为粒状黏土孔隙度有效应力系数，无量纲量；φ为孔隙度，无量纲量；Pp为孔隙压力，单位为MPa；Pc为围压，单位为MPa；
Figure DEST_PATH_IMAGE013
为黏土颗粒半径d对围压Pc的偏导数，
Figure DEST_PATH_IMAGE015
；
Figure DEST_PATH_IMAGE017
为黏土颗粒半径d对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE019
；
Figure DEST_PATH_IMAGE021
为孔隙内半径
Figure DEST_PATH_IMAGE022
对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE024
；
Figure DEST_PATH_IMAGE026
为孔隙内半径
Figure 880207DEST_PATH_IMAGE022
对围压Pc的偏导数，
Figure DEST_PATH_IMAGE028
；
Figure DEST_PATH_IMAGE030
为孔隙外半径b对孔隙压力Pp的偏导数，
Figure DEST_PATH_IMAGE032
；
Figure DEST_PATH_IMAGE034
为孔隙外半径b对围压Pc的偏导数，
Figure DEST_PATH_IMAGE036
；v c为黏土矿物的泊松比，无量纲量；v r为岩石骨架的泊松比，无量纲量；μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲量；
S304、对孔隙度定义式中的孔隙内半径
Figure 538590DEST_PATH_IMAGE004
、孔隙外半径b和黏土颗粒半径d求偏导数，所述孔隙度对孔隙内半径
Figure 219232DEST_PATH_IMAGE022
、孔隙外半径b和黏土颗粒半径d的偏导数分别为
Figure DEST_PATH_IMAGE038
、
Figure DEST_PATH_IMAGE040
、
Figure DEST_PATH_IMAGE042
；
S305、定义岩石骨架与黏土的刚度比关系式，所述岩石骨架与黏土的刚度比关系式为
Figure DEST_PATH_IMAGE044
，式中μc为黏土矿物拉梅系数，无量纲量；μr为岩石骨架拉梅系数，无量纲量；
S306、定义粒状黏土矿物含量关系式，所述粒状黏土矿物含量关系式为
Figure DEST_PATH_IMAGE046
；
S307、将粒状黏土矿物含量关系式、岩石骨架与黏土的刚度比关系式代入粒状黏土孔隙度有效应力系数展开式中，得到粒状黏土孔隙度有效应力系数方程、所述粒状黏土孔隙度有效应力系数方程为
Figure DEST_PATH_IMAGE048
，式中M、N为无量纲量，
Figure DEST_PATH_IMAGE050
、
Figure DEST_PATH_IMAGE052
；
S400、基于所建立粒状黏土物理模型，推导粒状黏土孔隙压缩系数数学模型；
S401、所述粒状黏土孔隙压缩系数数学模型为
Figure DEST_PATH_IMAGE054
，式中C pp为粒状黏土孔隙压缩系数，单位为MPa-1；
S500、将粒状黏土孔隙度有效应力系数方程、粒状黏土孔隙压缩系数数学模型代入孔隙度与有效应力关系方程，推导粒状黏土孔隙度应力敏感模型；
S501、所述孔隙度与有效应力关系方程为
Figure DEST_PATH_IMAGE056
，Pp,0为地面条件下的孔隙压力，单位为MPa；φ 0为地面条件下的孔隙度，无量纲量；
S502、所述粒状黏土孔隙度应力敏感模型，
Figure DEST_PATH_IMAGE058
，φ 0为地面条件下测得初始孔隙度，无量纲量；
S600、根据所得粒状黏土孔隙度应力敏感模型，将粒状黏土矿物在多孔介质中的尺寸参数及岩石参数代入粒状黏土孔隙度应力敏感模型，进行疏松砂岩储层孔隙度应力敏感分析。
2.根据权利要求1所述的粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法，其特征在于：所述方法除对富含黏土矿物的疏松砂岩适用外，当粒状黏土物理模型中黏土矿物含量F c为0时，对不含黏土矿物的疏松砂岩孔隙度应力敏感分析同样适用。
CN202110245188.6A 2021-03-05 2021-03-05 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法 Expired - Fee Related CN112595653B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110245188.6A CN112595653B (zh) 2021-03-05 2021-03-05 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110245188.6A CN112595653B (zh) 2021-03-05 2021-03-05 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法

Publications (2)

Publication Number Publication Date
CN112595653A CN112595653A (zh) 2021-04-02
CN112595653B true CN112595653B (zh) 2021-05-07

Family

ID=75210156

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110245188.6A Expired - Fee Related CN112595653B (zh) 2021-03-05 2021-03-05 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法

Country Status (1)

Country Link
CN (1) CN112595653B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117150183B (zh) * 2023-10-27 2024-02-13 西北工业大学青岛研究院 海底沉积物孔隙度计算方法、装置、电子设备及介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5684299A (en) * 1995-06-26 1997-11-04 Schlumberger Technology Corporation Method for determining porosity in an invaded gas reservoir
CN103776748A (zh) * 2014-02-14 2014-05-07 武汉科技大学 宾汉姆流体在多孔介质中的有效渗透率的预测方法
CN105719063A (zh) * 2016-01-15 2016-06-29 西南石油大学 一种页岩气储量品质分类的综合评价方法
CN107145696A (zh) * 2017-06-29 2017-09-08 中国石油大学(北京) 一种煤层气地上地下耦合求解的模拟方法
CA3041838A1 (en) * 2016-10-27 2018-05-03 Schlumberger Canada Limited Evaluation of formation utilizing wideband electromagnetic measurements
CN108804779A (zh) * 2018-05-24 2018-11-13 西安石油大学 一种基于数字岩心模型的稠油油藏储层解堵过程模拟方法
CN108875140A (zh) * 2018-05-24 2018-11-23 西安石油大学 一种基于数字岩心模型的稠油油藏沥青质沉积吸附损害模拟方法
CN109085112A (zh) * 2018-10-08 2018-12-25 中国石油天然气股份有限公司 致密岩样的渗透率测定方法及装置
CN109975189A (zh) * 2017-12-28 2019-07-05 中国石油天然气股份有限公司 孔隙型砂岩储层产能预测方法与装置
CN110617062A (zh) * 2019-10-14 2019-12-27 西南石油大学 一种储层剩余气分布特征的分析方法
CN112255159A (zh) * 2020-10-10 2021-01-22 西南石油大学 一种非均质页岩储层水化作用能力的综合评价方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6651741B2 (en) * 2001-10-13 2003-11-25 1407580 Ontario Inc. Method of increasing productivity of oil, gas and hydrogeological wells
CN102953718B (zh) * 2011-08-23 2015-03-25 中国石油化工股份有限公司 低渗透油层复合交替酸溶方法
CN106223918B (zh) * 2016-08-18 2019-08-20 西南石油大学 水力压裂破裂压力获得方法及装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5684299A (en) * 1995-06-26 1997-11-04 Schlumberger Technology Corporation Method for determining porosity in an invaded gas reservoir
CN103776748A (zh) * 2014-02-14 2014-05-07 武汉科技大学 宾汉姆流体在多孔介质中的有效渗透率的预测方法
CN105719063A (zh) * 2016-01-15 2016-06-29 西南石油大学 一种页岩气储量品质分类的综合评价方法
CA3041838A1 (en) * 2016-10-27 2018-05-03 Schlumberger Canada Limited Evaluation of formation utilizing wideband electromagnetic measurements
CN107145696A (zh) * 2017-06-29 2017-09-08 中国石油大学(北京) 一种煤层气地上地下耦合求解的模拟方法
CN109975189A (zh) * 2017-12-28 2019-07-05 中国石油天然气股份有限公司 孔隙型砂岩储层产能预测方法与装置
CN108804779A (zh) * 2018-05-24 2018-11-13 西安石油大学 一种基于数字岩心模型的稠油油藏储层解堵过程模拟方法
CN108875140A (zh) * 2018-05-24 2018-11-23 西安石油大学 一种基于数字岩心模型的稠油油藏沥青质沉积吸附损害模拟方法
CN109085112A (zh) * 2018-10-08 2018-12-25 中国石油天然气股份有限公司 致密岩样的渗透率测定方法及装置
CN110617062A (zh) * 2019-10-14 2019-12-27 西南石油大学 一种储层剩余气分布特征的分析方法
CN112255159A (zh) * 2020-10-10 2021-01-22 西南石油大学 一种非均质页岩储层水化作用能力的综合评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Claystone porosity and mechanical behavior vs.geologic burial stress;Russell Ewy;《Marine and Petroleum Geology》;20200711;第121卷;1-20 *
层状储层中孔隙度和粘土体积的确定;N Ruhowets;《国外地质勘探技术》;19910831(第8期);37-39 *

Also Published As

Publication number Publication date
CN112595653A (zh) 2021-04-02

Similar Documents

Publication Publication Date Title
Song et al. Multi-fractal distribution analysis for pore structure characterization of tight sandstone—A case study of the Upper Paleozoic tight formations in the Longdong District, Ordos Basin
Xiao et al. Impacts of clay on pore structure, storage and percolation of tight sandstones from the Songliao Basin, China: implications for genetic classification of tight sandstone reservoirs
CN106869896B (zh) 一种低渗透油气藏直井体积压裂储层改造体积预测方法
CN104990851A (zh) 一种新的页岩敏感性实验研究方法
CN112595653B (zh) 粒状黏土矿物的疏松砂岩储层孔隙度应力敏感分析方法
US20090228249A1 (en) Method of modelling a saturation dependant property in a sample
Li et al. Mechanical properties of the clayey-silty sediment-natural gas hydrate mixed system
CN111027818B (zh) 一种页岩油分类评价方法
CN104833789A (zh) 利用致密砂岩微观孔隙结构确定气水关系的方法
CN103674802A (zh) 岩石封闭孔隙度测定方法
Cai et al. Pore-throat structures of the Permian Longtan Formation tight sandstones in the South Yellow Sea Basin, China: A case study from borehole CSDP-2
Hu et al. The influence of inhomogeneous hydrate formation on permeability anisotropy of porous media
Sneider et al. Rock types, depositional history, and diagenetic effects, Ivishak reservoir, Prudhoe Bay field
Li et al. Pore‐Throat Combination Types and Gas‐Water Relative Permeability Responses of Tight Gas Sandstone Reservoirs in the Zizhou Area of East Ordos Basin, China
Elenbaas et al. A radial turbulent flow formula
Wang et al. Experimental study on fractal characteristics of fault filling medium in the tunnel and relationship between fractal dimension and permeability coefficient
CN112836944A (zh) 一种深井、超深井井壁稳定评价模型的建立方法
Song et al. Research on reservoir characteristics of Chang7 tight oil based on nano-CT
Zhou et al. An experimental investigation of the effect of acid stimulation on gas extraction from coal
CN115017779A (zh) 一种测量压裂缝几何参数的方法及***
Eggleston et al. The value of grain‐size hydraulic conductivity estimates: Comparison with high resolution in‐situ field hydraulic conductivity
Chen et al. A simple water retention model of dual-porosity soils based on self-similarity of bimodal pore size distribution
Zhao et al. Pore-scale hydraulic properties of virtual sandstone microstructures: spatial variations and voxel scale effects
Villarroel et al. Breakouts: physical and numerical modeling
Guo-Quan et al. Bimodal SWCC and Bimodal PSD of Soils with Dual‐Porosity Structure

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210507