CN111337976B - 一种ovt域横波速度方位各向异性反演实现方法 - Google Patents

一种ovt域横波速度方位各向异性反演实现方法 Download PDF

Info

Publication number
CN111337976B
CN111337976B CN202010254060.1A CN202010254060A CN111337976B CN 111337976 B CN111337976 B CN 111337976B CN 202010254060 A CN202010254060 A CN 202010254060A CN 111337976 B CN111337976 B CN 111337976B
Authority
CN
China
Prior art keywords
azimuth
sectors
inversion
transverse wave
gather
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
CN202010254060.1A
Other languages
English (en)
Other versions
CN111337976A (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.)
Beijing Elsiwave Technology Co ltd
Shanghai Borui Energy Technology Co ltd
Original Assignee
Shanghai Borui Energy Technology Co ltd
Beijing Elsiwave Technology Co ltd
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 Shanghai Borui Energy Technology Co ltd, Beijing Elsiwave Technology Co ltd filed Critical Shanghai Borui Energy Technology Co ltd
Priority to CN202010254060.1A priority Critical patent/CN111337976B/zh
Publication of CN111337976A publication Critical patent/CN111337976A/zh
Application granted granted Critical
Publication of CN111337976B publication Critical patent/CN111337976B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/303—Analysis for determining velocity profiles or travel times
    • 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/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack
    • 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
    • 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/626—Physical property of subsurface with anisotropy
    • 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/64—Geostructures, e.g. in 3D data cubes
    • G01V2210/642—Faults
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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

本发明公开了在勘探开发地球物理技术领域一种OVT域横波速度方位各向异性反演实现方法。该方法通过不同方位扇区划分两个原则和方法、道集优化原则即主要去除高角度噪声、每一个方位扇区开展没有井约束的地震反演等关键环节建立了一套完整的OVT域道集横波速度方位各向异性反演技术方法和流程。本发明解决了目前使用宽方位OVT域道集数据进行裂缝直接检测存在诸多难题，如不同方位扇区的划分原则和方法、道集优化处理原则、以及横波速度地震反演中否采用测井数据进行约束等。使用本发明的技术方法将显著提高利用宽方位地震数据直接检测裂缝的可靠性和精度。本发明的方法虽然是针对OVT域道集数据而设计，但很容易推广到分方位叠前偏移的道集数据。

Description

一种OVT域横波速度方位各向异性反演实现方法
技术领域
本发明属于勘探开发地球物理技术地震预测领域。
背景技术
横波分裂又叫横波双折射，1970年代最早被Crampin等人发现并证实。Crampin的研究小组利用人工地震地壳测深和天然地震资料，发现并证实了横波分裂现象的存在，而且快横波的偏振方向与应力场方向一致。
横波通过各向异性介质将分裂为速度各异的两个波。一个波偏振方向平行于裂缝走向，速度较快，称为快横波；另一个偏振方向垂直于裂缝走向，速度较慢，称为慢横波。
利用横波分裂的思想，首先对宽方位地震道集进行不同方位扇区划分，在每一个方位扇区内，我们将其近似看作各向同性，这样采用各向同性Zoeppritz线性方求解横波速度。最后通过椭圆拟合处理的方法，找到快横波的方向和大小，以及慢横波的大小，并进一步求出快横波速度与慢横波速度比值，即各向异性的强度。
裂缝地震直接检测主要是基于研究宽方位地震数据的方位各向异性，而横波速度的方位各向异性最强。通过椭圆拟合处理技术可以确定方位各向异性主方向(快方向即椭圆长轴方向)和各向异性强度(椭圆长轴和短轴之比)。以往方位各向异性反演主要基于叠后地震属性，例如振幅、衰减、P波阻抗等等，即便是使用叠前地震数据也几乎没有使用叠前反演的横波速度数据。
我们基于成熟的OVT(Offset Vector Tile炮检距向量片)处理技术和横波速度叠前反演技术，发明了一套基于横波分裂思想的横波速度方位各向异性反演方法和流程，其中的几个关键环节包括方位扇区划分原则和方法、方位扇区道集优化处理原则、方位扇区横波速度反射率无井约束地震反演等。
发明内容
本发明旨在解决基于OVT域道集方位各向异性反演完整的技术方法和流程。本发明采用以下技术解决方案。首先开展不同方位扇区划分，然后进行不同方位扇区道集优化处理、不同方位扇区横波速度弹性地震体反演、不同方位扇区无测井数据约束的横波速度反演，最后进行椭圆拟合处理并确定方位各向异性主方向(即裂缝方向)和各向异性强度(椭圆长轴/短轴比即裂缝密度)。
本发明步骤如下：
步骤一、根据OVT域道集数据中每一道偏移距数据的方位分布，开展方位扇区划分。方位扇区划分原则有两个：(1)根据单个扇区覆盖次数不低于35次的标准将360度方位OVT域道集划分为5～6个扇区，扇区太少不能保证椭圆拟合精度，扇区太多不能保证足够的覆盖次数从而导致叠前反演信噪比较低。(2)确保每个扇区都包含近、中、远偏移距数据，且各扇区最大偏移距基本相同，避免远偏移距不同造成不同扇区各向异性差异。
步骤二、各方位扇区叠前偏移后道集优化处理，主要去除高角度相干噪声(线性或双曲线多次波噪声)，尽量不做小时差噪声去除，因为这很可能抹杀道集上本来存在的方位各向异性。
步骤三、针对每一个不同方位扇区道集数据，利用公式(1)开展三参数AVO反演得到横波速度弹性地震体(反射率)，这里之所以使用横波速度是因为横波速度方位各向异性最强。
式中：Rp为反射系数；θ为入射角；vp为纵波速度；vs为横波速度；p为密度。
步骤四、针对每一个不同方位扇区横波速度弹性地震体，逐一开展横波速度无测井数据约束的地震反演。因为如果地震反演中使用测井数据约束，则会导致不同方位扇区横波速度趋于一致，这样就抹平了各方位扇区之间应有的横波速度的差异。方位各向异性反演就是通过不同方位扇区叠前反演找出弹性参数在不同方位扇区的差异。
步骤五、利用不同方位扇区(5个或6个)横波速度数据开展椭圆拟合处理，得到椭圆长轴(快横波)方位和大小、短轴(慢横波)大小并进一步得到长轴/短轴比(方位各向异性强度)。长轴方位就是各向异性的方向也即裂缝方向，长轴/短轴比就是各向异性强度也即裂缝密度。
本发明产生的效果包括提高方位各向异性反演的精度，使得裂缝直接检测更可靠，并将提高勘探成功率以及优化水平钻井和压裂设计。
附图说明
宽方位地震数据方位扇区划分示意图(图中每一根直线条代表该道数据偏移距大小和方位，线段长短代表偏移距大小，线段位置代表方位)。
(本示意图共划分了6个不同方位扇区，即扇区I：0°-40°：扇区II：44°-130°；扇区III：130°-180°；扇区IV：180°-224°；扇区V：224°-310°；扇区VI：310°-360°)
具体实施方式
以下结合附图对本发明的具体实施步骤进行说明。
本发明是这样实现的。
步骤一、根据OVT域道集数据中每一道偏移距数据的方位分布，开展方位扇区划分。方位扇区划分原则有两个：(1)根据单个扇区覆盖次数不低于35次的标准将360度方位OVT域道集划分为5～6个扇区，扇区太少不能保证椭圆拟合精度，扇区太多不能保证足够的覆盖次数从而导致叠前反演信噪比较低。(2)确保每个扇区都包含近、中、远偏移距数据，且各扇区最大偏移距基本相同，避免远偏移距不同造成不同扇区方位各向异性差异。如果地震数据总覆盖次数达到210次以上，可以划分为6个方位扇区，每个扇区覆盖次数不少于35次；如果地震数据覆盖次数低于210次，可以划分为5个方位扇区。
步骤二，各扇区叠前偏移后道集优化处理，可以采用诸如高分辨Radon变换技术主要去除高角度相干噪声(线性或双曲线多次波噪声)，尽量不做小时差噪声去除，因为这很可能抹杀道集上本来存在的方位各向异性。
步骤三，针对每一个不同方位扇区道集数据，利用公式(1)开展三参数AVO反演得到横波速度弹性地震体(反射率)。
步骤四，针对每一个不同方位扇区横波速度弹性地震体，逐一开展横波速度无测井数据约束的地震反演。因为如果地震反演中使用测井数据约束，则会导致不同方位扇区横波速度趋于一致，这样就抹平了各方位扇区之间应有的横波速度的差异。方位各向异性反演目标就是通过不同方位扇区叠前反演找出弹性参数在不同方位扇区的差异。
步骤五，利用不同方位扇区(5个或6个)横波速度数据开展椭圆拟合处理，得到椭圆长轴(快横波)方位和大小、短轴(慢横波)大小并进一步得到长轴/短轴比(方位各向异性强度)。长轴方位就是各向异性的方向也即裂缝方向，长轴/短轴比就是各向异性强度也即裂缝密度。
本发明的方法虽然是针对OVT域道集数据而设计，但很容易推广到分方位叠前偏移的道集数据，即所有的原则和方法都适用于分方位叠前偏移的道集数据。二者的区别仅是方位扇区划分的顺序，OVT域道集扇区划分是在叠前偏移之后，而分方位叠前偏移的道集数据是在叠前偏移前先划分好方位扇区。

Claims (1)

1.一种OVT域横波速度方位各向异性反演实现方法，其特征在于，方位扇区划分方法、道集优化只去除高角度噪声、每一个方位扇区开展没有井约束的地震反演；其实现步骤如下：
步骤一、根据OVT域道集数据中每一道偏移距数据的方位分布，开展方位扇区划分；
方位扇区划分原则和方法：(1)根据单个扇区覆盖次数不低于35次的标准将360度方位OVT域道集划分为5~6个扇区，扇区太少不能保证椭圆拟合精度，扇区太多不能保证足够的覆盖次数从而导致叠前反演信噪比较低；(2)确保每个扇区都包含近、中、远偏移距数据，且最大偏移距基本相同，避免远偏移距不同造成不同扇区各向异性差异；
步骤二、每个方位扇区道集优化处理，主要去除高角度相干噪声，尽量不做小角度或时差噪声去除，因为这很可能抹杀道集上本来存在的各向异性差异，所述高角度相干噪声包括线性或双曲线多次波噪声；
步骤三、针对每一个不同方位扇区道集数据，开展叠前横波速度反演得到横波速度弹性地震体；这里之所以使用横波速度是因为横波速度方位各向异性最强；
步骤四、针对每一个不同方位扇区横波速度弹性地震体，逐一开展横波速度没有井约束地震反演得到不同方位扇区横波速度绝对值数据；因为如果地震反演中使用井数据约束，则会导致不同方位扇区横波速度趋于一致，这样就失去了开展各向异性反演的意义；各向异性反演就是通过不同方位扇区叠前反演找出弹性参数在不同方位的差异；
步骤五、利用不同方位扇区横波速度数据开展椭圆拟合处理，得到椭圆长轴方位和大小、短轴大小并进一步得到长轴/短轴比；长轴方位就是方位各向异性的方向也即裂缝方向，长轴/短轴比就是各向异性强度也即裂缝密度。
CN202010254060.1A 2020-04-01 2020-04-01 一种ovt域横波速度方位各向异性反演实现方法 Active CN111337976B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010254060.1A CN111337976B (zh) 2020-04-01 2020-04-01 一种ovt域横波速度方位各向异性反演实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010254060.1A CN111337976B (zh) 2020-04-01 2020-04-01 一种ovt域横波速度方位各向异性反演实现方法

Publications (2)

Publication Number Publication Date
CN111337976A CN111337976A (zh) 2020-06-26
CN111337976B true CN111337976B (zh) 2023-08-25

Family

ID=71184659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010254060.1A Active CN111337976B (zh) 2020-04-01 2020-04-01 一种ovt域横波速度方位各向异性反演实现方法

Country Status (1)

Country Link
CN (1) CN111337976B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
CN107045143A (zh) * 2017-02-22 2017-08-15 中国石油天然气股份有限公司 一种预测裂缝发育的方法及装置
CN107728204A (zh) * 2016-08-11 2018-02-23 中国石油化工股份有限公司 基于叠前纵波各向异性的裂缝预测方法及系统
CN109655905A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 基于全方位角数据叠前avaz反演裂缝预测方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
CN107728204A (zh) * 2016-08-11 2018-02-23 中国石油化工股份有限公司 基于叠前纵波各向异性的裂缝预测方法及系统
CN107045143A (zh) * 2017-02-22 2017-08-15 中国石油天然气股份有限公司 一种预测裂缝发育的方法及装置
CN109655905A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 基于全方位角数据叠前avaz反演裂缝预测方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑多明等.基于OVT域五维道集碳酸盐岩叠前裂缝预测技术.《物探化探计算技术》.2020,第42卷(第1期),第10-12页. *

Also Published As

Publication number Publication date
CN111337976A (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
Pérez et al. Fracture detection in a carbonate reservoir using a variety of seismic methods
WO2016041189A1 (zh) 一种评价页岩气储层及寻找甜点区的方法
US20030167126A1 (en) Layer stripping converted reflected waveforms for dipping fractures
Wild Practical applications of seismic anisotropy
Angerer et al. Characterization of dipping fractures using PS mode-converted data
Diaz-Acosta et al. Investigation of fractured carbonate reservoirs by applying shear-wave splitting concept.
Lynn et al. Reflection shear-wave data along the principal axes of azimuthal anisotropy
CN111337976B (zh) 一种ovt域横波速度方位各向异性反演实现方法
Liao et al. An integrated study on fracture prediction using 3D P-wave field seismic data
CN110780345A (zh) 一种隧道超前地震探测地震数据的三维速度分析方法
ZHAO Application of multi-component seismic exploration in the exploration and production of lithologic gas reservoirs
Zeng* et al. Wide-azimuth seismic data processing technology for small fault and crack identification
Pei et al. Research and application of 5D seismic prediction technology
Wang et al. Application of PrestackFracture Prediction Technique in Deep Gas Exploration in Northern Songliao Basin
Deng et al. Application of combined imaging technology of multi azimuth seismic data in deepwater rugged seabed area
Zeng et al. Hti Medium Azimuthal Anisotropy Correction Technology and Application of Wide Azimuth Seismic Data
Zhu et al. Fracture prediction of buried hill based on wide azimuth seismic data in Bohai Bay
Larner et al. Investigation of azimuthal anisotropy from offset VSP data-a case study
Zhao et al. Application of seismic multi-attribute technique in fracture prediction of carbonate reservoir
Zeng et al. Research and application of integrated technology of tight sandstone gas acquisition processing and interpretation
Wu et al. Research and Application of Crack Prediction Technology for Pre-Stack Seismic Data
Wang et al. Fracture Prediction of Carbonate Karst Reservoirs Using Full-Azimuth Multi-wave Seismic Data
CN110687604A (zh) 一种无宽方位地震资料条件下纵、横波差异处理检测裂缝方法
Yi et al. Prestack depth migration processing and interpretation of 2D seismic data in Kedong tectonic zone, southwest of Tarim Basin
Xuming et al. The influence of the near-surface anisotropy

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
TA01 Transfer of patent application right

Effective date of registration: 20230725

Address after: No. 400, Floor 4, Building 3, No. A1, Dongwangzhuang, Wudaokou, Haidian District, Beijing 100080

Applicant after: BEIJING ELSIWAVE TECHNOLOGY Co.,Ltd.

Applicant after: Shanghai Borui Energy Technology Co.,Ltd.

Address before: No. 400, Floor 4, Building 3, No. A1, Dongwangzhuang, Wudaokou, Haidian District, Beijing 100083

Applicant before: BEIJING ELSIWAVE TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant