CN107167844A - 一种砂泥岩薄互层流体替换方法 - Google Patents
一种砂泥岩薄互层流体替换方法 Download PDFInfo
- Publication number
- CN107167844A CN107167844A CN201710398303.7A CN201710398303A CN107167844A CN 107167844 A CN107167844 A CN 107167844A CN 201710398303 A CN201710398303 A CN 201710398303A CN 107167844 A CN107167844 A CN 107167844A
- Authority
- CN
- China
- Prior art keywords
- fluid
- sand
- mrow
- mud interbed
- thin
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 172
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000006467 substitution reaction Methods 0.000 title claims abstract description 22
- 239000004576 sand Substances 0.000 claims abstract description 72
- 238000012986 modification Methods 0.000 claims abstract description 18
- 230000004048 modification Effects 0.000 claims abstract description 18
- 239000011435 rock Substances 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000000470 constituent Substances 0.000 claims description 10
- 238000010008 shearing Methods 0.000 claims description 5
- 230000007812 deficiency Effects 0.000 abstract description 4
- 230000008859 change Effects 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 25
- 239000004575 stone Substances 0.000 description 7
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000002637 fluid replacement therapy Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000006854 communication Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
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/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/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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/622—Velocity, density or impedance
- G01V2210/6224—Density
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种砂泥岩薄互层流体替换方法，它解决了现有技术中流体替换方法误差较大的不足，仅改变砂泥岩薄互层中砂层的流体类型，不改变砂泥岩薄互层中泥层的流体类型；并且既不需要已知砂泥岩薄互层中砂层和泥层所占的比例，也不需要对砂泥岩薄互层进行“升降尺度”处理，其技术方案为：根据已知砂泥岩薄互层的参数，利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩薄互层中纯净砂岩层流体替换后纵波柔度的改变量成正比的关系模型，计算砂泥岩薄互层流体替换后的弹性特征。
Description
技术领域
本发明涉及油气资源勘探技术领域，尤其涉及一种砂泥岩薄互层流体替换方法。
背景技术
流体替换是一种非常重要的岩石物理方法，通常用来定量预测将岩石孔隙中的流体类型 替换为另一种流体类型之后的岩石弹性特征差异，包括纵波速度、横波速度、密度以及各种 弹性模量等。利用这些岩石弹性特征差异可以校正泥浆侵入的测井曲线、分析流体类型对地 震振幅随偏移距变化(Amplitude variation with offset，AVO)的影响以及模拟四维地震(4D Seismic)的地震响应特征差异。
流体替换最常用的工具就是Gassmann方程。而Gassmann方程通常假设岩石为单一矿物 均匀多孔介质，这样才能保证弹性波在岩石传播过程中引起的孔隙压力可以在整个孔隙空间 得到充分均衡。而当Gassmann方程应用于砂泥岩薄互层时，这些假设条件就不能满足，会 导致错误的结果。这是因为砂泥岩薄互层中的泥岩会通过以下几种方法影响Gassmann方程 的应用效果：(1)由于泥岩或者泥质砂岩中泥质部分的低孔隙度和低渗透性，会阻止弹性波 在岩石传播过程中引起的孔隙压力无法得到均衡，(2)泥质砂岩中泥质部分违反了Gassmann 方程的单一矿物组成的假设(3)理想情况下的砂泥岩薄互层流体替换是只对薄互层中的砂岩 部分进行流体替换，而泥岩不进行流体替换，而现有的观测方法(包括地震、测井、VSP等) 分辨率无法达到区分砂泥岩薄互层的砂岩和泥岩部分。
针对泥质砂岩储层的流体替换方法有很多，例如：利用经验公式对Gassmann方程中的 干燥岩石骨架的弹性模量进行修正的方法；将岩石孔隙中充填的泥质作为岩石基质的一部分， 利用Voigt-Reuss-Hill平均方法计算岩石基质弹性模量，同时利用有效孔隙度用来代替总孔隙 度的方法。然而这些方法都不是砂泥岩薄互层提出的。一般砂泥岩薄互层的厚度在地震分辨 率以下，因此测量得到弹性特征都是多个砂层和泥层的平均值。如果直接应用Gassmann方 程进行流体替换，而不考虑砂泥岩薄互层薄互层的非均质性，就会出现误差。
常规的砂泥岩薄互层流体替换方法就是首先利用逆Backus平均对砂泥岩薄互层进行“降 尺度”，分别得到各个砂层和泥层的特征(孔隙度、等效弹性模量)，然后应用Gassmann方程 仅对各个砂层进行流体替换，最后再利用Backus平均方法对砂泥岩薄互层进行“降尺度”，即 对所有的砂层和泥层进行Backus平均，得到整个薄互层的弹性特征。当砂泥岩薄互层中砂层 和泥层的比例已知的情况下，这种“升降尺度”的方法是比较理想的。然而实际过程中，砂 泥岩薄互层中砂层和泥层的比例很难确定，因此可能导致较大的误差，特别是砂层的比例很 小的时候。
发明内容
为了克服现有技术的不足，本发明提供了一种砂泥岩薄互层流体替换方法，该方法仅仅 改变砂泥岩薄互层中砂层的流体类型，不改变砂泥岩薄互层中泥层的流体类型；并且既不需 要已知砂泥岩薄互层中砂层和泥层所占的比例，也不需要对砂泥岩薄互层进行“升降尺度”处 理。
本发明采用下述技术方案：
一种砂泥岩薄互层流体替换方法，根据已知砂泥岩薄互层的参数，利用砂泥岩薄互层流 体替换后纵波柔度的改变量与砂泥岩薄互层中纯净砂岩层流体替换后纵波柔度的改变量成正 比的关系模型，计算砂泥岩薄互层流体替换后的弹性特征。
步骤1：收集砂泥岩薄互层的纵波速度、密度、有效孔隙度和饱和度；
步骤2：根据砂泥岩薄互层的纵波速度VPsat1和密度ρsat1计算含有原始流体类型1的砂泥岩 薄互层的纵波模量Msat1和纵波柔度Csat1；
步骤3：根据纯净砂岩层的孔隙度φsand和含有原始流体类型1的纯净砂岩层的纵波速度 VPsand1和密度ρsand1，计算含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量Msand1和纵波柔度Csand1；
步骤4：根据岩石物理经验值选取砂泥岩薄互层中纯净砂岩层岩石基质的体积模量Ksand0和剪切模量μsand0计算砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0；
步骤5：根据原始孔隙流体类型1的饱和度、孔隙流体组成成分的体积模量计算孔隙流体 类型1的体积模量Kf1和密度ρf1；
步骤6：利用替换后的孔隙流体类型2的饱和度、孔隙流体组成成分的体积模量计算孔隙 流体类型2的体积模量Kf2和密度ρf2；
步骤7：利用Gassmann方程计算原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄 互层中纯净砂岩层的纵波模量Msand2和纵波柔度Csand2；
步骤8：根据砂泥岩薄互层的有效孔隙度φeff计算砂泥岩薄互层含有的原始孔隙流体类型1 被孔隙流体类型2替换后的纵波柔度Csat2和纵波模量Msat2；计算砂泥岩薄互层含有的原始孔隙 流体类型1被孔隙流体类型2替换后的密度ρsat2和纵波速度VPsat2。
进一步的，所述的步骤2中，假定砂泥岩薄互层原始流体类型1被流体类型2替换，其 相应参量下标分别以1、2加以区别；
进一步的，所述的步骤3中，假定砂泥岩薄互层中纯净砂岩层原始流体类型1被流体类 型2替换，其相应参量下标分别以1、2加以区别；选取砂泥岩薄互层中纯净砂岩层确定纯净 砂岩层的孔隙度φsand；
进一步的，所述的步骤4中，砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0：
进一步的，所述的步骤5中，孔隙流体类型1的体积模量Kf1和密度ρf1分别为：
ρf1＝ρhydro·(1-Sw1)+ρwater·Sw1 (7)
式中，Khydro和Kwater分别为碳氢化合物(油或气)、水的体积模量，ρhydro和ρwater分别为 碳氢化合物(油或气)、水的密度，Sw1为孔隙流体类型1的含水饱和度，1-Sw1为孔隙流体类型1的含碳氢化合物(油或气)饱和度。
进一步的，所述的步骤6中，假定原始孔隙流体类型1被孔隙流体类型2替换，孔隙流 体类型2的体积模量Kf2和密度ρf2分别为：
ρf2＝ρhydro·(1-Sw2)+ρwater·Sw2 (9)
式中，Sw2为为替换后的孔隙流体类型2的含水饱和度，1-Sw2为替换后的孔隙流体类型 2的含碳氢化合物(油或气)饱和度。
进一步的，所述的步骤7中，原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄 互层中纯净砂岩层的纵波模量Msand2和纵波柔度Csand2分别为：
进一步的，所述的步骤8中，利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩 薄互层中纯净砂岩层流体替换后纵波柔度的改变量成正比的关系模型，原始孔隙流体类型1 被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度Csat2和纵波模量Msat2分别为：
砂泥岩薄互层含有的原始孔隙流体类型1被孔隙流体类型2替换后的密度ρsat2和纵波速 度VPsat2分别为：
ρsat2＝ρsat1+φeff·(ρf2-ρf1) (14)
与现有技术相比，本发明的有益效果是：
(1)本发明利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩薄互层中纯净砂岩 层流体替换后纵波柔度的改变量成正比的关系模型，计算砂泥岩薄互层流体替换后的弹性特 征，具有普遍适用性，避免了常规的方法需要已知砂泥岩薄互层中砂层和泥层所占的比例的 缺陷；
(2)本发明提出的砂泥岩薄互层流体替换方法，直接对观测尺度的数据进行流体替换， 不需要对数据进行“升降尺度”处理，简单易用，弥补了常规方法在砂岩层所占比例很小的 情况下会出现很大的误差的不足。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实 施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1为砂泥岩薄互层流体替换的流程图；
图2为某油田含气砂泥岩薄互层的测井曲线；
图3为含有原始流体类型1的砂泥岩薄互层的纵波模量和纵波柔度；
图4为含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量和纵波柔度、纯净 砂岩层岩石基质的纵波模量；
图5为孔隙流体类型1的体积模量和密度；
图6为替换后的孔隙流体类型2的体积模量和密度；
图7为原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄互层中纯净砂岩层的纵 波模量和纵波柔度；
图8为原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度和纵波模 量；
图9为原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的密度和纵波速度。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指 明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的 相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申 请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图 包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其 指明存在特征、步骤、操作、器件、组件和/或它们的组合。
正如背景技术所介绍的，现有技术中存在常规替换方法误差较大的不足，为了解决如上 的技术问题，本申请提出了一种砂泥岩薄互层流体替换方法。
本申请的一种典型的实施方式中，如图1所示，提供了一种砂泥岩薄互层流体替换方法， 根据已知砂泥岩薄互层的参数，利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩薄 互层中纯净砂岩层流体替换后纵波柔度的改变量成正比的关系模型，计算砂泥岩薄互层流体 替换后的弹性特征。
步骤1：收集砂泥岩薄互层的纵波速度、密度、有效孔隙度和饱和度；
步骤2：根据砂泥岩薄互层的纵波速度VPsat1和密度ρsat1计算含有原始流体类型1的砂泥岩 薄互层的纵波模量Msat1和纵波柔度Csat1；
假定砂泥岩薄互层原始流体类型1被流体类型2替换，其相应参量下标分别以1、2加以 区别；
步骤3：根据纯净砂岩层的孔隙度φsand和含有原始流体类型1的纯净砂岩层的纵波速度 VPsand1和密度ρsand1，计算含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量Msand1和纵波柔度Csand1；
假定砂泥岩薄互层中纯净砂岩层原始流体类型1被流体类型2替换，其相应参量下标分 别以1、2加以区别；选取砂泥岩薄互层中纯净砂岩层确定纯净砂岩层的孔隙度φsand；
步骤4：根据岩石物理经验值选取砂泥岩薄互层中纯净砂岩层岩石基质的体积模量Ksand0和剪切模量μsand0计算砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0；
砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0为：
步骤5：根据原始孔隙流体类型1的饱和度、孔隙流体组成成分的体积模量计算孔隙流体 类型1的体积模量Kf1和密度ρf1；
孔隙流体类型1的体积模量Kf1和密度ρf1分别为：
ρf1＝ρhydro·(1-Sw1)+ρwater·Sw1 (7)
式中，Khydro和Kwater分别为碳氢化合物(油或气)、水的体积模量，ρhydro和ρwater分别为 碳氢化合物(油或气)、水的密度，Sw1为孔隙流体类型1的含水饱和度，1-Sw1为孔隙流体类 型1的含碳氢化合物(油或气)饱和度。
步骤6：利用替换后的孔隙流体类型2的饱和度、孔隙流体组成成分的体积模量计算孔隙 流体类型2的体积模量Kf2和密度ρf2；
假定原始孔隙流体类型1被孔隙流体类型2替换，孔隙流体类型2的体积模量Kf2和密度 ρf2分别为：
ρf2＝ρhydro·(1-Sw2)+ρwater·Sw2 (9)
式中，Sw2为为替换后的孔隙流体类型2的含水饱和度，1-Sw2为替换后的孔隙流体类型2 的含碳氢化合物(油或气)饱和度。
步骤7：利用Gassmann方程计算原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄 互层中纯净砂岩层的纵波模量Msand2和纵波柔度Csand2；
原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄互层中纯净砂岩层的纵波模量 Msand2和纵波柔度Csand2分别为：
步骤8：利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩薄互层中纯净砂岩层流 体替换后纵波柔度的改变量成正比的关系模型，根据砂泥岩薄互层的有效孔隙度φeff计算砂泥 岩薄互层含有的原始孔隙流体类型1被孔隙流体类型2替换后的纵波柔度Csat2和纵波模量 Msat2；计算砂泥岩薄互层含有的原始孔隙流体类型1被孔隙流体类型2替换后的密度ρsat2和纵 波速度VPsat2。
原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度Csat2和纵波模量 Msat2分别为：
砂泥岩薄互层含有的原始孔隙流体类型1被孔隙流体类型2替换后的密度ρsat2和纵波速 度VPsat2分别为：
ρsat2＝ρsat1+φeff·(ρf2-ρf1) (14)
某油田含气砂泥岩薄互层的测井曲线如图2所示，包括泥质含量、有效孔隙度、原始孔 隙流体类型1含水饱和度、孔隙流体类型1含水饱和度砂泥岩薄互层的密度、纵波速度等， 这些数据为实施步骤1的输入数据。
含有原始流体类型1的砂泥岩薄互层饱和岩石的纵波模量和纵波柔度如图3所示，它们 是由砂泥岩薄互层的纵波速度和密度根据步骤2中的式(1)和(2)计算得到。
含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量和纵波柔度、纯净砂岩层 岩石基质的纵波模量如图4所示，它们是由含有原始流体类型1的纯净砂岩层的纵波速度和 密度根据步骤3中的式(3)和(4)、步骤4中的式(5)计算得到。这里纯净砂岩层的孔隙 度、纵波速度和密度分别取0.13，3340m/s和2.43g/cm3，纯净砂岩层岩石基质的体积模量和 剪切模量分别取36GPa和45GPa。
孔隙流体类型1的体积模量和密度如图5所示，它们是由原始孔隙流体类型1的饱和度、 孔隙流体组成成分的体积模量等根据步骤5中的式(6)和(7)计算得到，这里孔隙流体组 成成分为水和气，孔隙流体类型1的含水饱和度即为图2中的含水饱和度，水的体积模量和 密度分别取2.29GPa和1.0g/cm3，气的体积模量和密度分别取0.0208GPa和0.00001g/cm3。
替换后的孔隙流体类型2的体积模量和密度如图6所示，它们是由孔隙流体类型2的饱 和度、孔隙流体组成成分的体积模量等根据步骤6中的式(8)和(9)计算得到，这里孔隙 流体组成成分为水，孔隙流体类型2的含水饱和度即为100％，水的体积模量和密度分别取 2.29GPa和1.0g/cm3。
原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄互层中纯净砂岩层的纵波模量 和纵波柔度如图7所示，它们是由原始孔隙流体类型1的体积模量、替换后的孔隙流体类型 2的体积模量、含原始孔隙流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量、纯净砂岩 层岩石基质纵波模量等根据计算步骤7中的式(10)和(11)得到。
原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度和纵波模量如图 8所示，它们是由原始孔隙流体类型1被孔隙流体类型2替换前后的砂泥岩薄互层中纯净砂 岩层的纵波模量和纵波柔度、纯净砂岩层孔隙度、砂泥岩薄互层的有效孔隙度等利用新发明 的砂泥岩薄互层流体替换后的纵波柔度改变量与砂泥岩薄互层中纯净砂岩层流体替换后的纵 波柔度改变量成正比的关系模型，根据步骤8中(12)和(13)计算得到。
原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的密度和纵波速度如图9所 示，它们是由原始孔隙流体类型1的密度、孔隙流体类型2的密度、含孔隙流体类型1砂泥 岩薄互层的密度、含孔隙流体类型2砂泥岩薄互层的纵波模量等根据步骤8中(14)和(15) 计算得到。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员 来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等 同替换、改进等，均应包含在本申请的保护范围之内。
Claims (10)
1.一种砂泥岩薄互层流体替换方法，其特征在于，根据已知砂泥岩薄互层的参数，利用砂泥岩薄互层流体替换后纵波柔度的改变量与砂泥岩薄互层中纯净砂岩层流体替换后纵波柔度的改变量成正比的关系模型，计算砂泥岩薄互层流体替换后的弹性特征。
2.根据权利要求1所述的一种砂泥岩薄互层流体替换方法，其特征在于，包括以下步骤：
步骤1：根据已知砂泥岩薄互层的纵波速度VPsat1和密度ρsat1，计算含有原始流体类型1的砂泥岩薄互层的纵波模量Msat1和纵波柔度Csat1；
步骤2：根据纯净砂岩层的孔隙度φsand和含有原始流体类型1的纯净砂岩层的纵波速度VPsand1和密度ρsand1，计算含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量Msand1和纵波柔度Csand1；
步骤3：利用Gassmann方程计算原始孔隙流体类型1被孔隙流体类型2替换后的砂泥岩薄互层中纯净砂岩层的纵波模量Msand2和纵波柔度Csand2；
步骤4：利用所述的关系模型，根据砂泥岩薄互层的有效孔隙度φeff，计算原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度Csat2、纵波模量Msat2、纵波速度VPsat2和密度ρsat2。
3.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤1中原始流体类型1的砂泥岩薄互层的纵波模量Msat1为砂泥岩薄互层密度ρsat1与其纵波速度VPsat1平方的乘积；原始流体类型1的砂泥岩薄互层的纵波柔度Csat1为其纵波模量Msat1的倒数。
4.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤2中原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波模量Msand1为含有原始流体类型1的纯净砂岩层的密度ρsand1与其纵波速度VPsand1平方的乘积；含有原始流体类型1的砂泥岩薄互层中纯净砂岩层的纵波柔度Csand1为其纵波模量Msand1的倒数。
5.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，还包括：选取砂泥岩薄互层中纯净砂岩层岩石基质的体积模量Ksand0和剪切模量μsand0，计算砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0。
6.根据权利要求5所述的一种砂泥岩薄互层流体替换方法，其特征在于，砂泥岩薄互层中纯净砂岩层岩石基质的纵波模量Msand0为其体积模量Ksand0与4/3倍剪切模量μsand0之和。
7.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤2中还包括：根据原始孔隙流体类型1的饱和度、孔隙流体组成成分的体积模量，计算孔隙流体类型1的体积模量Kf1和密度ρf1。
8.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤3中利用替换后的孔隙流体类型2的饱和度、孔隙流体组成成分的体积模量，计算孔隙流体类型2的体积模量Kf2和密度ρf2。
9.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤4中原始孔隙流体类型1被孔隙流体类型2替换后砂泥岩薄互层的纵波柔度Csat2和纵波模量Msat2分别为：
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>.</mo>
</mrow>
10.根据权利要求2所述的一种砂泥岩薄互层流体替换方法，其特征在于，所述的步骤4中砂泥岩薄互层含有的原始孔隙流体类型1被孔隙流体类型2替换后的密度ρsat2和纵波速度VPsat2分别为：
ρsat2＝ρsat1+φeff·(ρf2-ρf1)
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>P</mi>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<msub>
<mi>M</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
</msqrt>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710398303.7A CN107167844B (zh) | 2017-05-31 | 2017-05-31 | 一种砂泥岩薄互层流体替换方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710398303.7A CN107167844B (zh) | 2017-05-31 | 2017-05-31 | 一种砂泥岩薄互层流体替换方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107167844A true CN107167844A (zh) | 2017-09-15 |
| CN107167844B CN107167844B (zh) | 2018-06-19 |
Family
ID=59821274
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710398303.7A Active CN107167844B (zh) | 2017-05-31 | 2017-05-31 | 一种砂泥岩薄互层流体替换方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107167844B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109115987A (zh) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102288993A (zh) * | 2011-05-16 | 2011-12-21 | 中国石油大学(北京) | 基于砂岩油藏叠前地震数据的流体替换方法 |
| US20120232861A1 (en) * | 2009-11-30 | 2012-09-13 | Pengbo Lu | Adaptive Newton's Method For Reservoir Simulation |
| CN104570129A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气股份有限公司 | 碳酸盐岩岩石储层横波速度预测方法 |
| US20150240624A1 (en) * | 2013-08-16 | 2015-08-27 | Landmark Graphics Corporation | Determining reserve estimates for a reservoir |
| CN105487122A (zh) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | 用于确定储层流体识别因子的敏感性的方法及系统 |
| CN105974474A (zh) * | 2016-07-07 | 2016-09-28 | 中国石油大学(华东) | 一种基于孔隙替换的储层孔隙类型评价的方法 |
-
2017
- 2017-05-31 CN CN201710398303.7A patent/CN107167844B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120232861A1 (en) * | 2009-11-30 | 2012-09-13 | Pengbo Lu | Adaptive Newton's Method For Reservoir Simulation |
| CN102288993A (zh) * | 2011-05-16 | 2011-12-21 | 中国石油大学(北京) | 基于砂岩油藏叠前地震数据的流体替换方法 |
| US20150240624A1 (en) * | 2013-08-16 | 2015-08-27 | Landmark Graphics Corporation | Determining reserve estimates for a reservoir |
| CN105487122A (zh) * | 2014-09-18 | 2016-04-13 | 中国石油化工股份有限公司 | 用于确定储层流体识别因子的敏感性的方法及系统 |
| CN104570129A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气股份有限公司 | 碳酸盐岩岩石储层横波速度预测方法 |
| CN105974474A (zh) * | 2016-07-07 | 2016-09-28 | 中国石油大学(华东) | 一种基于孔隙替换的储层孔隙类型评价的方法 |
Non-Patent Citations (5)
| Title |
|---|
| 叶树刚: "含气煤层纵波反射波AVO响应特征研究", 《煤炭科学技术》 * |
| 李中元 等: "基于流体替换的横波速度反演", 《石油天然气学报》 * |
| 李景叶: "基于岩石物理的薄互层流体替换方法研究", 《石油学报》 * |
| 杨士明 等: "叠前弹性参数反演技术在识别白云岩储层中的应用", 《石油天然气学报》 * |
| 林凯等: "基于Gassmann方程的鲕滩储层流体替换模拟技术及其应用", 《石油物探》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109115987A (zh) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107167844B (zh) | 2018-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
| CN104712329B (zh) | 一种泥页岩油气饱和度的计算模型 | |
| CN104950331B (zh) | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 | |
| CN102854531B (zh) | 用于储层烃类检测的多尺度岩石物理图版方法及装置 | |
| Qadri et al. | Three-dimensional petrophysical modelling and volumetric analysis to model the reservoir potential of the Kupe Field, Taranaki Basin, New Zealand | |
| CN104570129A (zh) | 碳酸盐岩岩石储层横波速度预测方法 | |
| CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
| CN106125135B (zh) | 基于岩石物理模型的含气砂岩储层地震响应数值模拟方法 | |
| CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
| CN104007485A (zh) | 一种建立复杂孔隙泥页岩岩石物理模型方法 | |
| CN106337682A (zh) | 页岩储层脆性评价方法 | |
| Olea et al. | A methodology for the assessment of unconventional (continuous) resources with an application to the Greater Natural Buttes Gas Field, Utah | |
| CN105974474A (zh) | 一种基于孔隙替换的储层孔隙类型评价的方法 | |
| CN105840189A (zh) | 一种页岩气岩石物理与测井方法 | |
| CN105425280A (zh) | 一种矿物模量和孔隙结构预测方法 | |
| CN107194104A (zh) | 一种采用测井曲线计算页岩储层中无机质孔隙度、总孔隙度以及组分的方法 | |
| CN105931125A (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
| Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
| Acheampong et al. | Quantitative interpretation of time‐lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time‐lapse velocity responses | |
| CN111273341B (zh) | 根据裂缝空间分布的含裂缝储层岩石物理建模方法 | |
| CN109655936B (zh) | 一种碎屑岩岩性替换的弹性参数计算方法及系统 | |
| CN107167844A (zh) | 一种砂泥岩薄互层流体替换方法 | |
| CN110954942A (zh) | 粘土含量识别的图版制作方法及系统 | |
| CN112946739B (zh) | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 | |
| CN104792684B (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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liang Jie Inventor after: Zhang Jiajia Inventor after: Li Hongbing Inventor after: Zhao Wenzhi Inventor after: Yin Xingyao Inventor after: Zhang Guangzhi Inventor after: Zong Zhaoyun Inventor before: Zhang Jiajia Inventor before: Li Hongbing Inventor before: Zhao Wenzhi Inventor before: Yin Xingyao Inventor before: Zhang Guangzhi Inventor before: Zong Zhaoyun |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210202 Address after: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong. Patentee after: CHINA University OF PETROLEUM (EAST CHINA) Patentee after: PetroChina Company Limited Address before: 266580 No. 66 Changjiang West Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: China University of Petroleum (East China) |