CN113534290A - 一种部分饱和岩石声电性质联合模拟方法 - Google Patents

一种部分饱和岩石声电性质联合模拟方法 Download PDF

Info

Publication number
CN113534290A
CN113534290A CN202110812572.XA CN202110812572A CN113534290A CN 113534290 A CN113534290 A CN 113534290A CN 202110812572 A CN202110812572 A CN 202110812572A CN 113534290 A CN113534290 A CN 113534290A
Authority
CN
China
Prior art keywords
rock
conductivity
specific condition
partially saturated
saturated
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
Application number
CN202110812572.XA
Other languages
English (en)
Other versions
CN113534290B (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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202110812572.XA priority Critical patent/CN113534290B/zh
Publication of CN113534290A publication Critical patent/CN113534290A/zh
Application granted granted Critical
Publication of CN113534290B publication Critical patent/CN113534290B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • 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

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种部分饱和岩石声电性质联合模拟方法，其中该方法包括以下步骤：计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率；计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率；计算部分饱和岩石的纵波速度；计算岩石胶结指数；计算部分饱和岩石电导率；根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。本发明公开的部分饱和岩石声电性质联合模拟方法，在保证了声、电性质模拟过程中所构建的部分饱和岩石具有统一的微观结构的前提下，用以解决岩石骨架与孔隙不能同时连通的问题。

Description

一种部分饱和岩石声电性质联合模拟方法
技术领域
本发明涉及地球物理勘探技术领域，具体涉及一种部分饱和岩石声电性质联合模拟方法。
背景技术
储层流体饱和度的精确计算对于油气勘探非常重要。由于地层流体饱和度对岩石的弹性和电学性质的影响不可忽视，可以采用联合解释声波和电磁勘探数据的方式提高饱和度的预测精度，但是其关键前提是量化不同饱和度对声、电联合关系的影响。
目前，声学建模中所采取的岩石微观结构构建方法为在固体背景中加入孔隙包含物的方式；而电学建模中常利用导电流体中包含固体颗粒的方式。然而，建立岩石声电性质联合模拟方法的关键要素是声学模型与电学模型所面对的岩石具有一致的微观结构。由于现有的部分饱和岩石声、电模型不能保证岩石骨架与孔隙双方都相互连通，互相独立的部分饱和声电模型无法直接用来模拟部分饱和岩石的声电联合性质，致使目前建立的岩石物理模型还不能准确模拟饱和度对岩石声电联合性质的影响。
发明内容
本发明的目的在于提供一种部分饱和岩石声电性质联合模拟方法，用于解决目前建立的岩石物理模型还不能准确模拟饱和度对岩石声电联合性质的影响的问题。
本发明提供一种部分饱和岩石声电性质联合模拟方法，包括以下步骤：
步骤S1：计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
步骤S2：以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
步骤S3：向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
步骤S4：根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
步骤S5：向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数和孔隙水的电导率计算部分饱和岩石电导率；
步骤S6：根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
优选地，在所述步骤S1中，通过自洽模型计算孔隙度为40％所述特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率。
优选地，在所述步骤S2中，通过用微分等效介质模型计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，将特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率的结果推广至孔隙度为[0，1]全阈值任意孔隙度值下的干燥岩石体积模量和剪切模量以及饱和岩石的电导率。
优选地，在所述步骤S3中，将岩石孔隙中的部分饱和流体看作向特定条件或非特定条件下的干燥岩石加入由水包裹球状气泡的理想模型，使用斑块饱和模型计算部分饱和岩石纵波速度。
优选地，在所述步骤S4中，通过有效介质非对称电阻率模型计算得到所述岩石胶结指数。
优选地，所述步骤S5包括以下步骤：
首先，将饱和岩石中孔隙的内含物看作水包裹电导率为0的球状气泡的模型；
然后，加入气泡，使之变成部分饱和岩石，并统计气泡加入量；
其次，根据微分等效介质模型和气泡加入量计算该部分饱和岩石的气水混合物的电导率σmix；
最后，用气水混合物的电导率σmix代替孔隙水的电导率σw，结合岩石胶结指数m和矿物颗粒的电导率σg以及岩石孔隙度φ代入有效介质非对称电阻率模型计算部分饱和岩石电导率σ。
本发明还提供一种部分饱和岩石声电性质联合模拟装置，用于实现上述部分饱和岩石声电性质联合模拟方法，包括：
特定条件下处理单元，用于计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
非特定条件下处理单元，用于以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
部分饱和岩石的纵波速度处理单元，用于向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
岩石胶结指数处理单元，用于根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
部分饱和岩石电导率处理单元，用于向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数和孔隙水的电导率计算部分饱和岩石电导率；
部分饱和岩石声电性质联合关系处理单元，用于根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述部分饱和岩石声电性质联合模拟方法的步骤。
本发明还提供一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述部分饱和岩石声电性质联合模拟方法的步骤。
与现有技术相比，本发明的有益效果是：
本发明公开了一种部分饱和岩石声电性质联合模拟方法，该在保证了声、电性质模拟过程中所构建的部分饱和岩石具有统一的微观结构，还实现了岩石骨架与孔隙都各自均为相互连通的整体，把部分饱和岩石孔隙空间中所充填的流体视为水包含气泡的混合相流体，在保证了声、电性质模拟过程中所构建的部分饱和岩石具有统一的微观结构的前提下，用以解决岩石骨架与孔隙不能同时连通的问题。
附图说明
图1为本发明实施例1提供的部分饱和岩石声电联合性质的计算模型流程图；
图2为本发明实施例1提供的模型预测结果中纵波速度与饱和度的关系图；
图3为本发明实施例1提供的模型预测结果中电导率与饱和度的关系图。
具体实施方式
以下实施例用于说明本发明，但不用来限制本发明的范围。
岩石物理模型是研究饱和度对岩石声、电学性质影响规律的重要手段。White通过将部分饱和岩石孔隙中的流体假设为由水包含球状气泡的方式，建立了能定量表征部分饱和岩石声学性质的斑块饱和模型(又称White模型，英文全称Patchy saturation model，其原文出处为White,J.E.,1975,Computed seismic speeds and attenuation in rockswith partial gas saturation:Geophysics,40,224–232.)。Berg基于有效介质非对称电阻率模型(又称HB模型，英文全称Hanai-Berryman model，其原文出处为Hanai,T.,1961,Aremark on“Theory of the dielectric dispersion due to the interfacialpolarization and its application to emulsions”:Kolloid-Zeitschrift,175,61–62.)，建立了通过岩石电学性质表征岩石饱和度的增量模型，阐明了岩石电学性质与饱和度之间关系，这些理论模型的建立极大地推动了饱和度对岩石声电性质的影响机制和规律的发展。然而，以上成果都仅仅面向饱和度对岩石声学或电学性质等某一方面性质的影响，至于饱和度对岩石声电联合性质的影响，相关的岩石物理模拟方法仍然缺乏，导致针对地层流体饱和度评价的声电勘探资料联合解释缺少相应的理论支撑。
本发明公开了一种部分饱和岩石声电性质联合模拟方法，首先，联合自洽模型(又称为SC模型，英文全称Self-consistent model，其原文出处为Berryman,J.G.,1995,Mixture theories for rock properties,in T.J.Ahrens,ed.,Rock physics and phaserelations:American Geophysical Union,205–208.)和微分等效介质模型(又称DEM模型，英文全称Differential equivalent medium model，其原文出处为Asami,K.,2002,Characterization of heterogeneous systems by dielectric spectroscopy:Progressin Polymer Science,27,1617–1659.)计算得到干燥岩石的体积模量和剪切模量以及饱和岩石的电导率；然后，使用斑块饱和模型计算得到部分饱和岩石的纵波速度；再联合微分等效介质模型和有效介质非对称电阻率模型计算得到部分饱和岩石电导率；最后，以相同的含水饱和度为限定计算部分饱和岩石的纵波速度和电导率性质。
实施例1：一种部分饱和岩石声电性质联合模拟方法，
实施例1提供一种部分饱和岩石声电性质联合模拟方法，参考图1，包括以下步骤：
步骤S1：计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
具体地，干燥岩石的含水量为0，所含流体为天然气，其体积模量为0.101MPa，剪切模量为0；饱和岩石所含流体为35g/L的盐水，其电导率为4.69S/m。
通过自洽模型计算孔隙度为40％这一特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中，自洽模型具有以下公式：
(1-φ)(Km-KSC)Pm(αm)+φ(Kg-KSC)Ppore(αpore)＝0 (式1)
(1-φ)(μm-μSC)Qm(αm)+φ(μg-μSC)Qpore(αpore)＝0 (式2)
(1-φ)(σm-σSC)Rm(αm)+φ(σw-σSC)Rpore(αpore)＝0 (式3)
其中，Kg、μg和σw分别表示天然气的体积模量，剪切模量以及水的电导率；下标SC代表整体岩石的性质；下标m代表纵横比为αm的矿物颗粒的性质；下标g和w分别表示纵横比为αpore的孔隙中所充填的气或水的性质；Pm、Qm和Rm以及Ppore、Qpore和Rpore分别是与矿物颗粒纵横比和孔隙纵横比有关的参数组。
步骤S2：以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
具体地，通过用微分等效介质模型计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，将特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率的结果推广至孔隙度为[0，1]全阈值任意孔隙度值下的干燥岩石体积模量和剪切模量以及饱和岩石的电导率。
其中，微分等效介质模型具有以下公式：
Figure BDA0003168769520000051
Figure BDA0003168769520000052
Figure BDA0003168769520000053
式中，Ki是各组分(矿物颗粒或者孔隙)的体积模量，μi是各组分的剪切模量，fi是各组分的体积分数，Pi、Qi以及Ri是与矿物颗粒和孔隙纵横比有关的参数。
具体地，对于孔隙度小于40％的目标岩石，采用添加矿物颗粒的方式降低孔隙度；对于孔隙度大于40％的目标岩石，采用添加孔隙的方式增大孔隙度。
步骤S3：向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
将岩石孔隙中的部分饱和流体看作向特定条件或非特定条件下的干燥岩石加入由水包裹球状气泡的理想模型，使用斑块饱和模型计算部分饱和岩石纵波速度。
其中，干燥岩石的体积模量和剪切模量具有如下形式：
Figure BDA0003168769520000054
Figure BDA0003168769520000055
式中，Kpar和μpar分别表示所求部分饱和岩石的体积模量和剪切模量；K∞表示高频极限下不包含流体流动效应的特征单元的体积模量；a，b分别为气泡和水所占据区域的半径；ω为角频率；μ1和μ2分别代表干燥及饱和岩石的剪切模量，可以通过步骤S2获得；Sw表示含水饱和度；
K∞可由干燥岩石和饱和岩石的体积模量K1、K2以及饱和岩石的剪切模量μ2以及含水饱和度Sw表示：
Figure BDA0003168769520000061
其中K1、K2和μ1、μ2都可以通过步骤S2获得；
R1，R2，Q1，Q2，Z1，Z2是与岩石物理参数相关的函数，分别具有以下形式：
R1＝0 (式10)
Figure BDA0003168769520000062
Figure BDA0003168769520000063
Figure BDA0003168769520000064
Figure BDA0003168769520000065
Figure BDA0003168769520000066
其中κ为岩石渗透率；η1、η2分别为气和水的黏度；γ1、γ2分别为气和水所在区域的复传播常数，分别具有如下形式：
Figure BDA0003168769520000067
Figure BDA0003168769520000068
通过斑块饱和模型计算部分饱和岩石体积模量与剪切模量都为复数形式，分别取其实部带入公式18，可以进一步得到部分饱和岩石的纵波速度：
Figure BDA0003168769520000069
其中密度ρpar具有如下形式：
ρpar＝ρm(1-φ)+φ(1-Sw)ρg+φSwρw (式19)
步骤S4：根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
通过有效介质非对称电阻率模型计算得到岩石胶结指数，其中有效介质非对称电阻率模型的表达式：
Figure BDA0003168769520000071
其中，m代表岩石胶结指数，σ0代表饱和岩石的电导率，σw代表孔隙水的电导率，σg代表矿物颗粒的电导率，φ为岩石孔隙度。
步骤S5：向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数m和孔隙水的电导率σw代入有效介质非对称电阻率模型计算部分饱和岩石电导率；
首先，将饱和岩石中孔隙的内含物看作水包裹电导率为0的球状气泡的模型；
然后，加入气泡，使之变成部分饱和岩石，并统计气泡加入量；
其次，根据微分等效介质模型和气泡加入量计算该部分饱和岩石孔隙中的气水混合物的电导率σmix；
最后，用气水混合物的电导率σmix代替孔隙水的电导率σw，结合岩石胶结指数m和矿物颗粒的电导率σg以及岩石孔隙度φ代入有效介质非对称电阻率模型计算部分饱和岩石电导率σ；
步骤S6：根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
根据步骤S3中的纵波速度和步骤S5中的部分饱和岩石电导率σ，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
在保证了声、电性质模拟过程中所构建的部分饱和岩石具有统一的微观结构，还对岩石微观结构的刻画采取了“双连通结构”(bi-connected microstructure)，即岩石骨架与孔隙都各自均为相互连通的整体，把部分饱和岩石孔隙空间中所充填的流体视为水包含气泡的混合相流体，在保证了声、电性质模拟过程中所构建的部分饱和岩石具有统一的微观结构的前提下，用以解决岩石骨架与孔隙不能同时连通的问题。
该实施例分别对部分饱和岩石纵波速度和电导率的预测值进行了验证，取得了较好的验证效果，如图2和图3所示。其中模型的输入参数分别如下：矿物密度2.56g/cm3；孔隙度13.78％；渗透率20.37mD；斑块尺寸0.189mm；矿物颗粒纵横比0.88；孔隙纵横比0.34。图中模型预测曲线与实验测量值之间良好的拟合关系证明本发明提出的方法在部分饱和岩石声电性质联合模拟方面具有较高的精确度。
实施例2：一种部分饱和岩石声电性质联合模拟装置
实施例2提供一种部分饱和岩石声电性质联合模拟装置，用于实现实施例1的部分饱和岩石声电性质联合模拟方法，包括：
特定条件下处理单元，用于计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
非特定条件下处理单元，用于以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
部分饱和岩石的纵波速度处理单元，用于向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
岩石胶结指数处理单元，用于根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
部分饱和岩石电导率处理单元，用于向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数和孔隙水的电导率计算部分饱和岩石电导率；
部分饱和岩石声电性质联合关系处理单元，用于根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
实施例3：一种计算机可读存储介质
实施例3提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现实施例1的部分饱和岩石声电性质联合模拟方法的步骤。
实施例4：一种计算机设备
实施例4提供一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现实施例1的部分饱和岩石声电性质联合模拟方法的步骤。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (9)

1.一种部分饱和岩石声电性质联合模拟方法，其特征在于，包括以下步骤：
特定条件下处理：计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
非特定条件下处理：以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
部分饱和岩石的纵波速度处理：向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
岩石胶结指数处理：根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
部分饱和岩石电导率处理：向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数和孔隙水的电导率计算部分饱和岩石电导率；
部分饱和岩石声电性质联合关系处理：根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
2.如权利要求1所述的部分饱和岩石声电性质联合模拟方法，其特征在于，
在所述特定条件下处理中，通过自洽模型计算孔隙度为40％所述特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率。
3.如权利要求1所述的部分饱和岩石声电性质联合模拟方法，其特征在于，
在所述非特定条件下处理中，通过用微分等效介质模型计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，将特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率的结果推广至孔隙度为[0，1]全阈值任意孔隙度值下的干燥岩石体积模量和剪切模量以及饱和岩石的电导率。
4.如权利要求1所述的部分饱和岩石声电性质联合模拟方法，其特征在于，
在所述部分饱和岩石的纵波速度处理中，将岩石孔隙中的部分饱和流体看作向特定条件或非特定条件下的干燥岩石加入由水包裹球状气泡的理想模型，使用斑块饱和模型计算部分饱和岩石纵波速度。
5.如权利要求1所述的部分饱和岩石声电性质联合模拟方法，其特征在于，
在所述岩石胶结指数处理中，通过有效介质非对称电阻率模型计算得到所述岩石胶结指数。
6.如权利要求1所述的部分饱和岩石声电性质联合模拟方法，其特征在于，所述部分饱和岩石电导率处理包括以下步骤：
首先，将饱和岩石中孔隙的内含物看作水包裹电导率为0的球状气泡的模型；
然后，加入气泡，使之变成部分饱和岩石，并统计气泡加入量；
其次，根据微分等效介质模型和气泡加入量计算该部分饱和岩石的气水混合物的电导率σmix；
最后，用气水混合物的电导率σmix代替孔隙水的电导率σw，结合岩石胶结指数m和矿物颗粒的电导率σg以及岩石孔隙度φ代入有效介质非对称电阻率模型计算部分饱和岩石电导率σ。
7.一种部分饱和岩石声电性质联合模拟装置，其特征在于，包括：
特定条件下处理单元，用于计算特定条件下的干燥岩石的体积模量和剪切模量以及饱和岩石的电导率，其中该特定条件为孔隙度为40％；
非特定条件下处理单元，用于以特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率微分等效介质模型的基础，计算非特定条件下的干燥岩石的体积模量和剪切模量和饱和岩石的电导率，其中该非特定条件为孔隙度为40％以外；
部分饱和岩石的纵波速度处理单元，用于向特定条件或非特定条件下的干燥岩石加入由水包裹的球状气泡，使之变成部分饱和岩石，根据特定条件或非特定条件下的干燥岩石的体积模量和剪切模量，计算部分饱和岩石的纵波速度；
岩石胶结指数处理单元，用于根据特定条件或非特定条件下的饱和岩石电导率计算岩石胶结指数；
部分饱和岩石电导率处理单元，用于向特定条件或非特定条件下饱和岩石加入气泡，使之变成部分饱和岩石，计算水气混合物的电导率，用气水混合物的电导率代替孔隙水的电导率，根据岩石胶结指数和孔隙水的电导率计算部分饱和岩石电导率；
部分饱和岩石声电性质联合关系处理单元，用于根据所述部分饱和岩石的纵波速度和所述部分饱和岩石电导率，通过设置相同的含水饱和度，得到部分饱和岩石声电性质联合关系。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-6中任一项所述的部分饱和岩石声电性质联合模拟方法的步骤。
9.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-6任一所述的部分饱和岩石声电性质联合模拟方法的步骤。
CN202110812572.XA 2021-07-19 2021-07-19 一种部分饱和岩石声电性质联合模拟方法 Active CN113534290B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110812572.XA CN113534290B (zh) 2021-07-19 2021-07-19 一种部分饱和岩石声电性质联合模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110812572.XA CN113534290B (zh) 2021-07-19 2021-07-19 一种部分饱和岩石声电性质联合模拟方法

Publications (2)

Publication Number Publication Date
CN113534290A true CN113534290A (zh) 2021-10-22
CN113534290B CN113534290B (zh) 2023-05-16

Family

ID=78128663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110812572.XA Active CN113534290B (zh) 2021-07-19 2021-07-19 一种部分饱和岩石声电性质联合模拟方法

Country Status (1)

Country Link
CN (1) CN113534290B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
WO2012035036A1 (en) * 2010-09-14 2012-03-22 Statoil Asa Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
WO2016051282A2 (en) * 2014-09-30 2016-04-07 King Abdullah University Of Science And Technology Reservoir resistivity characterization incorporating flow dynamics
CN106202763A (zh) * 2016-07-18 2016-12-07 西南石油大学 双重介质储层岩石含水饱和度的计算方法
CN109613624A (zh) * 2019-01-09 2019-04-12 中国石油大学(华东) 一种储层岩石声电性质联合模拟方法
CN111090125A (zh) * 2019-12-30 2020-05-01 中国石油大学(北京) 确定致密砂岩储层的弹性模量的方法和装置
CN112946783A (zh) * 2021-01-29 2021-06-11 中国石油大学(北京) 一种水合物饱和度确定方法、装置及设备

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
WO2012035036A1 (en) * 2010-09-14 2012-03-22 Statoil Asa Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
WO2016051282A2 (en) * 2014-09-30 2016-04-07 King Abdullah University Of Science And Technology Reservoir resistivity characterization incorporating flow dynamics
US20180231681A1 (en) * 2014-09-30 2018-08-16 King Abdullah University Of Science And Technology Reservoir resistivity characterization incorporating flow dynamics
CN106202763A (zh) * 2016-07-18 2016-12-07 西南石油大学 双重介质储层岩石含水饱和度的计算方法
CN109613624A (zh) * 2019-01-09 2019-04-12 中国石油大学(华东) 一种储层岩石声电性质联合模拟方法
CN111090125A (zh) * 2019-12-30 2020-05-01 中国石油大学(北京) 确定致密砂岩储层的弹性模量的方法和装置
CN112946783A (zh) * 2021-01-29 2021-06-11 中国石油大学(北京) 一种水合物饱和度确定方法、装置及设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
TONGCHENG HAN 等: "Are self-consistent models capable of jointly modeling elastic velocity and electrical conductivity of reservoir sandstones?", 《GEOPHYSICS》 *
任舒波 等: "不同压力下部分饱和砂岩纵波衰减的理论及实验研究", 《地球物理学报》 *
池美瑶 等: "由岩石电性和弹性参数求取地层储层参数", 《科学技术与工程》 *

Also Published As

Publication number Publication date
CN113534290B (zh) 2023-05-16

Similar Documents

Publication Publication Date Title
Borgomano et al. Dispersion and attenuation measurements of the elastic moduli of a dual‐porosity limestone
Carcione et al. Rock-physics templates for clay-rich source rocks
CN108957542B (zh) 一种建立地震波衰减岩石物理图板的方法
EP1820137A2 (en) Integrated anisotropic rock physics model
Ogushwitz Applicability of the Biot theory. II. Suspensions
CN109613624B (zh) 一种储层岩石声电性质联合模拟方法
Ruiz et al. Predicting elasticity in nonclastic rocks with a differential effective medium model
Wang et al. Wave properties of gas-hydrate bearing sediments based on poroelasticity
CN113534290A (zh) 一种部分饱和岩石声电性质联合模拟方法
Yu et al. Rock physics modeling of heterogeneous carbonate reservoirs: porosity estimation and hydrocarbon detection
Zheng et al. Sound speed, attenuation, and reflection in gassy sediments
Sayar et al. Effective medium modeling of velocity dispersion and attenuation in isotropic rocks
CN105301642B (zh) 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置
CN111273341B (zh) 根据裂缝空间分布的含裂缝储层岩石物理建模方法
Markov et al. Novel approach for simulating the elastic properties of porous rocks including the critical porosity phenomena
Sharma Wave propagation in a dissipative poroelastic medium
Boadu Predicting oil saturation from velocities using petrophysical models and artificial neural networks
CN112649861A (zh) 裂缝储层岩石物理建模方法及系统
Zhan et al. Rock-physics modeling of ultrasonic P-and S-wave attenuation in partially frozen brine and unconsolidated sand systems and comparison with laboratory measurements
Wang et al. A new rock physics model for tight reservoirs
CN110967746A (zh) 流体饱和度地震反演方法及系统
Wang et al. A full-frequency band Kuster-Toksöz model and its application in velocity dispersion analysis
CN112649854A (zh) 基于双尺度模型地震波频散与衰减预测方法及系统
Thomsen Fluid dependence of rock compressibility, post Biot-Gassmann
CN117316329B (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