CN111337979A - 一种用解析方法确定地质界面真倾角的方法 - Google Patents

一种用解析方法确定地质界面真倾角的方法 Download PDF

Info

Publication number
CN111337979A
CN111337979A CN202010194564.9A CN202010194564A CN111337979A CN 111337979 A CN111337979 A CN 111337979A CN 202010194564 A CN202010194564 A CN 202010194564A CN 111337979 A CN111337979 A CN 111337979A
Authority
CN
China
Prior art keywords
true
angle
measuring
gamma
geological interface
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.)
Granted
Application number
CN202010194564.9A
Other languages
English (en)
Other versions
CN111337979B (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.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN202010194564.9A priority Critical patent/CN111337979B/zh
Publication of CN111337979A publication Critical patent/CN111337979A/zh
Application granted granted Critical
Publication of CN111337979B publication Critical patent/CN111337979B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/22Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring angles or tapers; for testing the alignment of axes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种用解析方法确定地质界面真倾角的方法,属于地震勘探领域。地震勘探要获得勘探面的真深度,必须知道真倾角。在勘探区域内,真倾角可由视倾角求得。目前用作图法确定地质界面真倾角比较直观,但误差较大,同时在野外施工作图多有不便。据此,提出用解析方法确定倾斜界面真倾角。采用本发明提出的用解析方法,需要计算两测线的视倾角,测定两测线在水平面内夹角,经计算机的迭代计算,就可以得到地质界面的真倾角。解析方法误差小、便于操作、比作图法更适合于野外施工。

Description

一种用解析方法确定地质界面真倾角的方法
技术领域
本发明属于地震勘探领域,具体涉及一种用解析方法确定地质界面真倾角的方法。
背景技术
目前地震勘探中要获得地质界面的真深度,必须在确定法线深度基础上,知道真倾角或视倾角和测线方位角。在勘探区域内,真倾角可由两条相交的测线视倾角求得。参考文献:陆基孟,王永刚.地震勘探原理[M].第三版,山东东营:中国石油大学出版社,2011,199-200页。
地震勘探中视倾角的确定方法如下:在上覆均匀介质情况下,假设叠加剖面上有一条倾斜的同相轴,如图1所示,在
Figure 61596DEST_PATH_IMAGE002
Figure 385261DEST_PATH_IMAGE004
两点的自激自收时间分别为
Figure 696157DEST_PATH_IMAGE006
Figure 747158DEST_PATH_IMAGE008
,对应的深度剖面如图2所示。设定界面的视倾角为ψ,界面上覆盖介质的波速为v,其它符号见图2,h1和h2分别为两个接收点法线深度,Δh为两者法线深度之差,Δx为两接收点之间距离,Δt0为两点自激自收时间差,由此可以得到:
sinψ
Figure 92689DEST_PATH_IMAGE010
(1)
得到视倾角后,目前用作图法获得真倾角,具体做法如下:已知两条可以相交的测线1和侧线2,用上述方法确定两者视倾角分别为ψ1和ψ2。用如图3中作图法确定真倾角:首先在平面上画出两条测线,点O为两条测线交点,从O点分别沿着两条测线方向画矢量
Figure 321676DEST_PATH_IMAGE012
Figure 588709DEST_PATH_IMAGE014
,选择适当比例尺,使得两个矢量的大小分别为sinψ1和sinψ2 ,然后分别沿两个矢量的端点A1和A2作测线的垂线,相交于B,连接OB,矢量
Figure 318768DEST_PATH_IMAGE016
的长度就等于界面的真倾角正弦sinγ,再进一步求出真倾角γ。
目前用作图法确定地质界面真倾角误差较大,同时在野外施工作图多有不便。为提高地震勘探野外工作中确定真倾角方法的实用性,提出一种用解析方法确定地质界面真倾角的方法。
发明内容
如图4,平面MDEN为水平面,平面FDEG是地面以下的倾斜地质界面,其走向是DE方向。DA和EA是地面上的两条测线,测线平移不改变视倾角,所以设两条测线相交于A点,两条测线的视倾角分别为ψ1和ψ2,AC垂直于倾斜地质界面FDEG,交于C点,在水平面MDEN中,作AB垂直于DE,由于AC垂直于倾斜地质界面FDEG,由于DE垂直于AB,还垂直于AC,所以DE垂直于平面ABC,可以判断出DE垂直于CB,所以有:
在直角三角形ACD中, sinψ1=AC/AD,所以AD=AC/sinψ1
在直角三角形ACE中, sinψ2=AC/AE,所以AE=AC/sinψ2
在直角三角形ACB中, sinγ=AC/AB, 所以AB=AC/sinγ;
在三角形DAE中,测定得到∠DAE=θ,∠θ=∠DAB+∠EAB;
Figure 643439DEST_PATH_IMAGE018
也就是:
θ=arccos(sinψ1/sinγ)+arccos(sinψ2/ sinγ)
上式中隐含γ,无法求出显式解,可以通过迭代方法计算出。为了进行迭代计算,令:
Δθ=|θ-arccos(sinψ1/sinγ)-arccos(sin ψ2/ sinγ) | (2)
取ψ1和ψ2中较大者为γ下限,取π/2为γ上限,按一定步长Δγ增加γ,用公式(2)计算,Δγ尽量取得小一些,例如取Δγ=π/7200,以便得到更精确的真倾角γ。建立数据文件,记录 Δθ和γ,获取Δθ最小值,并取Δθ最小值对应的γ,即得到了真倾角γ。
一种用解析方法确定地质界面真倾角的方法,具体流程如下:
步骤1:获取以弧度为单位的两测线1和测线2的视倾角ψ1和ψ2,画出表示测线1和测线2的直线,使之相交,交点为O;
步骤2:用量角器测量图3中的∠A2OA1,也就是图4中的θ,并转换成以弧度为单位;
步骤3:取ψ1和ψ2中较大者为γ下限,取π/2为γ上限,以Δγ=π/7200为步长,用公式(2)计算,建立数据文件,记录Δθ和γ;
Δθ=|θ-arccos(sinψ1/sinγ)-arccos(sin ψ2/ sinγ) | (2)
步骤4:在上述数据文件中,拾取Δθ的最小值,其对应的γ值就是真倾角。
有益技术效果:
目前地震勘探计算倾斜面真倾角采用绘图法,误差大,且不适合在野外实施;采用本发明提出的解析方法,仅仅需要测定两测线在水平面内夹角、计算得到两测线的视倾角,经计算机的迭代计算,就可以得到倾斜界面的真倾角。解析方法误差小、便于操作、更适合于野外施工。
附图说明
图1 为倾斜反射面上的同相轴示意图;
图2 为目前用绘图法确定视倾角示意图;
图3 为本发明实施例的确定真倾角示意图;
图4 为本发明实施例的用解析法确定真倾角示意图;
图5为本发明实施例的一种用解析方法确定地质界面真倾角的方法流程图。
具体实施方式
下面结合附图和具体实施实例对发明做进一步说明,一种用解析方法确定地质界面真倾角的方法,如图5所示,具体流程如下:
步骤1:获取以弧度为单位的两测线1和测线2的视倾角ψ1和ψ2,画出表示测线1和测线2的直线,使之相交,交点为O,在实施例中,ψ1=π/4,ψ2=π/6,如图3;
步骤2:用量角器测量图3中的∠A2OA1,也就是图4中的θ,并转换成以弧度为单位,在实施例中θ=π/3;
步骤3:取ψ1和ψ2中较大者为γ下限,即π/4为下限,取π/2为γ上限,以Δγ=π/7200为步长,用公式(2)计算,建立数据文件,记录Δθ和γ;
Δθ=|θ-arccos(sinψ1/sinγ)-arccos(sin ψ2/ sinγ) | (2)
步骤4:在上述数据文件中,拾取Δθ的最小值,其对应的γ值就是真倾角,本实施例中计算得到γ=46.65°。

Claims (1)

1.一种用解析方法确定地质界面真倾角的方法,属于地震勘探领域,其特征在于,具体流程如下:
步骤1:获取以弧度为单位的两测线1和测线2的视倾角ψ1和ψ2,画出表示测线1和测线2的直线,使之相交;
步骤2:用量角器测量两测线间夹角θ,并转换成以弧度为单位;
步骤3: 取ψ1和ψ2中较大者为真倾角γ下限,取π/2为γ上限,以Δγ=π/7200为步长,用公式计算,建立数据文件,记录Δθ和γ;
所述公式为:
Δθ=|θ-arccos(sinψ1/sinγ)-arccos(sin ψ2/ sinγ)|
步骤4:在上述数据文件中,拾取Δθ的最小值,其对应的γ值就是真倾角。
CN202010194564.9A 2020-03-19 2020-03-19 一种用解析方法确定地质界面真倾角的方法 Expired - Fee Related CN111337979B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010194564.9A CN111337979B (zh) 2020-03-19 2020-03-19 一种用解析方法确定地质界面真倾角的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010194564.9A CN111337979B (zh) 2020-03-19 2020-03-19 一种用解析方法确定地质界面真倾角的方法

Publications (2)

Publication Number Publication Date
CN111337979A true CN111337979A (zh) 2020-06-26
CN111337979B CN111337979B (zh) 2022-08-16

Family

ID=71186197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010194564.9A Expired - Fee Related CN111337979B (zh) 2020-03-19 2020-03-19 一种用解析方法确定地质界面真倾角的方法

Country Status (1)

Country Link
CN (1) CN111337979B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4092779A (en) * 1976-01-30 1978-06-06 Bozidar M Stokic Device for determining angle relationships of or between planes
CN1517720A (zh) * 2003-01-13 2004-08-04 大庆石油管理局 二维地震叠偏剖面成图空校方法
CN101950030A (zh) * 2010-08-30 2011-01-19 钟世航 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
CN105093307A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 下古生界斜地层真厚度求取方法
CN105783849A (zh) * 2016-05-09 2016-07-20 辽宁工程技术大学 一种岩层真倾角和真走向的测量系统及方法
CN106560826A (zh) * 2016-05-31 2017-04-12 孙庆先 利用煤层赋存要素和断层要素确定断煤交线位置的方法
WO2017160273A1 (en) * 2016-03-14 2017-09-21 Schlumberger Technology Corporation Dip estimation via a modified structure tensor

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4092779A (en) * 1976-01-30 1978-06-06 Bozidar M Stokic Device for determining angle relationships of or between planes
CN1517720A (zh) * 2003-01-13 2004-08-04 大庆石油管理局 二维地震叠偏剖面成图空校方法
CN101950030A (zh) * 2010-08-30 2011-01-19 钟世航 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
CN105093307A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 下古生界斜地层真厚度求取方法
WO2017160273A1 (en) * 2016-03-14 2017-09-21 Schlumberger Technology Corporation Dip estimation via a modified structure tensor
CN105783849A (zh) * 2016-05-09 2016-07-20 辽宁工程技术大学 一种岩层真倾角和真走向的测量系统及方法
CN106560826A (zh) * 2016-05-31 2017-04-12 孙庆先 利用煤层赋存要素和断层要素确定断煤交线位置的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YIN, JZ等: "A new method of restraining the Fabry-Perot resonance", 《OPTICAL AND QUANTUM ELECTRONICS》 *
应秀芝,康寅巧: "从两相交测线的视倾角求取", 《石油地球物理勘探》 *

Also Published As

Publication number Publication date
CN111337979B (zh) 2022-08-16

Similar Documents

Publication Publication Date Title
US20140081614A1 (en) Geological Exploration Method for Making Plan and Elevation Drawings Directly by Rotational Tin Network and Non-Profiling Method
Qi et al. Volumetric aberrancy to map subtle faults and flexures
CN105783919A (zh) H型标量传感器阵列对磁性目标的追踪定位方法
CN105674945A (zh) 一种基于mems传感器的海底滑坡监测装置及监测方法
CN210719199U (zh) 一种水下机器人多设备组合导航系统
CN110806209A (zh) 一种水下机器人多设备组合导航系统及方法
CN112540371A (zh) 一种近底多波束坐标转换处理方法
CN108507534A (zh) Rtk倾斜测量方法和系统
WO2017071436A1 (zh) 一种检波器的质控方法、装置、计算机可读存储介质及设备
CN111337979B (zh) 一种用解析方法确定地质界面真倾角的方法
CN102419457B (zh) 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法
CN102721966A (zh) 相干类测深声纳正下方高精度测深方法及系统
CN101900831A (zh) 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN103591936B (zh) 一种倾斜角的测量方法
US20210270607A1 (en) Calibration of a Magnetometer in a Towed Object Telemetry Unit Based on Turn Data
WO2017071435A1 (zh) 一种地震全向矢量旋度检波器
CN112684498A (zh) 一种基于宽方位地震数据的储层裂缝预测方法及系统
CN101290223A (zh) 海底地形测量波束角效应改正方法
CN109507730B (zh) 一种改进的确定折射面法线深度的方法
Ostashev et al. Ray tracing in a stratified moving atmosphere: Azimuthal deviation in sound propagation
CN111505715A (zh) 一种基于深度域中心入射角计算裂缝参数的方法
CN110275145A (zh) 探地雷达测量误差计算方法及装置
Yu et al. Arctangent function‐based third derivative attribute for characterisation of faults
CN107990896B (zh) 一种基于xyz格式地图的路径生成方法
Benites et al. A ground-motion transfer function matrix between two nearby rock and soil sites: a system identification problem

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

Granted publication date: 20220816

CF01 Termination of patent right due to non-payment of annual fee