CN111337979B - 一种用解析方法确定地质界面真倾角的方法 - Google Patents
一种用解析方法确定地质界面真倾角的方法 Download PDFInfo
- Publication number
- CN111337979B CN111337979B CN202010194564.9A CN202010194564A CN111337979B CN 111337979 B CN111337979 B CN 111337979B CN 202010194564 A CN202010194564 A CN 202010194564A CN 111337979 B CN111337979 B CN 111337979B
- Authority
- CN
- China
- Prior art keywords
- true
- angle
- gamma
- measuring
- psi
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004458 analytical method Methods 0.000 title claims abstract description 14
- 238000010276 construction Methods 0.000 abstract description 4
- 238000013507 mapping Methods 0.000 abstract description 2
- 239000013598 vector Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/22—Measuring 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
Abstract
一种用解析方法确定地质界面真倾角的方法,属于地震勘探领域。地震勘探要获得勘探面的真深度,必须知道真倾角。在勘探区域内,真倾角可由视倾角求得。目前用作图法确定地质界面真倾角比较直观,但误差较大,同时在野外施工作图多有不便。据此,提出用解析方法确定倾斜界面真倾角。采用本发明提出的用解析方法,需要计算两测线的视倾角,测定两测线在水平面内夹角,经计算机的迭代计算,就可以得到地质界面的真倾角。解析方法误差小、便于操作、比作图法更适合于野外施工。
Description
技术领域
本发明属于地震勘探领域,具体涉及一种用解析方法确定地质界面真倾角的方法。
背景技术
目前地震勘探中要获得地质界面的真深度,必须在确定法线深度基础上,知道真倾角或视倾角和测线方位角。在勘探区域内,真倾角可由两条相交的测线视倾角求得。参考文献:陆基孟,王永刚.地震勘探原理[M].第三版,山东东营:中国石油大学出版社,2011,199-200页。
地震勘探中视倾角的确定方法如下:在上覆均匀介质情况下,假设叠加剖面上有一条倾斜的同相轴,如图1所示,在,两点的自激自收时间分别为和,对应的深度剖面如图2所示。设定界面的视倾角为ψ,界面上覆盖介质的波速为v,其它符号见图2,h1和h2分别为两个接收点法线深度,Δh为两者法线深度之差,Δx为两接收点之间距离,Δt0为两点自激自收时间差,由此可以得到:
得到视倾角后,目前用作图法获得真倾角,具体做法如下:已知两条可以相交的测线1和测线 2,用上述方法确定两者视倾角分别为ψ1和ψ2。用如图3中作图法确定真倾角:首先在平面上画出两条测线,点O为两条测线交点,从O点分别沿着两条测线方向画矢量和,选择适当比例尺,使得两个矢量的大小分别为sinψ1和sinψ2 ,然后分别沿两个矢量的端点A1和A2作测线的垂线,相交于B,连接OB,矢量的长度就等于界面的真倾角正弦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;
θ=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:在上述数据文件中,拾取Δθ的最小值,其对应的γ值就是真倾角。
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 CN111337979A (zh) | 2020-06-26 |
CN111337979B true 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)
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 |
-
2020
- 2020-03-19 CN CN202010194564.9A patent/CN111337979B/zh not_active Expired - Fee Related
Patent Citations (7)
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)
Title |
---|
A new method of restraining the Fabry-Perot resonance;Yin, JZ等;《Optical and Quantum Electronics》;20030701;第35卷(第9期);第873-878页 * |
从两相交测线的视倾角求取;应秀芝,康寅巧;《石油地球物理勘探》;19941230;第29卷;第189-190页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111337979A (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4357660A (en) | Formation dip and azimuth processing technique | |
US10520632B2 (en) | Rotating magnetic field range finder and its measurement method for relative separation between drilling wells | |
CN105783919A (zh) | H型标量传感器阵列对磁性目标的追踪定位方法 | |
EP2690468B1 (en) | A streamer for seismic prospection comprising tilt compensation of directional sensors | |
CN112558006A (zh) | 一种用于海底缆线的双三轴电磁探测及路由定位方法 | |
CN105674945A (zh) | 一种基于mems传感器的海底滑坡监测装置及监测方法 | |
CN108507534A (zh) | Rtk倾斜测量方法和系统 | |
WO2017071440A1 (zh) | 一种地震全向矢量检波器 | |
CN111337979B (zh) | 一种用解析方法确定地质界面真倾角的方法 | |
WO2017071436A1 (zh) | 一种检波器的质控方法、装置、计算机可读存储介质及设备 | |
CN102419457B (zh) | 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法 | |
CN102564653B (zh) | 高地应力地区存在岩芯饼化现象时测量地应力的方法 | |
CN102721966A (zh) | 相干类测深声纳正下方高精度测深方法及系统 | |
CN103591936B (zh) | 一种倾斜角的测量方法 | |
CN104965232A (zh) | 低纬度地区磁构造格架自动提取方法 | |
CN100559400C (zh) | 位场多方向多尺度边缘检测方法 | |
WO2017071435A1 (zh) | 一种地震全向矢量旋度检波器 | |
US8340912B2 (en) | Seismic attributes for structural analysis | |
US11041720B2 (en) | Calibration of a magnetometer in a towed object telemetry unit based on turn data | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
CN109752762A (zh) | 单发多收观测装置瞬变电场数据动校正方法及装置 | |
CN114895367A (zh) | 岩体产状信息测量方法 | |
CN109507730B (zh) | 一种改进的确定折射面法线深度的方法 | |
CN110398777B (zh) | 一种井中检波器定位方法 | |
CN111505715A (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 | ||
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 |