CN106842311A - 粘弹性流体因子直接反演方法 - Google Patents
粘弹性流体因子直接反演方法 Download PDFInfo
- Publication number
- CN106842311A CN106842311A CN201510883703.8A CN201510883703A CN106842311A CN 106842311 A CN106842311 A CN 106842311A CN 201510883703 A CN201510883703 A CN 201510883703A CN 106842311 A CN106842311 A CN 106842311A
- Authority
- CN
- China
- Prior art keywords
- factor
- viscoelastic
- wave
- fluid factor
- viscoelastic fluid
- 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.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000014509 gene expression Effects 0.000 claims abstract description 17
- 238000010008 shearing Methods 0.000 claims abstract description 13
- 239000011435 rock Substances 0.000 claims description 13
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 238000005562 fading Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 241001269238 Data Species 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000035945 sensitivity Effects 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/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
- 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
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
本发明公开了一种粘弹性流体因子直接反演方法，以地下介质的粘弹性性质为理论基础，从理论上推导用粘弹性地震纵波反射近似表达式，建立地震纵波反射系数与粘弹性流体因子、剪切模量和密度等参数的理论关系，利用弹性波阻抗反演方法，通过叠前不同角度道集资料的联合，直接反演粘弹性流体因子。实践证明使用本发明方法明显提高储层流体识别的成功率。
Description
技术领域：
该发明属于油气勘探地震资料处理领域,具体是一种粘弹性流体因子直接反演方法。
现有技术：
地震勘探关注的最终目标是地下岩石中的流体，通过不同类型流体对地震波响应的差异来对流体进行区分。实际上，许多介质参数都与流体有密切的联系，如纵波速度、密度、纵波阻抗等，但是这些参数对流体响应的敏感程度不同。常用的流体因子反演方法是利用叠前地震资料直接反演这些弹性流体因子。现阶段油气预测的难点是如何从粘弹性介质理论出发，从叠前地震资料中直接反演粘弹性流体因子，提高对流体识别的精度和可靠性。
发明内容
本发明的目的是提出一种从叠前地震资料中直接反演粘弹性流体因子，提高对流体识别的精度和可靠性的粘弹性流体因子直接反演方法。
本发明的技术方案是：以地下介质的粘弹性性质为理论基础，从理论上推导用粘弹性地震纵波反射近似表达式，建立地震纵波反射系数与粘弹性流体因子、剪切模量和密度等参数的理论关系，利用弹性波阻抗反演方法，通过叠前不同角度道集资料的联合，直接反演粘弹性流体因子。
本发明粘弹性流体因子直接反演方法进一步包括:
(1)粘弹性流体因子纵波反射系数方程建立：
弹性介质的用流体因子表示的反射系数方程
其中，
θ表示角度，为弹性介质中纵波反射系数；f,μ,ρ分别表示弹性介质的流体因子，剪切模量和密度；Δf,Δμ,Δρ分别表示弹性介质的流体因子，剪切模量和密度的差值；表示弹性介质中干岩石的纵横波速度比的平方；表示弹性介质中水饱和岩石的纵横波速度比的平方；Vp,Vs表示表示弹性介质中的纵波速度 和横波速度。
代入粘弹性流体因子表达式
fela＝fane-ΔfQ (2)
其中，fela表示弹性介质的流体因子，fane表示粘弹性介质的流体因子，表示弹性介质与粘弹性介质间流体因子的差值，ω为介质频率，ωr为参考频率，QP为纵波品质因子。
经过一定的数学推导，得到用粘弹性流体因子表示的纵波反射系数方程
其中
Rpp表示纵波反射系数；表示粘弹性介质中水饱和岩石的纵横波速度比的平方；β,α表示粘弹性介质中的复纵波速度和复横波速度；fane分别表示粘弹性介质的流体因子；Δfane分别表示粘弹性介质的流体因子的差值。
假设干岩石纵横波速度比不受衰减影响。将水饱和岩石的纵横波速度比代入反射系数表达式，得
忽略反射系数中的虚部项，则上式可以进一步简化为
该表达式即建立了地震纵波反射系数与粘弹性流体因子、剪切模量和密度等三个参数之间的理论关系。基于该表达式，即可利用叠前地震资料直接反演粘弹性流体因子。
(2)粘弹性波阻抗方程建立：
按照弹性波阻抗公式，粘弹性波阻抗QEI进一步表达为含粘弹性流体因子的粘弹性阻抗
其中，
(3)粘弹性波阻抗反演：
叠前地震粘弹性流体因子直接反演方法是以弹性阻抗反演方法为基本框架，对多个角度数据体进行反演。一般而言，QEI反演同样需要测井资料和地质模型作为约束，减少反演的不确定性。QEI反演的流程可以简单概括为以下几个部分：角度部分叠加道集提取，粘弹性弹性阻抗曲线计算，角度子波提取，粘弹性弹性阻抗体反演。
从叠前地震资料中得到至少三个不同角度的粘弹性弹性阻抗体，即：QEI(θ1)、QEI(θ2)、QEI(θ3)，利用这三个粘弹性弹性阻抗体便可提取粘弹性流体因子、剪切模量和密度。由于方程(6)的非线性，若直接计算，势必带来不少麻烦，因此将模型进行线性变换。将方程两边取对数，可以得到如下模型：
ln(QEI)＝a(θ)ln(fane)+b(θ)ln(μ)+c(θ)ln(ρ) (7)
为了得到ln(fane)、ln(μ)和ln(ρ)，需要三个不同角度的粘弹性阻抗体。将三个角度值分别代入理论模型表达式(7)可以得到如下的方程组：
方程(8)写成矩阵的形式为：
对于方程(9)的求解，由于角度已知，只需要知道三个相互独立的粘弹性阻抗数据体，用常规的求解方法便可以求得。
发明效果
本发明粘弹性流体因子直接反演方法，填补了现有技术的空白，经实践证明，可以大幅提高储层流体识别的精度和成功率。
附图说明
图1弹性流体因子剖面；
图2粘弹性流体因子剖面。
具体实施方式
下面通过具体实施例对本发明组进一步说明。
一种利用叠前地震资料直接反演粘弹性流体因子的粘弹性流体因子直接反演方法。具体包括：
(1)粘弹性流体因子纵波反射系数方程建立：
弹性介质的用流体因子表示的反射系数方程
其中，
θ表示角度，为弹性介质中纵波反射系数；f,μ,ρ分别表示弹性介质的流体因子，剪切模量和密度；Δf,Δμ,Δρ分别表示弹性介质的流体因子，剪切模量和密度的差值；表示弹性介质中干岩石的纵横波速度比的平方；表示弹性介质中水饱和岩石的纵横波速度比的平方；Vp,Vs表示表示弹性介质中的纵波速度和横波速度。
代入粘弹性流体因子表达式
fela＝fane-ΔfQ (2)
其中，fela表示弹性介质的流体因子，fane表示粘弹性介质的流体因子，表示弹性介质与粘弹性介质间流体因子的差值，ω为介质频率，ωr为参考频率，QP为纵波品质因子。
经过一定的数学推导，得到用粘弹性流体因子表示的纵波反射系数方程
其中
Rpp表示纵波反射系数；表示粘弹性介质中水饱和岩石的纵横波速度比 的平方；β,α表示粘弹性介质中的复纵波速度和复横波速度；fane分别表示粘弹性介质的流体因子；Δfane分别表示粘弹性介质的流体因子的差值。
假设干岩石纵横波速度比不受衰减影响。将水饱和岩石的纵横波速度比代入反射系数表达式，得
忽略反射系数中的虚部项，则上式可以进一步简化为
该表达式即建立了地震纵波反射系数与粘弹性流体因子、剪切模量和密度等三个参数之间的理论关系。基于该表达式，即可利用叠前地震资料直接反演粘弹性流体因子。
(2)粘弹性波阻抗方程建立：
按照弹性波阻抗公式，粘弹性波阻抗QEI进一步表达为含粘弹性流体因子的粘弹性阻抗
其中，
(3)粘弹性波阻抗反演：
叠前地震粘弹性流体因子直接反演方法是以弹性阻抗反演方法为基本框架，对多个角度数据体进行反演。一般而言，QEI反演同样需要测井资料和地质模型作为约束，减少反演的不确定性。QEI反演的流程可以简单概括为以下几个部分：角度部分叠加道集提取，粘弹性弹性阻抗曲线计算，角度子波提取，粘弹性弹性阻抗体反演。
从叠前地震资料中得到至少三个不同角度的粘弹性弹性阻抗体，即：QEI(θ1)、QEI(θ2)、QEI(θ3)，利用这三个粘弹性弹性阻抗体便可提取粘弹性流体因子、剪切模量和密度。由于方程(6)的非线性，若直接计算，势必带来不少麻烦，因此将模型进行线性变换。将方程两边取对数，可以得到如下模型：
ln(QEI)＝a(θ)ln(fane)+b(θ)ln(μ)+c(θ)ln(ρ) (7)
为了得到ln(fane)、ln(μ)和ln(ρ)，需要三个不同角度的粘弹性阻抗体。将三个角度值分别代入理论模型表达式(7)可以得到如下的方程组：
方程(8)写成矩阵的形式为：
对于方程(9)的求解，由于角度已知，只需要知道三个相互独立的粘弹性阻抗数据体，用常规的求解方法便可以求得。
利用本方法对埕岛地区馆陶组地层进行测试。图1和图2是过埕北255井的剖面反演得到的弹性参数流体因子和粘弹性流体因子剖面，可以看出，在图中所示油水同层位置粘弹性流体因子比弹性流体因子吻合效果更好。
Claims (2)
1.粘弹性流体因子直接反演方法，其特征是：以地下介质的粘弹性性质为基础，推导出粘弹性地震纵波反射近似表达式，建立地震纵波反射系数与粘弹性流体因子、剪切模量和密度的关系，利用弹性波阻抗反演方法，通过叠前不同角度道集资料的联合，直接反演粘弹性流体因子。
2.根据权利要求1所述的粘弹性流体因子直接反演方法，其特征是，
(1)粘弹性流体因子纵波反射系数方程建立包括：
弹性介质用流体因子表示的反射系数方程
其中，
代入粘弹性流体因子表达式
fela＝fane-ΔfQ (2)
其中，
经过数学推导，得到用粘弹性流体因子表示的纵波反射系数方程：
其中
假设干岩石纵横波速度比不受衰减影响，将饱和岩石的纵横波速度比代入反射系数表达式，得
忽略反射系数中的虚部项，则上式可以进一步简化为
该表达式即建立了地震纵波反射系数与粘弹性流体因子、剪切模量和密度三个参数之间的理论关系，基于该表达式，即可利用叠前地震资料直接反演粘弹性流体因子；
(2)粘弹性波阻抗方程建立包括：
按照弹性波阻抗公式，粘弹性波阻抗QEI进一步表达为含粘弹性流体因子的粘弹性阻抗：
其中，
(3)粘弹性波阻抗反演包括：
从叠前地震资料中得到至少三个不同角度的粘弹性弹性阻抗体，即：QEI(θ1)、QEI(θ2)、QEI(θ3)，利用这三个粘弹性弹性阻抗体提取粘弹性流体因子、剪切模量和密度，将方程(6)的非线性模型进行线性变换得到如下模型
ln(QEI)＝a(θ)ln(fane)+b(θ)ln(μ)+c(θ)ln(ρ) (7)
为了得到ln(fane)、ln(μ)和ln(ρ)，需要三个不同角度的粘弹性阻抗体，将三个角度值分别代入理论模型表达式(7)得到如下的方程组：
方程(8)写成矩阵的形式为：
求解方程(9)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510883703.8A CN106842311A (zh) | 2015-12-04 | 2015-12-04 | 粘弹性流体因子直接反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510883703.8A CN106842311A (zh) | 2015-12-04 | 2015-12-04 | 粘弹性流体因子直接反演方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106842311A true CN106842311A (zh) | 2017-06-13 |
Family
ID=59148986
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510883703.8A Pending CN106842311A (zh) | 2015-12-04 | 2015-12-04 | 粘弹性流体因子直接反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106842311A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107656315A (zh) * | 2017-08-23 | 2018-02-02 | 中国石油天然气股份有限公司 | 储层流体因子的确定方法和装置 |
| CN109471165A (zh) * | 2018-12-03 | 2019-03-15 | 中国石油化工股份有限公司 | 基于包含敏感岩性识别因子为变量的avo近似式叠前反演方法 |
| CN110501745A (zh) * | 2019-08-27 | 2019-11-26 | 中海石油(中国)有限公司上海分公司 | 一种固液双相剥离烃类检测方法 |
| CN111060964A (zh) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
| CN111366974A (zh) * | 2018-12-26 | 2020-07-03 | 中国石油天然气股份有限公司 | 基于频变固液解耦流体因子的叠前地震反演方法和装置 |
| CN112764106A (zh) * | 2021-01-28 | 2021-05-07 | 中海石油深海开发有限公司 | 一种流体类型的确定方法、装置、设备及存储介质 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104316966A (zh) * | 2014-11-12 | 2015-01-28 | 中国石油大学(华东) | 一种流体识别方法及系统 |
-
2015
- 2015-12-04 CN CN201510883703.8A patent/CN106842311A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104316966A (zh) * | 2014-11-12 | 2015-01-28 | 中国石油大学(华东) | 一种流体识别方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| 刘浩杰等: "《地层吸收衰减反演方法》", 30 September 2015, 北京：石油工业出版社 * |
| 郑静静等: "曲波域瞬时衰减能量属性在储层识别中的应用研究", 《石油物探》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107656315A (zh) * | 2017-08-23 | 2018-02-02 | 中国石油天然气股份有限公司 | 储层流体因子的确定方法和装置 |
| CN107656315B (zh) * | 2017-08-23 | 2019-09-06 | 中国石油天然气股份有限公司 | 储层流体因子的确定方法和装置 |
| CN111060964A (zh) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
| CN111060964B (zh) * | 2018-10-16 | 2022-05-10 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
| CN109471165A (zh) * | 2018-12-03 | 2019-03-15 | 中国石油化工股份有限公司 | 基于包含敏感岩性识别因子为变量的avo近似式叠前反演方法 |
| CN111366974A (zh) * | 2018-12-26 | 2020-07-03 | 中国石油天然气股份有限公司 | 基于频变固液解耦流体因子的叠前地震反演方法和装置 |
| CN111366974B (zh) * | 2018-12-26 | 2022-10-04 | 中国石油天然气股份有限公司 | 基于频变固液解耦流体因子的叠前地震反演方法和装置 |
| CN110501745A (zh) * | 2019-08-27 | 2019-11-26 | 中海石油(中国)有限公司上海分公司 | 一种固液双相剥离烃类检测方法 |
| CN112764106A (zh) * | 2021-01-28 | 2021-05-07 | 中海石油深海开发有限公司 | 一种流体类型的确定方法、装置、设备及存储介质 |
| CN112764106B (zh) * | 2021-01-28 | 2023-09-29 | 中海石油深海开发有限公司 | 一种流体类型的确定方法、装置、设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106842311A (zh) | 粘弹性流体因子直接反演方法 | |
| CN104453873B (zh) | 页岩油气经济有效层段的评价方法 | |
| Rutqvist et al. | Coupled reservoir-geomechanical analysis of CO2 injection and ground deformations at In Salah, Algeria | |
| US20170145793A1 (en) | Method For Modeling Stimulated Reservoir Properties Resulting From Hydraulic Fracturing In Naturally Fractured Reservoirs | |
| Shi et al. | Assessment of carbon dioxide plume behaviour within the storage reservoir and the lower caprock around the KB-502 injection well at In Salah | |
| CN104483703B (zh) | 一种井震联合古地层压力预测方法 | |
| CN104360389B (zh) | 一种致密砂岩储层岩石弹性模量计算方法 | |
| CN105759310B (zh) | 复杂非均质储层介质中地震波衰减及速度频散预测方法 | |
| CN104898161B (zh) | 一种基于测井响应模拟体的有效砂岩预测方法 | |
| WO2010065209A3 (en) | Multiple anisotropic parameter inversion for a tti earth model | |
| CN104635261B (zh) | 山前带重磁电震联合分体建模方法 | |
| CN104570127B (zh) | 一种利用地震波速度同时反演孔隙度和泥质含量的方法 | |
| CN102339325A (zh) | 一种分析离散裂缝性油藏流体流动的方法 | |
| CN104500017A (zh) | 一种优化水平井分段压裂位置的方法 | |
| CN105319585A (zh) | 一种利用薄层干涉振幅恢复识别油气藏的方法 | |
| CN106979006A (zh) | 地层压力的确定方法和装置 | |
| CN104570065B (zh) | 一种利用地震波阻抗定量反演孔隙度的方法 | |
| CN104597491A (zh) | 一种基于方位弹性阻抗差异的裂缝储层参数反演方法 | |
| CN106526679A (zh) | 一种利用旋回分析法恢复不整合剥蚀厚度的方法 | |
| CN102914799A (zh) | 非等效体波场正演模拟方法及装置 | |
| CN104614764A (zh) | 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法 | |
| Meckel | Shelf margin deltas: the key to big reserves | |
| CN104422955A (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
| CN104484573A (zh) | 一种地层的刚性系数的计算新方法 | |
| CN106886045A (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |