CN104198516B - 应用核磁共振t2几何平均值预测横波速度的方法 - Google Patents
应用核磁共振t2几何平均值预测横波速度的方法 Download PDFInfo
- Publication number
- CN104198516B CN104198516B CN201410453724.1A CN201410453724A CN104198516B CN 104198516 B CN104198516 B CN 104198516B CN 201410453724 A CN201410453724 A CN 201410453724A CN 104198516 B CN104198516 B CN 104198516B
- Authority
- CN
- China
- Prior art keywords
- msub
- transverse wave
- geometric mean
- wave velocity
- mrow
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 26
- 239000011435 rock Substances 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 239000011148 porous material Substances 0.000 abstract description 12
- 230000000694 effects Effects 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 235000019738 Limestone Nutrition 0.000 description 2
- 102000011990 Sirtuin Human genes 0.000 description 2
- 108050002485 Sirtuin Proteins 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及预测横波速度领域,特别涉及一种应用核磁共振T2几何平均值预测横波波速的方法,是一种利用岩心刻度测井建立T2几何平均值与横波波速关系模型预测横波速度的方法。本发明的预测方法不仅包含孔隙度的影响,而且还考虑了孔隙尺度、比表面等孔隙结构的影响,因此,该方法岩石物理基础更完善,物理意义更明确,预测精度更高。本发明模型建立简单,预测效果能达到地球物理应用对横波速度应用的要求,相对理论模型预测法更简单实用。
Description
技术领域
本发明涉及预测横波速度领域,特别涉及一种应用核磁共振T2几何平均值预测横波速度的方法。
背景技术
横波是传播方向与质点震动方向垂直的波,横波速度是用于研究地层性质重要参数,但其测量比较困难,因此,借助其他手段进行地层横波速度预测具有重要意义。目前横波速度预测方法分为两大类:第一种是回归经验公式预测法,通过对大量的横波速度数据和其他数据进行统计分析回归得到横波速度经验关系式。比较常用的经验公式有,Castagna等给出了不同条件下的横波速度经验公式;han等,基于75块砂岩岩样在不同压力下的测量结果,给出了纵波速度和横波速度经验公式。Castagna和Backus对横波速度经验公式进行了总结,分别给出饱含水石灰岩、白云岩、砂岩以及泥岩的横波速度回归公式。经验公式预测法大部分涉及纵横波速、孔隙度等之间的关系,未考虑孔隙结构,缺乏明确的物理意义,往往精度有限。第二种是理论模型预测法,1995年Xu和White联合K-T模型,微分等效介质理论,以及Gassmann理论,提出了砂、泥岩混合岩石的理论模型。该模型考虑了孔隙形状对纵波、横波速度的影响,具有较高的精度。为了减小Xu-White模型参数提取误差对预测精度的影响,乔悦东、高云峰等在理论研究的基础上,通过引进虚拟孔隙度变量和非线性优化算法,形成了一种基于Xu-White模型的优化测井横波速度预测技术。孙福利等采用Biot-Gassmann低频速度模型,以及Pride公式建立起基质弹性模量与骨架弹性模量关系,通过迭代方式计算出合适的固结系数,进而得到横波速度大小。理论模型具有较好的物理基础和明确的物理意义,但理论模型通常比较复杂且计算量较大。
现有技术的缺点:
(1)预测地层横波速度经验公式,如利用纵波速度、密度、泥质含量等参数来预测横波速度,它们共同的切入点是基于岩石孔隙度,这是因为这些参数均与孔隙度有关,需要指出的是孔隙度只是一个相对体积概念(孔隙体积占岩石总体积的比例),它不存在孔隙形态、尺度的影响,但是横波是一种剪切波,横波速度不仅要考虑孔隙度,还要考虑孔隙尺度、比表面等孔隙结构参数影响。因此,目前常见的利用纵波速度、泥质含量、密度等预测横波速度方法,虽然计算效率较快,但其岩石物理基础存在明显缺欠,精度较差。
(2)理论模型计算考虑了孔隙形状等对横波速度的影响,计算精度较高,但是由于计算参 数较多或者要求给定基质矿物弹性参数等,数据不充分、结果多解性较强、横波速度预测计算量较大、计算效率低。
发明内容
本发明针对现有技术的不足,提供一种计算简单、预测准确的,利用岩心刻度测井建立T2几何平均值与横波速度关系模型预测横波速度的方法。
本发明的技术方案是:
一种应用核磁共振T2几何平均预测值横波速度方法,所述的方法如下:
步骤1:T2几何平均值的获取
对工区的岩石样品在进行核磁共振测量,得到核磁共振T2谱;通过得到的核磁共振T2谱,得到T2分布点和对应孔隙度分量,结合公式(1)得到T2几何平均值;
式中:T2gm为几何平均值,φ为岩石总孔隙度,T2i为各T2分布点,φi为各T2分布点对应的孔隙度分量;
其中,核磁共振测量方法为:使用主频为2MHz的核磁共振采集仪器,采用回波间隔0.37ms,等待时间6000ms,累加次数128次,回波个数4096的采集参数进行核磁共振信号采集;采用SIRT反演方法,迭代次数5000进行T2谱反演。
步骤2:横波速度测量
对工区的岩石样品进行横波速度测量,得到相应的横波速度;
其中,横波速度测量方法为:使用计算机脉冲信号发生接收器产生和接受脉冲信号;应用Olympus横波探头测量横波信号;使用混合信号示波器采集横波波形,根据横波波形判断横波到时,并由此计算横波速度。
步骤3:建立T2几何平均值和横波速度的关系模型
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=AT2gm -a (2)
其中,vs为横波速度,T2gm为几何平均值,A为系数,a为指数;
其中,
(3)式是在理想的岩石物理体积模型子下得出的,A为化简后的系数。(2)式是真实地层中T2几何平均值与横波速度关系式,系数A以及T2的指数大小a,与孔隙形状因子、 比面、岩性等有关,需要根据不同地区实验数据拟合得出。在真实地层中,孔隙结构一般较复杂,故T2的指数的绝对值一般小1。
步骤4:横波速度预测
将步骤3建立的T2几何平均值和相应的横波速度关系模型加载到测井解释软件平台中,带入计算得到的T2几何平均值,得到预测横波速度曲线。
本发明的有益效果是:
利用核磁共振T2几何平均值预测横波速度,其优势在于它不仅包含孔隙度的影响,而且还考虑了孔隙尺度、比表面等孔隙结构的影响,因此,该方法岩石物理基础更完善,物理意义更明确,预测精度更高。本发明模型建立简单,预测效果能达到地球物理应用对横波速度应用的要求,相对理论模型预测法更简单实用。
利用本发明的方法,根据核磁共振进行横波速度预测处理,并将其同偶极横波测井实测横波速度进行对比,结果如图5,显示利用本方法预测的横波速度同实际情况非常吻合。
附图说明
附图1为本发明具体实施例1的Vs-T2gm关系曲线;
附图2为本发明具体实施例2的Vs-T2gm关系曲线;
附图3为本发明具体实施例3的Vs-T2gm关系曲线;
附图4为本发明具体实施例4的Vs-T2gm关系曲线;
附图5为本发明所采取的华北某井砂岩段横波速度预测成果图。
具体实施方式
本发明的以华北某砂岩为例,具体实施方式如下:
实施例1:砂岩
步骤1:T2几何平均值的获取
对工区的岩石样品在进行核磁共振测量,得到核磁共振T2谱;通过得到的核磁共振T2谱,得到T2分布点和对应孔隙度分量,结合公式(1)得到T2几何平均值;
式中:T2gm为几何平均值,φ为岩石总孔隙度,T2i为各T2分布点,φi为各T2分布点对应的孔隙度分量;
其中,核磁共振测量方法为:使用主频为2MHz的核磁共振采集仪器,采用回波间隔 0.37ms,等待时间6000ms,累加次数128次,回波个数4096的采集参数进行核磁共振信号采集;采用SIRT反演方法,迭代次数5000进行T2谱反演。
步骤2:横波速度测量
对工区的岩石样品进行横波速度测量,得到相应的横波速度;
其中,横波速度测量方法为:使用计算机脉冲信号发生接收器产生和接受脉冲信号;应用Olympus横波探头测量横波信号;使用混合信号示波器采集横波波形,根据横波波形判断横波到时,并由此计算横波速度。
实验数据如下:
岩样编号 | T2几何平均值(ms) | 横波速度(m/s) |
1 | 2.6522 | 2339.241 |
2 | 2.66579 | 2271.774 |
3 | 7.1105 | 2012.291 |
4 | 2.74471 | 2412.414 |
5 | 17.57908 | 1927.041 |
6 | 28.788 | 1894.083 |
7 | 7.91401 | 2132.749 |
8 | 19.48719 | 1868.297 |
9 | 3.42875 | 2547.712 |
10 | 2.35947 | 2487.5 |
11 | 7.59156 | 2067.021 |
12 | 2.45175 | 2538.983 |
13 | 2.87476 | 2670.548 |
14 | 7.04889 | 2024.51 |
15 | 6.47566 | 2049.189 |
16 | 5.90399 | 2117.614 |
17 | 6.12803 | 2269.737 |
18 | 39.2503 | 1793.103 |
19 | 15.0718 | 2016.548 |
20 | 11.67429 | 2136.449 |
21 | 15.07156 | 1967.624 |
步骤3:建立T2几何平均值和横波速度的关系模型(如图1)
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=2796.6T2gm -0.125
其中,vs为横波速度,T2gm为几何平均值。
步骤4:横波速度预测
将步骤3建立的T2几何平均值和相应的横波速度关系模型加载到测井解释软件平台中,带入计算得到的T2几何平均值,得到预测横波速度曲线。
实施例2:深灰色油斑含砾不等粒砂岩:横波速度预测过程同实施例1。实验数据如下:
岩样编号 | T2几何平均值(ms) | 横波速度(m/s) |
1 | 11.58283 | 1948.01388 |
2 | 15.31018 | 1849.36783 |
3 | 6.93642 | 2039.66712 |
4 | 5.20351 | 2345.05009 |
5 | 21.73941 | 1740.66274 |
6 | 18.70473 | 1701.51255 |
7 | 20.03811 | 1639.06962 |
8 | 3.6627 | 2481.98277 |
9 | 10.17297 | 1851.66382 |
10 | 7.69804 | 1843.87433 |
11 | 15.21751 | 1743.67129 |
建立T2几何平均值和横波速度的关系模型(如图2)
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=3128.3T2gm -0.207
实施例3:灰岩:横波速度预测过程同实施例1。
实验数据如下:
岩样编号 | T2几何平均值(ms) | 横波速度(m/s) |
1 | 2.603 | 2796.2963 |
2 | 2.26 | 2665.83072 |
3 | 2.356 | 2766.76602 |
4 | 52.475 | 1420.04 |
5 | 18.488 | 1876.33588 |
6 | 10.361 | 2205.19 |
7 | 12.896 | 1942.18 |
8 | 2.494 | 2809.7166 |
9 | 2.55 | 3002.11864 |
10 | 1.921 | 3060.31746 |
建立T2几何平均值和横波速度的关系模型(如图3)
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=3437.8T2gm -0.215
实施例4:云岩:横波速度预测过程同实施例1。
实验数据如下:
岩样编号 | T2几何平均值(ms) | 横波速度(m/s) |
1 | 5.551 | 3126.82927 |
2 | 34.702 | 3059.88858 |
3 | 196.753 | 2658.31296 |
4 | 195.11 | 2582.0802 |
5 | 264.045 | 2549.25373 |
6 | 351.368 | 2552.49267 |
7 | 357.229 | 2646.30422 |
8 | 20.65 | 3113.74738 |
9 | 21.155 | 3418.81919 |
10 | 229.855 | 2692.69207 |
11 | 120.079 | 2783.06092 |
12 | 148.149 | 2657.45192 |
13 | 534.107 | 2514.62687 |
14 | 409.972 | 2502.1565 |
建立T2几何平均值和横波速度的关系模型(如图4)
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=3776T2gm -0.066
如图5所示,对四个层段的砂岩通过该方法进行横波预测,将计算的横波速度曲线与实际测得的横波速度曲线进行比较,两条曲线基本吻合,预测结果能够满足实际工作需要。在岩性变化较大或井径扩井的地方预测效果会受到影响。在缺少测井横波资料的地区,利用该技术能够提供储层岩性、流体识别等技术研究和应用所需要的测井横波速度数据。
需要注意的是,不同的岩性或相同岩性不同的地区建立的T2几何平均值与横波速度计算模型不同,如图1-图4。因此,对工区建立模型时,应划分不同的岩性,并结合不同的流体性质进行建立,以提高横波速度预测精度。另外,井孔扩径段或仪器高度偏心,核磁共振测得的回波信号会受到井眼泥浆信号的影响,造成NMR孔隙度大于地层真实孔隙度,此时NMR测井资料失真,利用T2几何平均值预测横波速度会受到影响,需要进行必要的修正。
Claims (1)
1.一种应用核磁共振T2几何平均预测值横波速度方法,所述的方法如下:
步骤1:T2几何平均值的获取
对工区的岩石样品在进行核磁共振测量,得到核磁共振T2谱;通过得到的核磁共振T2谱,得到T2分布点和对应孔隙度分量,结合公式(1)得到T2几何平均值;
式中:T2gm为几何平均值,φ为岩石总孔隙度,T2i为各T2分布点,φi为各T2分布点对应的孔隙度分量;
步骤2:横波速度测量
对工区的岩石样品进行横波速度测量,得到相应的横波速度;
步骤3:建立T2几何平均值和横波速度的关系模型
根据不同岩性和层段建立T2几何平均值和相应的横波速度关系模型:
vs=AT2gm -a (2)
其中,vs为横波波速,T2gm为几何平均值,A为系数,a为指数;
步骤4:横波速度预测
将步骤3建立的T2几何平均值和相应的横波速度关系模型加载到测井解释软件平台中,带入计算得到的T2几何平均值,得到预测横波速度曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410453724.1A CN104198516B (zh) | 2014-09-05 | 2014-09-05 | 应用核磁共振t2几何平均值预测横波速度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410453724.1A CN104198516B (zh) | 2014-09-05 | 2014-09-05 | 应用核磁共振t2几何平均值预测横波速度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104198516A CN104198516A (zh) | 2014-12-10 |
CN104198516B true CN104198516B (zh) | 2015-06-03 |
Family
ID=52083841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410453724.1A Active CN104198516B (zh) | 2014-09-05 | 2014-09-05 | 应用核磁共振t2几何平均值预测横波速度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104198516B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866160B (zh) * | 2016-06-16 | 2017-09-01 | 中国石油大学(华东) | 一种基于纵波约束的核磁共振t2截止值计算方法 |
CN107843611B (zh) * | 2016-09-20 | 2020-07-14 | 中国石油化工股份有限公司 | 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法 |
CN110057853B (zh) * | 2019-04-11 | 2022-02-18 | 中国石油大学(华东) | 一种基于低场核磁共振响应的岩石杨氏模量计算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5265016A (en) * | 1989-07-19 | 1993-11-23 | Amoco Corporation | Method of shear wave velocity estimation |
JP2001004606A (ja) * | 1999-06-25 | 2001-01-12 | Taisei Corp | せん断波速度の測定方法 |
CN102042011B (zh) * | 2010-10-13 | 2014-10-01 | 中国石油化工集团公司 | 利用常规测井资料构建伪核磁t2谱的方法 |
CN103424772B (zh) * | 2012-05-24 | 2016-09-14 | 中国石油化工股份有限公司 | 一种基于岩石物理的碳酸盐储层横波速度预测方法 |
CN103576196B (zh) * | 2013-10-28 | 2016-05-11 | 西北大学 | 一种随压力变化的孔隙介质横波速度预测方法 |
-
2014
- 2014-09-05 CN CN201410453724.1A patent/CN104198516B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104198516A (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108072902B (zh) | 一种碳酸盐岩岩石物理模型的建模方法 | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
CN102937721B (zh) | 利用初至波走时的有限频层析成像方法 | |
CN109490965B (zh) | 一种定量评价地层非均匀性的方法及装置 | |
CN103293563B (zh) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 | |
CN110348135B (zh) | 一种随钻声波测井评价地层渗透率的方法 | |
CN106855636A (zh) | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 | |
CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN103954999A (zh) | 一种适用于低孔隙度砂泥岩地层的横波速度预测新方法 | |
CN103728659A (zh) | 一种提高地下岩溶探测精度的方法 | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
CN104345346A (zh) | 一种获取裂缝宽度的方法 | |
CN105629308A (zh) | 相控非均质力学参数地应力方法 | |
CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN101285381B (zh) | 一种泄漏模式波反演软地层横波速度的方法 | |
WO2020215170A1 (zh) | 基于地震岩石物理实验分析的测井与地震速度匹配方法 | |
CN105158796A (zh) | 确定toc含量的方法和装置 | |
CN104198516B (zh) | 应用核磁共振t2几何平均值预测横波速度的方法 | |
CN111368245B (zh) | 一种考虑多影响因素的定量计算岩石Biot系数方法 | |
EP4089446A1 (en) | A method for seismic frequency resonance exploration technology | |
CN102914797B (zh) | 一种获得地层各向异性系数的方法及装置 | |
CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
CN104007466A (zh) | 一种基于纵波振幅的无井约束叠前反演实现的储层与流体预测方法 |
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 |