CN109001803B - 采用拟合微电位曲线提高地震分辨率的方法 - Google Patents

采用拟合微电位曲线提高地震分辨率的方法 Download PDF

Info

Publication number
CN109001803B
CN109001803B CN201811127695.4A CN201811127695A CN109001803B CN 109001803 B CN109001803 B CN 109001803B CN 201811127695 A CN201811127695 A CN 201811127695A CN 109001803 B CN109001803 B CN 109001803B
Authority
CN
China
Prior art keywords
micro
data
new
potential
seismic
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
Application number
CN201811127695.4A
Other languages
English (en)
Other versions
CN109001803A (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.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
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 China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201811127695.4A priority Critical patent/CN109001803B/zh
Publication of CN109001803A publication Critical patent/CN109001803A/zh
Application granted granted Critical
Publication of CN109001803B publication Critical patent/CN109001803B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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

本发明提供一种采用拟合微电位曲线提高地震分辨率的方法，包括：步骤1，对地震剖面进行精细地层解释；步骤2，根据相邻两口井同一层位、同一位置的微电位数据点RN1、RN2，得到新层面的数据体RN’的位置；步骤3，根据相邻两井的波阻抗（I1、I2），对同一层位、同一位置的微电位数据（RN1、RN2）进行变换，且按新数据点RN’距两井的距离（D1、D2）设定变形的权重；步骤4，采用步骤3中的计算方法，算出每个精细层面的全部微电极新数据RN’；步骤5，形成拟合微电位RN’新三维数据。该采用拟合微电位曲线提高地震分辨率的方法大大提高地震剖面的分辨率，对地震资料升级、更精细的地层解释提供了新的方法。

Description

采用拟合微电位曲线提高地震分辨率的方法
技术领域
本发明涉及地球物理技术领域，特别是涉及到一种采用拟合微电位曲线提高地震分辨率的方法。
背景技术
目前提高地震分辨率的方法，几乎全部是对波阻抗曲线进行合成分解，而不是把地震剖面的全套波阻抗数据，全部换成其它高精度测井曲线的数据，因此地震剖面精度一般很低，模糊程度达几十米左右。然而，微电位是分辨率最高的测井曲线，精度一般为0.125米一个数据点，但如果用该曲线代替波阻抗地震剖面，仍缺乏合适的方法。
与之类似的是，用低分辨率照片，过去认为是无法生成新的高清晰度照片的，但这没有考虑人的不断丰富的经验的影响，目前人脸复原技术，就是凭借丰富的经验和一些关键节点，将极其模糊的低分辨率照片，还原为清晰绘图的过程。对于经验丰富的人员，采用该方法识别的准确度目前为100％，可见低分辨率图像向高分辨率图像的过度，是可以凭借大量的经验来提高的，是可以实现的。
过去未能用微电位替换波阻抗，一是如果要建立微电极RN数据体的话，该数据体看似跟相邻两井的距离有关，即离井最近一根RN曲线，应该基本和该井的微电极相同，但远离井的数据点无法得出；二是如果将两井之间的某点的数据，仅仅对微电极进行简单的距离加权平均的话，等于浪费了全套地震数据；三是微电极与波阻抗有一定相关性，但目前找不到简便合适的计算式；四是新的微电极数据的位置也很难找准；五是不做归一化的微电位数据，在不同深度段，连数值范围都是不一样的。总之，由于找不到方法，用微电极替换地震剖面从未实现过。为此我们发明了一种新的采用拟合微电位曲线提高地震分辨率的方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种用目前分辨率最高的微电位曲线，对地震波阻抗数据进行改善和替换，大幅提高地震数据的精度的采用拟合微电位曲线提高地震分辨率的方法。
本发明的目的可通过如下技术措施来实现：采用拟合微电位曲线提高地震分辨率的方法，该采用拟合微电位曲线提高地震分辨率的方法包括：步骤1，对地震剖面进行精细地层解释，从而得到多个地层层面；步骤2，根据相邻两口井同一层位、同一位置的微电位数据点RN1、RN2，得到新层面的数据体RN’的位置；步骤3，根据相邻两井的波阻抗(I1、I2)，对同一层位、同一位置的微电位数据(RN1、RN2)进行变形，且按新数据点RN’距两井的距离(D1、D2)设定变形的权重；步骤4，采用步骤3中的计算方法，算出每个精细层面的全部微电极新数据RN’；步骤5，根据步骤2中新微电极数据RN’的位置，步骤3中微电极新数据RN’的大小，形成拟合微电位RN’新三维数据。
本发明的目的还可通过如下技术措施来实现：
在步骤2中，将相邻两口井同一层位、同一位置的微电位数据点RN1、RN2，沿着层面延伸，制定出整个新层面的数据体RN’的位置。
在步骤2中，新层面的新数据体的位置，为步骤3的变形后的数据RN’预留的位置。
在步骤3中，根据相邻两井的波阻抗(I1、I2)，将同一层位、同一位置的微电位数据(RN1、RN2)变换为RN1*I/I1、RN2*I/I2。
在步骤3中，按新数据点RN’距两井的距离(D1、D2)设定变形的权重，具体计算式为：
RN’＝RN1*I/I1*[D2/(D1+D2)]+RN2*I/I2*[D1/(D1+D2)]。
本发明中的采用拟合微电位曲线提高地震分辨率的方法，采用了一个简便易行的算式，对新数据进行求取，数据点的空间位置通过精细地震解释确定，解决了波阻抗替换成高精度微电位曲线的问题，从而大幅提高了地震数据的精度。该方法以同一井位地震波形与微电位作为起点，对紧邻一根波阻抗曲线进行不同深度的不规则缩放，如此用新建的虚拟微电位曲线，对低分辨率地震波阻抗曲线，替换，实现微电位剖面图的输出。在一定程度上借鉴了用低分辨率照片转变为人脸高精度图像的成功经验性，抓住几个关键点的信息(层面位置、距离影响、波阻抗的影响)，用目前分辨率最高的微电位曲线，对地震波阻抗数据进行改善和替换，大幅提高地震数据的精度。
附图说明
图1为本发明的采用拟合微电位曲线提高地震分辨率的方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的采用拟合微电位曲线提高地震分辨率的方法的流程图。
步骤1，对地震剖面进行精细地层解释，从而得到多个地层层面；
步骤2，将相邻两口井同一层位、同一位置的微电位数据点RN1、RN2，沿着层面延伸，制定出整个新层面的数据体RN’的位置；新层面的新数据体的位置，为步骤3的变形后的数据RN’预留的位置。
步骤3，用相邻两井的波阻抗(I1、I2)，对同一层位、同一位置的微电位数据(RN1、RN2)变换RN1*I/I1、RN2*I/I2，且按新数据点RN’距两井的距离(D1、D2)设定变形的权重，具体计算式为：
RN’＝RN1*I/I1*[D2/(D1+D2)]+RN2*I/I2*[D1/(D1+D2)]；
步骤4，采用该算式，算出每个精细层面的全部微电极新数据RN’；
步骤5，根据步骤2中新微电极数据RN’的位置，步骤3中微电极新数据RN’的大小，形成拟合微电位RN’新数据。
本发明中的采用拟合微电位曲线提高地震分辨率的方法，用微电位替换波阻抗数据，既考虑了距离加权因素，又用波阻抗对微电位进行了扭曲变换。该算法大大提高地震剖面的分辨率，将其由几十米的地震分辨率，提高到每0.125m一个RN'数据点，对地震资料升级、更精细的地层解释提供了新的方法。
以上实例的说明只是用于帮助理解本发明的核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (2)

1.采用拟合微电位曲线提高地震分辨率的方法，其特征在于，该采用拟合微电位曲线提高地震分辨率的方法包括：
步骤1，对地震剖面进行精细地层解释，从而得到多个地层层面；
步骤2，根据相邻两口井同一层位、同一位置的微电位数据点A、B，得到新层面的新数据体RN’的位置；
步骤3，根据相邻两井的波阻抗IA、IB，对同一层位、同一位置的微电位数据RNA、RNB进行变形，且按新层面的新数据体RN’距两井的距离DA、DB设定变形的权重；
步骤4，采用步骤3中的计算方法，算出每个精细层面的全部微电位新数据RNAB；
步骤5，根据步骤2中微电位新数据RNAB的位置、步骤3中微电位新数据RNAB的大小，形成拟合微电位RNAB新三维数据；
在步骤2中，新层面的新数据体的位置，为步骤3的变形后微电位新数据RNAB预留的位置；
在步骤3中，根据相邻两井的波阻抗IA、IB，将同一层位、同一位置的微电位数据RNA、RNB变形为RNA*IAB/IA、RNB*IAB/IB；按新数据点距两井的距离DA、DB设定变形的权重，具体求值公式为：
RNAB＝RNA*IAB/IA*[DB/(DA+DB)]+RNB*IAB/IB*[DA/(DA+DB)]。
2.根据权利要求1所述的采用拟合微电位曲线提高地震分辨率的方法，其特征在于，在步骤2中，将相邻两口井同一层位、同一位置的微电位数据点A、B，沿着层面延伸，制定出整个新层面的新数据体RN’的位置。
CN201811127695.4A 2018-09-26 2018-09-26 采用拟合微电位曲线提高地震分辨率的方法 Active CN109001803B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811127695.4A CN109001803B (zh) 2018-09-26 2018-09-26 采用拟合微电位曲线提高地震分辨率的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811127695.4A CN109001803B (zh) 2018-09-26 2018-09-26 采用拟合微电位曲线提高地震分辨率的方法

Publications (2)

Publication Number Publication Date
CN109001803A CN109001803A (zh) 2018-12-14
CN109001803B true CN109001803B (zh) 2021-02-23

Family

ID=64589412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811127695.4A Active CN109001803B (zh) 2018-09-26 2018-09-26 采用拟合微电位曲线提高地震分辨率的方法

Country Status (1)

Country Link
CN (1) CN109001803B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2444953B (en) * 2006-12-19 2009-07-22 Westerngeco Seismic Holdings Method for obtaining an image of a subsurface by regularizing irregularly sampled seismic data
CN101887132A (zh) * 2009-05-15 2010-11-17 中国石油天然气股份有限公司 一种井震结合定量预测砂岩储层流体饱和度的方法
CN103033846B (zh) * 2011-10-10 2016-06-08 中国石油化工股份有限公司 地质相控制的地震反演系统和地震反演方法
US9229127B2 (en) * 2013-02-21 2016-01-05 Saudi Arabian Oil Company Methods program code, computer readable media, and apparatus for predicting matrix permeability by optimization and variance correction of K-nearest neighbors
CN104991274B (zh) * 2015-07-03 2017-10-03 中国石油大学(华东) 多种地质因素定量约束下的单一圈闭级的有利区优选方法
CN106908842B (zh) * 2017-03-09 2019-06-11 中国石油天然气股份有限公司 一种确定单井沉积微相曲线的方法及装置
CN107065011B (zh) * 2017-06-22 2018-11-02 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法

Also Published As

Publication number Publication date
CN109001803A (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
CN107622104B (zh) 一种文字图像识别标注方法及系统
CN107688201B (zh) 基于rbm地震叠前信号聚类方法
CN102353985B (zh) 基于非下采样Contourlet变换的拟声波曲线构建方法
CN110276402B (zh) 一种基于深度学习语义边界增强的盐体识别方法
CN106022355B (zh) 基于3dcnn的高光谱图像空谱联合分类方法
CN107607992B (zh) 基于卷积神经网络的多波匹配方法
CN106897545B (zh) 一种基于深度置信网络的肿瘤预后预测系统
CN108830421B (zh) 致密砂岩储层的含气性预测方法及装置
CN111815773B (zh) 适用于机器学习算法的三维复杂地质模型标签制作方法
CN109061764A (zh) 一种分频融合波阻抗反演方法
CN103488971B (zh) 生物礁储层的几何形态识别方法
CN106886047A (zh) 一种接收函数和重力联合反演地壳厚度和波速比的方法
CN104375172A (zh) 火山岩下伏地层的构造形态正确成像的方法
CN112502701A (zh) 一种低渗透储层综合地质-工程的分类评价方法
CN115690908A (zh) 一种基于拓扑感知的三维手势姿态估计方法
CN109001803B (zh) 采用拟合微电位曲线提高地震分辨率的方法
CN109598068B (zh) 古构造约束建模方法、装置和设备
CN105334533B (zh) 一种微古地貌恢复方法
CN104714250B (zh) 实用的内幕小层自动解释方法
CN106504319B (zh) 井间储层三维对比图的生成方法及装置
CN107191184B (zh) 一种改进的埋藏史恢复方法
Bunton et al. Estimation of vocal tract area functions in children based on measurement of lip termination area and inverse acoustic mapping
CN106308758A (zh) 一种基于体温数据曲线的筛查方法及系统
CN112800518B (zh) 一种利用相邻地层互相关约束的地层面模型修正方法
CN107817521B (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
GR01 Patent grant
GR01 Patent grant