CN109582922B - 一种水合物饱和度的现场快速判别方法及处理终端 - Google Patents

一种水合物饱和度的现场快速判别方法及处理终端 Download PDF

Info

Publication number
CN109582922B
CN109582922B CN201811316046.9A CN201811316046A CN109582922B CN 109582922 B CN109582922 B CN 109582922B CN 201811316046 A CN201811316046 A CN 201811316046A CN 109582922 B CN109582922 B CN 109582922B
Authority
CN
China
Prior art keywords
hydrate
value
hydrate concentration
quick discrimination
concentration
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.)
Expired - Fee Related
Application number
CN201811316046.9A
Other languages
English (en)
Other versions
CN109582922A (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.)
Tianjin Branch Of China Petroleum Group Logging Co Ltd
Guangzhou Marine Geological Survey
Original Assignee
Tianjin Branch Of China Petroleum Group Logging Co Ltd
Guangzhou Marine Geological Survey
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 Tianjin Branch Of China Petroleum Group Logging Co Ltd, Guangzhou Marine Geological Survey filed Critical Tianjin Branch Of China Petroleum Group Logging Co Ltd
Priority to CN201811316046.9A priority Critical patent/CN109582922B/zh
Publication of CN109582922A publication Critical patent/CN109582922A/zh
Application granted granted Critical
Publication of CN109582922B publication Critical patent/CN109582922B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种水合物饱和度的现场快速判别方法及处理终端,所述方法包括依次进行的如下步骤:步骤1:预设多元函数,并根据所述多元函数在平面直角坐标系下建立对应的多元函数交会图,得到水合物饱和度快速判别图版;步骤2:获取由从天然气水合物钻探井位采集的测井资料中得到的深电阻率RT值和声波时差AC值,并将RT值和AC值投影到水合物饱和度快速判别图版中,得到投影位置,根据投影位置得到水合物饱和度的值。本发明具有以下优点:只需要从现场的测井资料中得到声波时差AC和深电阻率RT的值即可快速判别出水合物饱和度,并采用水合物饱和度快速判别图版,不需要专业化的软件进行数据处理,即可在现场第一时间判别出水合物饱和度,节约时间,效率极高。

Description

一种水合物饱和度的现场快速判别方法及处理终端
技术领域
本发明涉及水合物饱和度计算技术领域,具体是一种水合物饱和度的现场快速判别方法及处理终端。
背景技术
在天然气水合物勘探开发过程中,水合物饱和度的确定是必不可少的重要参数。而对水合物饱和度的确定,常规的方法是基于Archie公式来计算的,需要将现场测井资料由测井解释人员利用专业化的软件进行处理后才能获取,不能在作业现场第一时间确定,而使得在勘探开发施工现场,不能及时依据水合物饱和度计算调整和完善方案部署,浪费了大量的时间,导致整体的水合物勘探作业效率低下。针对这样的情况,需要有能够在现场对水合物饱和度快速判别的方法,提高效率。
发明内容
针对现有技术的不足,本发明的目的之一提供一种水合物饱和度的现场快速判别方法,其能够解决现场对水合物饱和度快速判别的问题;
本发明的目的之二提供一种处理终端,能够解决现场对水合物饱和度快速判别的问题。
实现本发明的目的之一的技术方案为:一种水合物饱和度的现场快速判别方法,包括依次进行的如下步骤:
步骤1:预设多元函数,并根据所述多元函数在平面直角坐标系下建立对应的多元函数交会图,得到水合物饱和度快速判别图版,所述多元函数为:
式中,RT表示深电阻率,RT表示深电阻率基线,AC表示声波时差,AC表示声波时差基线,Sh表示水合物饱和度;
步骤2:获取由从天然气水合物钻探井位的测井资料中得到的RT值和AC值,并将RT值和AC值投影到水合物饱和度快速判别图版中,得到投影位置,根据投影位置得到水合物饱和度的值。
进一步地,所述水合物饱和度快速判别图版为电子件或纸质件。
实现本发明的目的之二的技术方案为:一种处理终端,其包括,
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行水合物饱和度的现场快速判别方法的步骤。
本发明的有益效果为:本发明具有以下优点:
1)只需要从天然气水合物钻探井位的测井资料中得到声波时差AC和深电阻率RT的值即可快速判别出水合物饱和度,不需要其他数据;
2)采用平面直角坐标系下建立的水合物饱和度快速判别图版,不需要专业化的软件进行数据处理,即可在现场第一时间判别出水合物饱和度,节约时间,效率极高;
3)水合物饱和度快速判别图版使用简单,不需要专业知识,任何人员都可以广泛使用,适用人群范围更广。
附图说明
图1为本发明对3口井进行拟合系数的示意图;
图2为采用双参数饱和度模型计算出的水合物饱和度和由实验确定出的水合物饱和度的对比图;
图3为本发明一种处理终端的结构示意图;
图4为本发明较佳实施例的流程图;
图5为本发明的水合物饱和度快速判别图版的示意图;
图6为天然气水合物饱和度快速判别图版应用示例图。
具体实施方式
下面,结合附图以及具体实施方式,对本发明做进一步描述:
首先介绍本发明的相关基础技术:如图1和图2所示,为了避免Archie公式对水合物饱和度计算的不准确问题,提高水合物饱和度的计算准确度和计算精度,可采用深电阻率+声波时差的双参数水合物饱和度计算方法,采用双参数计算水合物饱和度是基于以下思路得到的:
在长期对水合物物理特征的研究过程中,以及众多的地球物理测井资料中,深电阻率和声波时差这两个参数受水合物影响最为明显,相对高的深电阻率和相对低声波时差是水合物储层典型的测井响应特征,因此可以利用深电阻率和声波时差这两个参数来计算出水合物饱和度,并因此建立了计算公式(1)的数学模型来计算水合物饱和度Sh,也即是建立了基于深电阻率+声波时差双参数的数学模型,称之为双参数饱和度模型,其公式(1)如下:
式中,RT表示深电阻率,可以通过深电阻率测井得到,RT表示深电阻率基线,是目的层段非水合物地层深电阻率,可以从深电阻率测井曲线上得到,AC表示声波时差,可以通过声波测井得到,AC表示声波时差基线,是目的层段非水合物地层声波时差,可以从声波测井曲线上得到,a、b为双参数饱和度模型的系数,均为常数。
并将公式(1)两边同时除以lgR,得到公式(2):
Sh/lgR=a+blg(lgΔt/lgR) (2)
公式(2)属于线性函数y=a+bx的形式,因此可以通过线性回归拟合出系数a和b的值,进而确定饱和度计算公式。拟合出系数a和b的值的具体过程如下:
首先,选取南海神狐海域某3口井的水合物饱和度实验数据进行拟合,本实施例选取南海神狐海域3口井的水合物饱和度实验数据,实际情况也可以选择其他数量井的水合物饱和度实验数据和其他含有水合物区域的井,但选取的井为需要待计算水合物饱和度目标区域所在的井,通常来说,选取井的数量越多,所拟合出的系数a和b越精确,其中深电阻率RT可使用随钻相位深电阻率或不随钻测得的深电阻率,声波时差AC可使用随钻声波时差或不随钻测得的声波时差,也即RT、RT、AC和AC是对应3口井实测得到的测井值,Sh为通过对应井的岩心实验数据,实测得到的测井值和岩心实验数据统称为饱和度岩心数据,其中饱和度岩心数据如图1(包括图1a)至图1c))所示,所述饱和度岩心数据至少包括Sh求取得到3口井对应的系数后进行求算数平均,最终得到系数a和b的值分别为0.2和2.6,从而确定a和b的值,因此得到公式(3):
对公式(3)进行验证,采用公式(3)重新计算出上述3口井的水合物饱和度,并与通过实验确定该3口井以及其它井的水合物饱和度进行比较,以此验证公式(3)的适应性,其结果如图2和表1所示,图2中,包络曲线表示采用公式(3)计算出的水合物饱和度,离散点表示通过实验确定出的水合物饱和度,可以看出采用基于双参数的公式(3)计算出的水合物饱和度与实验确定的水合物饱和度一致性非常好,可见公式(3)能够避免Archie公式对水合物饱和度计算的不准确问题,有效提高了水合物饱和度的计算精度,降低了地层含天然气的影响。
表1是分别采用公式(3)计算出的水合物饱和度、采用Archie公式计算出的水合物饱和度和通过实验确定出的水合物饱和度的数据对比表格,其中岩电参数选取参考已发表文献数据,从表格中可以看出,采用本发明的双参数饱和度模型计算结果与岩心实验结果吻合较好,相比常规阿尔奇公式计算方法具有较高的精度。
表1双参数模型计算方法适应性分析对比表
如图4至图6所示,综上,从以上的基础技术资料,可见采用公式(3)计算水合物饱和度是准确且可行的,由公式(3)得到公式(5):
从公式(5)可以看出,由于RT在待测的目标区域内的各井的值是一致的、AC在待测的目标区域内的各井的测井值也是一致的,因此,RT和AC均可视为常数,因此可以将公式(5)看成是一个多元函数,因此可利用公式(5)在平面直角坐标系下建立对应的多元函数交会图,如图5所示,其中RT为纵坐标,AC为横坐标,多元函数交会图中的曲线表示Sh的值;
如图6所示,其为天然气水合物饱和度快速判别图版的一个应用示例图。利用多元函数交会图可以现场快速判别出水合物饱和度的值,只需要从天然气水合物钻探井位采集得到的测井资料中获得RT的值和AC的值,而RT的值和AC的值是可以在天然气水合物钻探井位的测井资料得到的数据,例如,在XX井某深度地层声波时差AC是130μs/ft,深电阻率RT是5ohm.m,通过AC=130,RT=5投影到图5的多元函数交会图中,找出其投影位置即可得到水合物饱和度的值,例如图6中投影位置为图中的“圆点”处,可以直接读出该深度下的水合物饱和度介于50%-60%之间,而且可以进一步估读出该深度下的水合物饱和度约为53%,由于公式(3)计算出的水合物饱和度是准确的,因此可以在现场利用公式(5)建立的多元函数交会图快速判别水合物饱和度的值,并将利用公式(5)建立的多元函数交会图称之为水合物饱和度快速判别图版。
在实际使用时,为了操作方便,除了将水合物饱和度快速判别图版做成电子版,还可以提前建立水合物饱和度快速判别图版并做成纸质件,然后根据天然气水合物钻探井位采集获得的AC值和RT值快速判别出水合物饱和度Sh的值,以上所有的操作都可以在现场完成,且不需要专业化的软件进行数据处理即可得到,操作非常简单。
而为了能够提高判别水合物饱和度的精确值,可以提高多元函数交会图的分辨率,即在图中增多Sh的曲线,比如图5中只有6条Sh的曲线(其值分别为0.1、0.2、0.3、0.4、0.5和0.6),可以增加至51条Sh的曲线例如分别为0.01,0.02,……,0.1,0.11,0.12,……,0.2,0.21,……,0.6,从而有效提高了分辨率,得到的水合物饱和度也更加精确。
另外,如图3所示,本发明还涉及一种实现以上方法的实体装置的处理终端100,其包括,
存储器101,用于存储程序指令;
处理器102,用于运行所述程序指令,以执行基水合物饱和度的现场快速判别方法的步骤,具体步骤与前面所描述的水合物饱和度的现场快速判别方法的步骤相同,在此不再进行赘述。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。

Claims (3)

1.一种水合物饱和度的现场快速判别方法,其特征在于:包括依次进行的如下步骤:
步骤1:预设多元函数,并根据所述多元函数在平面直角坐标系下建立对应的多元函数交会图,得到水合物饱和度快速判别图版,所述多元函数为:
式中,RT表示深电阻率,RT表示深电阻率基线,是目的层段非水合物地层深电阻率,从深电阻率测井曲线上得到,AC表示声波时差,AC表示声波时差基线,是目的层段非水合物地层声波时差,从声波测井曲线上得到,Sh表示水合物饱和度;
步骤2:获取由天然气水合物钻探井位采集的测井资料中得到的RT值和AC值,并将RT值和AC值投影到水合物饱和度快速判别图版中,得到投影位置,根据投影位置得到水合物饱和度的值。
2.根据权利要求1所述的水合物饱和度的现场快速判别方法,其特征在于:所述水合物饱和度快速判别图版为电子件或纸质件。
3.一种处理终端,其特征在于:其包括,
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行如权利要求1-2任意一项所述的水合物饱和度的现场快速判别方法的方法的步骤。
CN201811316046.9A 2018-11-06 2018-11-06 一种水合物饱和度的现场快速判别方法及处理终端 Expired - Fee Related CN109582922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811316046.9A CN109582922B (zh) 2018-11-06 2018-11-06 一种水合物饱和度的现场快速判别方法及处理终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811316046.9A CN109582922B (zh) 2018-11-06 2018-11-06 一种水合物饱和度的现场快速判别方法及处理终端

Publications (2)

Publication Number Publication Date
CN109582922A CN109582922A (zh) 2019-04-05
CN109582922B true CN109582922B (zh) 2019-10-18

Family

ID=65921681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811316046.9A Expired - Fee Related CN109582922B (zh) 2018-11-06 2018-11-06 一种水合物饱和度的现场快速判别方法及处理终端

Country Status (1)

Country Link
CN (1) CN109582922B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110927256B (zh) * 2019-04-22 2022-01-11 广东石油化工学院 一种基于Wood波阻抗法的天然气水合物饱和度计算方法
CN110057715B (zh) * 2019-04-23 2020-05-19 青岛海洋地质研究所 一种实验与数值模拟过程中水合物饱和度的计算分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103527172A (zh) * 2013-10-16 2014-01-22 中国石油集团川庆钻探工程有限公司 可变岩电耦合指数含水饱和度计算方法
CN103867197A (zh) * 2014-04-04 2014-06-18 中国石油集团川庆钻探工程有限公司 复杂岩性天然气层声波时差判别法
CN106556687A (zh) * 2016-11-21 2017-04-05 中国石油大学(华东) 弱胶结未成岩水合物声学和饱和度同步测试装置及方法
CN107578343A (zh) * 2017-07-26 2018-01-12 中国石油天然气股份有限公司 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7983846B2 (en) * 2004-09-14 2011-07-19 The Trustees Of Columbia University In The City Of New York Systems and methods for determining in-situ gas hydrate saturation
CN105158796B (zh) * 2015-09-16 2017-09-01 中国石油天然气股份有限公司 确定toc含量的方法和装置
CN107762483B (zh) * 2017-11-02 2020-10-13 中国石油天然气集团公司 一种基于测井曲线的相关系数及包络面积的流体识别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103527172A (zh) * 2013-10-16 2014-01-22 中国石油集团川庆钻探工程有限公司 可变岩电耦合指数含水饱和度计算方法
CN103867197A (zh) * 2014-04-04 2014-06-18 中国石油集团川庆钻探工程有限公司 复杂岩性天然气层声波时差判别法
CN106556687A (zh) * 2016-11-21 2017-04-05 中国石油大学(华东) 弱胶结未成岩水合物声学和饱和度同步测试装置及方法
CN107578343A (zh) * 2017-07-26 2018-01-12 中国石油天然气股份有限公司 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Doug Murray.etc.."Formation evaluation of gas hydrate reservoirs".《SPWLA 46th Annual Logging Symposium》.2005, *
杨胜雄等."南海北部神狐海域天然气水合物成藏特征及主控因".《地学前缘 ( 中国地质大学 ( 北京 ); 北京大学 )》.2017,第24卷(第4期), *
汪成芳等."王场油田水淹层测井识别与定量评价方法".《江汉石油职工大学学报》.2014, *

Also Published As

Publication number Publication date
CN109582922A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
CN104453874B (zh) 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法
Narr et al. A method for estimating subsurface fracture density in core
Lu et al. A multiresolution approach to reservoir parameter estimation using wavelet analysis
CN108303752A (zh) 砂砾岩有效储层常规测井定量识别方法
CN104863574B (zh) 一种适用于致密砂岩储层的流体识别方法
CN109582922B (zh) 一种水合物饱和度的现场快速判别方法及处理终端
CA2869825C (en) Three-dimensional multi-modal core and geological modeling for optimal field development
CN108804728B (zh) 水平井地层储层分级分析方法及计算机可读存储介质
Jin et al. Elliptical fracture network modeling with validation in Datong Mine, China
CN109582923A (zh) 一种基于双参数计算水合物饱和度的方法及处理终端
CN112861890B (zh) 一种储层评价模型构建方法及储层识别方法
CN105804737A (zh) 一种基于迭代算法求解地层孔隙度的方法
Jia et al. A multidisciplinary approach in cover-collapse sinkhole analyses in the mantle karst from Guangzhou City (SE China)
CN104459790A (zh) 含油气性盆地有效储层的分析方法和装置
Trenhaile The morphology of a drumlin field
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
Lin et al. Study on the equivalence between gas hydrate digital rocks and hydrate rock physical models
CN114060015B (zh) 一种致密砂岩含气性的评价方法及装置
Li et al. Efficient fracture network characterization in enhanced geothermal reservoirs by the integration of microseismic and borehole logs data
CN108733875A (zh) 一种Zoophycos遗迹化石分形拓扑模型及其构建方法和应用
CN108548765A (zh) 一种变粘土骨架的孔隙度计算方法
CN111624652B (zh) 确定地层中地质参数的方法
CN111983707B (zh) 一种致密砂岩气层可动水识别方法及系统
CN107818185A (zh) 一种获取致密砂岩地层因素的方法
Shen et al. Streamline-Based Production-Data Integration in Fractured-Reservoir Characterization and Its Application

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191018

Termination date: 20211106