CN102508296A - 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 - Google Patents
一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 Download PDFInfo
- Publication number
- CN102508296A CN102508296A CN2011103601853A CN201110360185A CN102508296A CN 102508296 A CN102508296 A CN 102508296A CN 2011103601853 A CN2011103601853 A CN 2011103601853A CN 201110360185 A CN201110360185 A CN 201110360185A CN 102508296 A CN102508296 A CN 102508296A
- Authority
- CN
- China
- Prior art keywords
- phi
- centerdot
- rho
- modulus
- wave
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000011435 rock Substances 0.000 claims abstract description 74
- 239000012530 fluid Substances 0.000 claims abstract description 65
- 230000004044 response Effects 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 34
- 238000004458 analytical method Methods 0.000 claims description 23
- 239000012071 phase Substances 0.000 claims description 23
- 238000010008 shearing Methods 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 16
- 230000035699 permeability Effects 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 7
- 239000007791 liquid phase Substances 0.000 claims description 6
- 239000007790 solid phase Substances 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000009432 framing Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 230000008450 motivation Effects 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 2
- 108010014172 Factor V Proteins 0.000 claims 1
- 239000000463 material Substances 0.000 claims 1
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 5
- 239000011707 mineral Substances 0.000 abstract description 5
- 238000005553 drilling Methods 0.000 abstract description 2
- 230000035515 penetration Effects 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 21
- 238000010586 diagram Methods 0.000 description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 12
- 239000000243 solution Substances 0.000 description 9
- 238000011160 research Methods 0.000 description 4
- 239000004575 stone Substances 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000012797 qualification Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001421775 Thereus Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- JLQUFIHWVLZVTJ-UHFFFAOYSA-N carbosulfan Chemical compound CCCCN(CCCC)SN(C)C(=O)OC1=CC=CC2=C1OC(C)(C)C2 JLQUFIHWVLZVTJ-UHFFFAOYSA-N 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007790 scraping Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种非饱和双重孔隙介质地震波频散衰减分析方法及装置,该方法包括:步骤1:获取包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,并根据所述地质数据生成储层岩石干骨架模型;步骤2:获取包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据所述测量数据及Van derWaals方程生成储层流体模型;步骤3:根据所述储层岩石干骨架模型、储层流体模型及描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程求解平面波解,生成纵波、横波的相速度及逆品质因子;步骤4:根据所述纵波、横波的相速度及逆品质因子生成速度、衰减、波阻抗、AVO响应特征等与频率之间的关系。
Description
技术领域
本发明是关于地球物理勘探中的地震岩石物理模型及地震波数值模拟技术,特别是关于一种非饱和双重孔隙介质地震波频散衰减分析方法及装置。
背景技术
上地壳沉积岩在形成过程中普遍具有一定的孔隙度,随着距地表深度的增加,在压实、固结与充填作用下,岩石的孔隙度会显著降低,但是,局部高孔隙度发育的砂层仍然会为气态烃的储集创造条件。在含水层以下,岩石孔隙中同时饱和气、水两相流体的情况非常普遍,因此,有必要深入探索非饱和岩石中的弹性、电阻率、磁性等物理性质,为实际的地球探测与矿产勘查工作提供依据。近年来,中国天然气供应的需求急剧增加,为满足人民生产生活的需要,生产部门为天然气田的勘测与开发提出了更高的技术要求,为达到提高气层检测精度、逐步实现定量识别的目的,开展深入研究以探索含气、水两相流体非饱和岩石中的地震波传播与响应规律的工作势在必行。
早期针对含流体多孔介质中的弹性波传播规律的研究,在理论上主要基于宏观均匀性假设,近似认为流体与固体均匀分布于岩石内部的每一个力学微元中,每一个力学单元内部仅含有一个流体相与一个固体相,而介质内部的所有力学微元都是相同的、均一的。这样的做法明显忽略了岩石内部的局部非均匀性,在实际的应用中也被发现明显低估了地震频段内(几十到几百Hz)岩石内部的波速频散与能量衰减现象。
基于宏观均一的Biot理论,White等人在1975年首次引入了岩石内部小尺度非均匀性的概念,考虑了含水岩石内部局部分布的气泡对地震波传播的影响。在实际岩石中,局部气体集中形成的非浸润状态的气泡在地震波的挤压作用下会发生耗散性的振动,即局部流体流动效应。在低频端,局部流体流动的充分发生会造成地震波能量的大量损失,并使得岩石变得松弛,而在高频地震波的激励下,由于流体在单个振荡周期内不能完成完整的局部荡动,岩石会呈现出比较“硬”的等效状态。Dutta和Sheriff在1979年对White理论进行了改进,改善了White理论的低频预测结果与Gassmann理论零频极限的吻合状况。Shapiro和Muller在1999年调查了局部流体流动中的流体压力与岩石参数的关系,得出中观尺度非均匀性诱发的局部流体流动是导致低频地震波能量强衰减的主因结论,2006年,Carcione等的研究显示,若渗透率下降、流体黏性上升或地层厚度增加,则地震波的衰减与频散曲线有向低频段移动的趋势。
聂建新等人在2004年基于同时包含Biot摩擦机制与喷射流机制的非饱和多孔隙BISQ模型,利用小生境遗传算法实现了储层参数(孔隙度、渗透率、含气饱和度等)的反演。刘炯等人在2009年用孔隙介质力学的方法,研究了嵌入有球状气泡的含水岩石模型中的地震波传播规律。聂建新、杨顶辉与巴晶在2010年借鉴等效介质的思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含非饱和流体孔隙介质中波传播问题的黏弹性BISQ模型。目前,研究者们在相关问题上已基本达到一致的认识,即局部流体流动是导致孔隙岩石中地震波速度频散与能力衰减的决定性因素,然而,至今为止,文献中还没有出现一个能够普遍被各种研究与产业部门接受的描述局域流的基础理论与定量技术。
综上所述,现有技术的研究中主要存在以下问题:
(1)基于黏弹性理论与BISQ理论的相关研究中,由于引入了一些并不具有明确物理意义或并不易于直接物理实现的参数或系数(如各类松弛时间、黏弹性系数或特征喷射流长度等),使得这些理论虽然在进行岩石波速现象的描述与预测中非常有效,但相关理论与核心参数的数学基础与物理内涵却难以得到实现或验证,因此即使抛开实验验证,这类方法在实际工程中很难得到推广应用。
(2)在基于精细模型的数值模拟与地震响应分析的相关研究中,由于需要建立过于详细的岩石模型并设置边界条件,使得这类方法虽然能够较准确的给出实际岩石中的地震波响应,但在实际工程中却往往由于缺乏足够的先验信息,并且由于相关方法本身建模的高复杂性与算法的高消耗性,使得此类方法也无法在实际工程应用中得到很好的实现。
(3)某些地震岩石模型在动力学方程组的推导过程中采用了统计热力学关系,使得方程的一些核心参数在计算过程中仅能通过频率域的关系式进行上、下限的估算,而无法做到更准确的地震响应的定量预测。
发明内容
本发明提供一种非饱和双重孔隙介质地震波频散衰减分析方法及系统,以分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度、含气饱和度的变化规律。
为了实现上述目的,本发明提供一种非饱和双重孔隙介质地震波频散衰减分析方法,该方法包括:步骤1:获取包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,并根据所述地质数据生成储层岩石干骨架模型;步骤2:获取包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据所述测量数据及Van der Waals方程生成储层流体模型;步骤3:根据所述储层岩石干骨架模型、储层流体模型及描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程求解平面波解,生成纵波、横波的相速度及逆品质因子;步骤4:根据所述纵波、横波的相速度及逆品质因子生成速度、衰减、波阻抗、AVO响应特征等与频率之间的关系。
Km,μm——骨架基质的体积模量和剪切模量;
Ki,μi——第i个包体的体积模量和剪切模量;
xi——每种包体的体积比率;
N——包体的个数;
ζm=μm(9Km+8μm)/6(Km+2μm);
Pmi,Qmi——骨架基质m中加入包体材料i后的弹性模量。
进一步地,在所述步骤2中,生成储层流体模型的步骤包括采用所述Vander Waals方程计算超临界状态下甲烷气体的密度和体积模量:
P——地层压力;
T——地层温度;
a——气体有关的参数,对于甲烷a=0.225Pa(m3/mol)2=879.9Pa(m3/kg)2;
b——气体有关的参数,对于甲烷b=42.7cm3/mol=2.675×10-3m3/kg;
ρg——气体密度;
Kg——气体体积模量;
R——气体常数,R=8.31J/(mol°K)=519.4J/(kg°K)。
进一步地,在所述步骤3中,描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程如下:
ui、Ui——分别为固相和液相的位移;
e、ξ——分别为固相、液相的体应变;
ζ——地震波激励过程中产生的一种局域性的流体变形增量;
ρf1——背景相流体的密度;
R0——嵌入体(气泡)半径;
bi——(i=1,2)分别表示两种流体饱和区域中各自的耗散系数;
κ10——表示主骨架(背景相)中的渗透率;
η1——表示背景相中的流体黏度;
φ——岩石中总的孔隙度,φ=φ1+φ2=v1φ10+v2φ20,v1与v2分别表示两种流体所占据的骨架占总骨架的体积比率,φ10与φ20分别表示两个区域内部的局部孔隙度;
A、N、Q1、R1、Q2、R2——双孔介质中的六个Biot弹性参数;
N=μb
ks,kb,μs,μb——分别为岩石基质和骨架的体积模量、剪切模量;
ρ11、ρ12、ρ13、ρ22、ρ33——表示双孔介质中的五个密度参数,满足下式:
(1-φ1-φ2)ρs=ρ11+ρ12+ρ13
进一步地,在所述步骤3中,所述生成纵波、横波的相速度及逆品质因子,包括:以位移场ui和Ui表示在非饱和双重孔隙介质中传播的弹性平面波的通解形式,将位移场平面波解 j=x,y,z.代入到所述Biot-Rayleigh方程,得到Christoffel方程,在A和有非零解的情况下,根据Christoffel方程的行列式等于零得到纵波和横波的频散关系,最后由复波数可求得弹性波的相速度Vi和逆品质因子
为实现上述目的,本发明还提供一种非饱和双重孔隙介质地震波频散衰减分析装置,该装置包括:岩石骨架弹性模量生成模块,用于根据包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,采用Voigt-Reuss-Hill平均模型生成骨架基质的等效弹性模量,并应用等效介质理论或者接触理论计算干燥岩石的等效体积模量和剪切模量;储层流体模型建立模块,用于根据包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据流体相态分析,应用Van der Waals方程或者流体试验测量结果建立储层流体模型;非饱和双重孔隙介质地震岩石物理建模模块,用于输入岩石骨架模型和流体模型,应用双重孔隙介质Biot-Rayleigh理论建立非饱和储层岩石物理模型;频散特性分析模块,用于分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度及含气饱和度的变化规律。
本发明实施例的有益效果在于,本发明基于Rayleigh理论描述纵波激励下非饱和岩石中气泡的局域流体流动,从经典力学的哈密顿原理导出双重孔隙介质中的波传播方程,即Biot-Rayleigh方程,该方程格式简洁,所有相关参数物理可测,具有较好的物理可实现性,可用于分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度、含气饱和度的变化规律。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例非饱和双重孔隙介质地震波频散衰减分析方法流程图;
图2为本发明实施例岩石干骨架建模流程图;
图3为本发明实施例储层流体建模流程图;
图4为本发明实施例非饱和双重孔隙介质地震波频散衰减计算流程图;
图5为本发明实施例非饱和双重孔隙介质地震波频散衰减分析装置结构框图;
图6及图7分别为三种理论预测纵波速度、衰减的结果对比示意图;
图8及图9分别为频率为0.1~10MHz,饱和度为0~100%时纵波速度和衰减变化示意图;
图10为实验室观测数据和Biot-Rayleigh理论预测数据在地震频率,纵波速度随饱和度变化趋势图;
图11为Biot-Rayleigh理论预测的含70%甲烷水砂岩纵波速度示意图;
图12为预测的横波速度示意图;
图13为泥岩+70%甲烷水砂岩的AVO特征曲线示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
如图1所示,本实施例提供一种非饱和双重孔隙介质地震波频散衰减分析方法,该方法包括:
步骤S101:获取包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,并根据所述地质数据生成储层岩石干骨架模型。
图2为本发明岩石干骨架建模(建立储层岩石干骨架模型)流程图,其具体实现原理如下:
已知岩石中各种矿物成分与其弹性模量,而矿物间的组合方式未知的情况下,Voigt界限和Reuss界限给出了岩石弹性模量的上下限。
Voigt界限又被称为等应变平均,其计算式如下:
Reuss界限又被称为等应力平均,其计算式如下:
(1)、(2)两个公式中,N表示构成成分的个数,fi为第i个介质的体积分量,Mi为第i个介质的弹性模量,MV表示Voigt上限,MR表示Reuss下限。M可以代表任何一种模量:K,μ,E等。
Voigt-Reuss-Hill平均模型是Voigt上限和Reuss下限二者的算术平均,其表达式如下:
(3)式中,MVRH即为Voigt-Reuss-Hill平均模量,可以表示任何一种弹性模量,也是其后计算岩石干骨架模量时所需的岩石骨架基质模量。
Km,μm——骨架基质的体积模量和剪切模量;
Ki,μi——第i个包体的体积模量和剪切模量;
xi——每种包体的体积比率;
N——包体的个数;
ζm=μm(9Km+8μm)/6(Km+2μm);
Pmi,Qmi——骨架基质m中加入包体材料i后的弹性模量。
步骤S102:获取包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据所述测量数据及Van der Waals方程生成储层流体模型。
图3为本发明实施例储层流体建模流程图,在步骤S102中,采用Van derWaals方程计算目的层温度、压力下的甲烷气体的密度和体积模量:
其中,公式(6)、(7)P——地层压力;
T——地层温度;
a——气体有关的参数,对于甲烷a=0.225Pa(m3/mol)2=879.9Pa(m3/kg)2;
b——气体有关的参数,对于甲烷b=42.7cm3/mol=2.675×10-3m3/kg;
ρg——气体密度;
Kg——气体体积模量;
R——气体常数,R=8.31J/(mol°K)=519.4J/(kg°K)。
步骤S103:根据所述储层岩石干骨架模型、储层流体模型及描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程求解平面波解,生成纵波、横波的相速度及逆品质因子。
图4为本发明实施例非饱和双重孔隙介质地震波频散衰减计算流程图,描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程组如下:
(8)式中,ui、Ui——分别为固相和液相的位移;
e、ξ——分别为固相、液相的体应变;
ζ——地震波激励过程中产生的一种局域性的流体变形增量;
ρf1——背景相流体的密度;
R0——嵌入体(气泡)半径;
bi——(i=1,2)分别表示两种流体饱和区域中各自的耗散系数;
κ10——表示主骨架(背景相)中的渗透率;
η1——表示背景相中的流体黏度;
φ——岩石中总的孔隙度,φ=φ1+φ2=v1φ10+v2φ20,v1与v2分别表示两种流体所占据的骨架占总骨架的体积比率,φ10与φ20分别表示两个区域内部的局部孔隙度;
A、N、Q1、R1、Q2、R2——双孔介质中的六个Biot弹性参数;
N=μb
ks,kb,μs,μb——分别为岩石基质和骨架的体积模量、剪切模量;
ρ11、ρ12、ρ13、ρ22、ρ23——表示双孔介质中的五个密度参数,满足下式:
(1-φ1-φ2)ρs=ρ11+ρ12+ρ13 (14)
接下来求解Biot-Rayleigh方程组的平面波解,得到纵波和横波的相速度Vi和逆品质因子
步骤S104:根据所述纵波、横波的相速度及逆品质因子生成速度、衰减、波阻抗、AVO响应特征等与频率之间的关系,根据上式可以分析得到非饱和岩石中的纵波频散与衰减随频率、孔隙度、含气饱和度的变化规律。
如图5所示,本发明提供一种非饱和双重孔隙介质地震波频散衰减分析装置,该非饱和双重孔隙介质地震波频散衰减分析装置包括:岩石骨架弹性模量生成模块501,储层流体模型建立模块502,非饱和双重孔隙介质地震岩石物理建模模块503及频散特性分析模块504。
岩石骨架弹性模量生成模块501用于根据包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,采用Voigt-Reuss-Hill平均模型生成骨架基质的等效弹性模量,并应用等效介质理论或者接触理论计算干燥岩石的等效体积模量和剪切模量。
根据目的区的文献调研、地质报告、测井资料、岩芯照片与电镜分析得到矿物组分、矿物体积比率、渗透率、孔隙率、泥质含量等,采用Voigt-Reuss-Hill平均模型计算骨架基质的等效弹性模量,应用等效介质理论或者接触理论计算干燥岩石的等效体积模量和剪切模量,如上述公式(1)至(5)所示。
储层流体模型建立模块502用于根据包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据流体相态分析,应用Van der Waals方程或者流体试验测量结果建立储层流体模型。
根据钻井、录井、测井资料和文献调研等得到目的区的温度、压力信息,根据流体相态分析,应用Van der Waals方程或者流体试验测量结果建立储层流体模型,如上述公式(6)、(7)所示.
非饱和双重孔隙介质地震岩石物理建模模块503用于输入岩石骨架模型和流体模型,应用双重孔隙介质Biot-Rayleigh理论建立非饱和储层岩石物理模型,如公式(8)至(21)所示。
频散特性分析模块504用于分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度及含气饱和度的变化规律。
下面以具体的实例详细说明本发明:
实例1
Biot-Rayleigh理论、Johnson理论和White理论预测纵波速度、衰减的结果对比:
为检查拓展到非饱和岩石状况的Biot-Rayleigh理论的合理性,设计一个非饱和岩石物理模型,应用Biot-Rayleigh理论、Johnson理论和White三种理论进行了纵波速度和能量衰减的预测,对预测结果进行了对比分析。岩石参数为:基质体积模量35GPa,骨架体积模量7Gpa,骨架剪切模量9Gpa,水体积模量2.25GPa,空气体积模量0.1MPa,水黏度0.001Pa*s,空气黏度0.00001Pa*s,基质平均密度2.65g/cm3,水密度0.99g/cm3,空气密度0.1g/cm3,孔隙度为0.15,渗透率0.1D,平均气泡尺寸0.25m,空气饱和度为0.95。图6及图7分别为三种理论预测纵波速度、衰减的结果对比示意图。
实例2
全频带非饱和岩石纵波的速度与衰减的定量预测:
设计一个非饱和岩石物理模型,基本的岩石物理参数为:石英体积模量39GPa,砂岩固结系数15,水体积模量2.51×109Pa,气体积模量1.44×105Pa,水黏度0.001Pa*s,气黏度0.000022Pa*s,砂岩颗粒平均密度2.69g/cm3,水密度1.04g/cm3,气密度0.01g/cm3,平均气泡尺寸lmm,孔隙度为0.11,渗透率0.5mD。图8及图9分别为频率为0.1~10MHz,饱和度为0~100%时纵波速度和衰减变化示意图。
实例3
北海高孔砂岩宽频带非饱和岩石波速观测数据对比:
为验证本文理论成果的适用性,采用了Adam与Batzle于2006年发表的宽频带非饱和岩石波速观测数据与本发明Biot-Rayleigh理论预测的结果进行对比分析。高孔砂岩来自北海,主要成分为石英,孔隙度为0.35,渗透率为8.7D,气泡平均尺寸为5mm,图10为实验室观测数据和Biot-Rayleigh理论预测数据在地震频率(<100Hz)、测井频率(50KHz)、超声频率下(500kHz),纵波速度随饱和度变化趋势图。
实例4
含非饱和甲烷砂岩纵波速度、衰减及AVO特征:
以南海渐新统陵水组长石岩屑质石英砂岩为例设计一个非饱和甲烷岩石物理模型,基本的岩石物理参数为:石英体积模量37GPa,剪切模量44Gpa,密度2.65g/cm3;长石体积模量37.5GPa,剪切模量15Gpa,密度2.62g/cm3;粘土体积模量21GPa,剪切模量7Gpa,密度2.6g/cm3;水体积模量2.5×109Pa,水密度1.04g/cm3,埋深3.5Km下甲烷体积模量0.081GPa,密度0.17g/cm3;孔隙度为0.15,甲烷饱和度为70%,盖层为泥岩,纵波速度3400m/s,横波速度1600m/s。图11为Biot-Rayleigh理论预测的含70%甲烷水砂岩纵波速度示意图,图12为预测的横波速度示意图,图13为泥岩+70%甲烷水砂岩的AVO特征曲线示意图。
本发明实施例的有益效果在于,本发明基于Rayleigh理论描述纵波激励下非饱和岩石中气泡的局域流体流动,从经典力学的哈密顿原理导出双重孔隙介质中的波传播方程,即Biot-Rayleigh方程,该方程格式简洁,所有相关参数物理可测,具有较好的物理可实现性,可用于分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度、含气饱和度的变化规律。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种非饱和双重孔隙介质地震波频散衰减分析方法,其特征在于,所述的方法包括:
步骤1:获取包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,并根据所述地质数据生成储层岩石干骨架模型;
步骤2:获取包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据所述测量数据及Van der Waals方程生成储层流体模型;
步骤3:根据所述储层岩石干骨架模型、储层流体模型及分析非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程求解平面波解,生成纵波、横波的相速度及逆品质因子;
步骤4:根据所述纵波、横波的相速度及逆品质因子生成速度、衰减、波阻抗、AVO响应特征等与频率之间的关系。
3.根据权利要求1所述的方法,其特征在于:
在所述步骤2中,生成储层流体模型的步骤包括采用所述Van der Waals方程计算超临界状态下甲烷气体的密度和体积模量:
P——地层压力;
T——地层温度;
a——气体有关的参数,对于甲烷a=0.225Pa(m3/mol)2=879.9Pa(m3/kg)2;
b——气体有关的参数,对于甲烷b=42.7cm3/mol=2.675×10-3m3/kg;
ρg——气体密度;
Kg——气体体积模量;
R——气体常数,R=8.31J/(mol°K)=519.4J/(kg°K)。
4.根据权利要求1所述的方法,其特征在于:
在所述步骤3中,描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程如下:
ui、Ui——分别为固相和液相的位移;
e、ξ——分别为固相、液相的体应变;
ζ——地震波激励过程中产生的一种局域性的流体变形增量;
ρf1——背景相流体的密度;
R0——嵌入体(气泡)半径;
bi——(i=1,2)分别表示两种流体饱和区域中各自的耗散系数;
κ10——表示主骨架(背景相)中的渗透率;
η1——表示背景相中的流体黏度;
φ——岩石中总的孔隙度,φ=φ1+φ1=v1φ10+v2φ20,v1与v2分别表示两种流体所占据的骨架占总骨架的体积比率,φ10与φ10分别表示两个区域内部的局部孔隙度;
A、N、Q1、R1、Q2、R2——双孔介质中的六个Biot弹性参数;
N=μb
ks,kb,μs,μb——分别为岩石基质和骨架的体积模量、剪切模量;
ρ11、ρ12、ρ13、ρ22、ρ23——表示双孔介质中的五个密度参数,满足下式:
(1-φ1-φ2)ρs=ρ11+ρ12+ρ13
6.一种非饱和双重孔隙介质地震波频散衰减分析装置,其特征在于,所述装置包括:
岩石骨架弹性模量生成模块,用于根据包括渗透率、孔隙率、泥质含量及矿物成分的地质数据,采用Voigt-Reuss-Hill平均模型生成骨架基质的等效弹性模量,并应用等效介质理论或者接触理论计算干燥岩石的等效体积模量和剪切模量;
储层流体模型建立模块,用于根据包括钻井数据、录井数据、测井信息及流体实验测量数据的测量数据,根据流体相态分析,应用Van der Waals方程或者流体试验测量结果建立储层流体模型;
非饱和双重孔隙介质地震岩石物理建模模块,用于输入岩石骨架模型和流体模型,应用双重孔隙介质Biot-Rayleigh理论建立非饱和储层岩石物理模型;
频散特性分析模块,用于分析非饱和岩石中的纵横波频散与衰减随频率、孔隙度及含气饱和度的变化规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103601853A CN102508296B (zh) | 2011-11-14 | 2011-11-14 | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103601853A CN102508296B (zh) | 2011-11-14 | 2011-11-14 | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508296A true CN102508296A (zh) | 2012-06-20 |
CN102508296B CN102508296B (zh) | 2013-07-31 |
Family
ID=46220401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103601853A Active CN102508296B (zh) | 2011-11-14 | 2011-11-14 | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508296B (zh) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103163553A (zh) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | 基于多重孔隙介质模型的地震烃类检测方法及装置 |
CN103616716A (zh) * | 2013-12-04 | 2014-03-05 | 中国石油天然气集团公司 | 一种双孔隙物理模型及其制作方法与应用 |
CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
CN103792574A (zh) * | 2012-10-30 | 2014-05-14 | 中国石油化工股份有限公司 | 一种储层频变含气性检测方法 |
CN104237936A (zh) * | 2013-06-06 | 2014-12-24 | 中国石油化工股份有限公司 | 一种油气检测的频变反演方法 |
CN104360383A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种地震波衰减预测的方法及系统 |
CN104573150A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种针对富含有机质页岩的各向异性岩石物理建模方法 |
CN105607123A (zh) * | 2016-01-15 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种计算随机孔隙介质模型地震波特征信息的方法和装置 |
CN105759310A (zh) * | 2014-12-16 | 2016-07-13 | 中国石油化工股份有限公司 | 复杂非均质储层介质中地震波衰减及速度频散预测方法 |
CN106324672A (zh) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及系统 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107942375A (zh) * | 2017-11-17 | 2018-04-20 | 河海大学 | 一种基于声波方程的非线性优化隐式时空域有限差分数值模拟方法 |
CN108957542A (zh) * | 2018-07-10 | 2018-12-07 | 河海大学 | 一种建立地震波衰减岩石物理图板的方法 |
CN109116420A (zh) * | 2018-10-16 | 2019-01-01 | 河海大学 | 一种含裂隙的孔隙介质纵波速度与衰减预测方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
CN109521090A (zh) * | 2018-10-18 | 2019-03-26 | 天津大学 | 一种激光无损表征薄膜杨氏模量的优化方法 |
CN110219643A (zh) * | 2019-06-05 | 2019-09-10 | 中国石油天然气集团有限公司 | 油藏含水饱和度确定方法及装置 |
CN110456412A (zh) * | 2019-07-09 | 2019-11-15 | 河海大学 | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 |
CN111562613A (zh) * | 2020-04-15 | 2020-08-21 | 东华理工大学 | 薄储层或互储层模型地震波频变反射系数分析方法 |
CN111597763A (zh) * | 2020-04-09 | 2020-08-28 | 东华理工大学 | 十字型多尺度流孔隙介质全频段弹性波频散衰减分析方法 |
CN111856571A (zh) * | 2020-06-15 | 2020-10-30 | 中国石油大学(北京) | 测井频段速度向地震频段速度转换方法、装置及设备 |
CN112799127A (zh) * | 2020-12-30 | 2021-05-14 | 成都理工大学 | 考虑裂缝孔隙岩石渗流非均匀差异的地震波频散和衰减数值计算方法 |
CN112946783A (zh) * | 2021-01-29 | 2021-06-11 | 中国石油大学(北京) | 一种水合物饱和度确定方法、装置及设备 |
CN113406697A (zh) * | 2020-03-17 | 2021-09-17 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN114114403A (zh) * | 2021-12-22 | 2022-03-01 | 东北石油大学 | 一种基于分数阶拉氏算子的各向异性衰减介质模拟方法 |
CN115774286A (zh) * | 2022-12-06 | 2023-03-10 | 上海电子信息职业技术学院 | 一种含黏弹性流体孔隙介质纵波速度预测方法及其应用 |
CN117055114A (zh) * | 2023-10-09 | 2023-11-14 | 中国石油大学(华东) | 储层沉积物游离气饱和度定量分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005119276A2 (en) * | 2004-05-27 | 2005-12-15 | Exxonmobil Upstream Research Company | Method for predicting lithology and porosity from seismic reflection data |
FR2949884A1 (fr) * | 2009-09-09 | 2011-03-11 | Ronald Tai | Procede de calcul deterministe des parametres d'ecoulements en milieu poreux et des proprietes de ces milieux poreux par la generation d'un milieu poreux artificiel equivalent |
CN102053259A (zh) * | 2009-10-28 | 2011-05-11 | 中国石油化工股份有限公司 | 一种油气藏识别方法 |
-
2011
- 2011-11-14 CN CN2011103601853A patent/CN102508296B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005119276A2 (en) * | 2004-05-27 | 2005-12-15 | Exxonmobil Upstream Research Company | Method for predicting lithology and porosity from seismic reflection data |
FR2949884A1 (fr) * | 2009-09-09 | 2011-03-11 | Ronald Tai | Procede de calcul deterministe des parametres d'ecoulements en milieu poreux et des proprietes de ces milieux poreux par la generation d'un milieu poreux artificiel equivalent |
CN102053259A (zh) * | 2009-10-28 | 2011-05-11 | 中国石油化工股份有限公司 | 一种油气藏识别方法 |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792574A (zh) * | 2012-10-30 | 2014-05-14 | 中国石油化工股份有限公司 | 一种储层频变含气性检测方法 |
CN103163553B (zh) * | 2013-02-01 | 2015-09-09 | 中国石油天然气股份有限公司 | 基于多重孔隙介质模型的地震烃类检测方法及装置 |
CN103163553A (zh) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | 基于多重孔隙介质模型的地震烃类检测方法及装置 |
CN104237936B (zh) * | 2013-06-06 | 2017-03-15 | 中国石油化工股份有限公司 | 一种油气检测的频变反演方法 |
CN104237936A (zh) * | 2013-06-06 | 2014-12-24 | 中国石油化工股份有限公司 | 一种油气检测的频变反演方法 |
CN104573150A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种针对富含有机质页岩的各向异性岩石物理建模方法 |
CN103616716B (zh) * | 2013-12-04 | 2016-05-11 | 中国石油天然气集团公司 | 一种双孔隙物理模型及其制作方法与应用 |
CN103616716A (zh) * | 2013-12-04 | 2014-03-05 | 中国石油天然气集团公司 | 一种双孔隙物理模型及其制作方法与应用 |
CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
CN104360383A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种地震波衰减预测的方法及系统 |
CN104360383B (zh) * | 2014-11-12 | 2017-03-22 | 中国石油大学(华东) | 一种地震波衰减预测的方法及系统 |
CN105759310B (zh) * | 2014-12-16 | 2018-05-08 | 中国石油化工股份有限公司 | 复杂非均质储层介质中地震波衰减及速度频散预测方法 |
CN105759310A (zh) * | 2014-12-16 | 2016-07-13 | 中国石油化工股份有限公司 | 复杂非均质储层介质中地震波衰减及速度频散预测方法 |
CN106324672A (zh) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及系统 |
CN106324672B (zh) * | 2015-07-03 | 2018-02-13 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及系统 |
CN105607123B (zh) * | 2016-01-15 | 2017-11-10 | 中国石油天然气股份有限公司 | 一种计算随机孔隙介质模型地震波特征信息的方法和装置 |
CN105607123A (zh) * | 2016-01-15 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种计算随机孔隙介质模型地震波特征信息的方法和装置 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107479090B (zh) * | 2016-06-08 | 2019-11-12 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107942375A (zh) * | 2017-11-17 | 2018-04-20 | 河海大学 | 一种基于声波方程的非线性优化隐式时空域有限差分数值模拟方法 |
CN107942375B (zh) * | 2017-11-17 | 2019-04-30 | 河海大学 | 一种基于声波方程的非线性优化隐式时空域有限差分数值模拟方法 |
CN108957542A (zh) * | 2018-07-10 | 2018-12-07 | 河海大学 | 一种建立地震波衰减岩石物理图板的方法 |
CN108957542B (zh) * | 2018-07-10 | 2020-04-03 | 河海大学 | 一种建立地震波衰减岩石物理图板的方法 |
CN109116420A (zh) * | 2018-10-16 | 2019-01-01 | 河海大学 | 一种含裂隙的孔隙介质纵波速度与衰减预测方法 |
CN109521090A (zh) * | 2018-10-18 | 2019-03-26 | 天津大学 | 一种激光无损表征薄膜杨氏模量的优化方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
CN110219643A (zh) * | 2019-06-05 | 2019-09-10 | 中国石油天然气集团有限公司 | 油藏含水饱和度确定方法及装置 |
CN110456412A (zh) * | 2019-07-09 | 2019-11-15 | 河海大学 | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 |
CN113406697B (zh) * | 2020-03-17 | 2023-09-26 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN113406697A (zh) * | 2020-03-17 | 2021-09-17 | 中国石油天然气股份有限公司 | 储层的横波品质因子获取方法和装置 |
CN111597763A (zh) * | 2020-04-09 | 2020-08-28 | 东华理工大学 | 十字型多尺度流孔隙介质全频段弹性波频散衰减分析方法 |
CN111597763B (zh) * | 2020-04-09 | 2023-06-27 | 东华理工大学 | 十字型多尺度流孔隙介质全频段弹性波频散衰减分析方法 |
CN111562613A (zh) * | 2020-04-15 | 2020-08-21 | 东华理工大学 | 薄储层或互储层模型地震波频变反射系数分析方法 |
CN111562613B (zh) * | 2020-04-15 | 2023-03-31 | 东华理工大学 | 薄储层或互储层模型地震波频变反射系数分析方法 |
CN111856571B (zh) * | 2020-06-15 | 2022-04-08 | 中国石油大学(北京) | 测井频段速度向地震频段速度转换方法、装置及设备 |
CN111856571A (zh) * | 2020-06-15 | 2020-10-30 | 中国石油大学(北京) | 测井频段速度向地震频段速度转换方法、装置及设备 |
CN112799127B (zh) * | 2020-12-30 | 2022-04-12 | 成都理工大学 | 考虑裂缝孔隙岩石渗流非均匀差异的地震波频散和衰减数值计算方法 |
CN112799127A (zh) * | 2020-12-30 | 2021-05-14 | 成都理工大学 | 考虑裂缝孔隙岩石渗流非均匀差异的地震波频散和衰减数值计算方法 |
CN112946783A (zh) * | 2021-01-29 | 2021-06-11 | 中国石油大学(北京) | 一种水合物饱和度确定方法、装置及设备 |
CN112946783B (zh) * | 2021-01-29 | 2024-02-06 | 中国石油大学(北京) | 一种水合物饱和度确定方法、装置及设备 |
CN114114403A (zh) * | 2021-12-22 | 2022-03-01 | 东北石油大学 | 一种基于分数阶拉氏算子的各向异性衰减介质模拟方法 |
CN114114403B (zh) * | 2021-12-22 | 2023-06-27 | 东北石油大学 | 一种基于分数阶拉氏算子的各向异性衰减介质模拟方法 |
CN115774286A (zh) * | 2022-12-06 | 2023-03-10 | 上海电子信息职业技术学院 | 一种含黏弹性流体孔隙介质纵波速度预测方法及其应用 |
CN115774286B (zh) * | 2022-12-06 | 2023-09-05 | 上海电子信息职业技术学院 | 一种含黏弹性流体孔隙介质纵波速度预测方法及其应用 |
CN117055114A (zh) * | 2023-10-09 | 2023-11-14 | 中国石油大学(华东) | 储层沉积物游离气饱和度定量分析方法 |
CN117055114B (zh) * | 2023-10-09 | 2023-12-29 | 中国石油大学(华东) | 储层沉积物游离气饱和度定量分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102508296B (zh) | 2013-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508296B (zh) | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 | |
Kuster et al. | Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations | |
Li et al. | Permeability evolution of shale under anisotropic true triaxial stress conditions | |
CN103760081B (zh) | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 | |
CN103323352B (zh) | 天然气水合物沉积物动三轴力学-声学-电学同步测试的实验装置及方法 | |
CN103984027B (zh) | 基于椭球体双重孔隙模型的岩石纵波速度预测方法 | |
Candela et al. | Laboratory evidence for particle mobilization as a mechanism for permeability enhancement via dynamic stressing | |
Zimmerman | Compressibility of sandstones | |
Dvorkin et al. | Squirt flow in fully saturated rocks | |
CN103293563B (zh) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 | |
CN103412336B (zh) | 一种非均质油藏中岩石系统的纵波速度预测方法 | |
CN106154351A (zh) | 一种低孔渗储层渗透率的估算方法 | |
CN103256046A (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
Ishibashi et al. | GeoFlow: A novel model simulator for prediction of the 3‐D channeling flow in a rock fracture network | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
CN104573150A (zh) | 一种针对富含有机质页岩的各向异性岩石物理建模方法 | |
Sheng et al. | Extended finite element modeling of multi-scale flow in fractured shale gas reservoirs | |
King | Recent developments in seismic rock physics | |
CN105301642B (zh) | 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置 | |
CN116738794A (zh) | 孔裂隙介质的岩石物理数值模拟方法、装置、设备及介质 | |
Mogensen et al. | Numerical study of long-time growth of hydraulic fractures in a line drive | |
Du et al. | Migration mechanism of fine particles in aquifer during water injection | |
Yang et al. | A model to calculate fracture conductivity considering proppant transport and embedment depth | |
Dande et al. | Effect of Fluids on the Elastic Properties of 3D-Printed Anisotropic Rock Models | |
Popkov et al. | On one of precise solutions to integrate hydrocarbon mobilization rhythms, dense reservoir development stages and sand recovery geological evolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |