CN109581503B - 一种解决沙漠探区层析静校正边界问题的方法及系统 - Google Patents

一种解决沙漠探区层析静校正边界问题的方法及系统 Download PDF

Info

Publication number
CN109581503B
CN109581503B CN201710899715.9A CN201710899715A CN109581503B CN 109581503 B CN109581503 B CN 109581503B CN 201710899715 A CN201710899715 A CN 201710899715A CN 109581503 B CN109581503 B CN 109581503B
Authority
CN
China
Prior art keywords
sand dune
point
sdb
static correction
dune
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
CN201710899715.9A
Other languages
English (en)
Other versions
CN109581503A (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
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710899715.9A priority Critical patent/CN109581503B/zh
Publication of CN109581503A publication Critical patent/CN109581503A/zh
Application granted granted Critical
Publication of CN109581503B publication Critical patent/CN109581503B/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
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • 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/53—Statics correction, e.g. weathering layer or transformation to a datum

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)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提出了一种解决沙漠探区层析静校正边界问题的方法及系统，该方法包括：将接收点分成属于有炮点区和属于无炮点区两部分；计算接收点的沙丘低速带校正量和沙丘厚度；计算接收点的沙丘底部面下低速带静校正量和沙丘时间；利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；利用拟沙丘曲线函数计算属于无炮点区接收点的静校正量。本发明的方法应用于新疆沙漠地区多块三维地震勘探工区的资料处理，都很好地解决了层析静校正存在的边界问题，效果非常明显。

Description

一种解决沙漠探区层析静校正边界问题的方法及系统
技术领域
本发明属于地震勘探资料数据处理领域，尤其涉及沙漠工区地震勘探静校正计算技术。
背景技术
基于初至时间的层析静校正技术是解决沙漠工区静校正问题的有效手段。初至层析近地表速度模型反演的精度取决于近道初至信息，在有近道初至信息的有炮点区反演的近地表速度模型可靠，相应地计算的静校正量也可靠。而在工区四周无炮点区，由于缺乏近道初至信息，反演的近地表速度模型精度明显下降，计算的接收点静校正量精度随之下降。所有基于初至层析的静校正技术都有这样的边界问题，只是程度上有差异。
现有的解决方法多采用把有炮点区近地表速度模型或静校正量进行外推，外推的精度决定了静校正精度，工区四周无炮点区接收点静校正量的精度难以令人满意。
发明内容
本发明提出一种“拟沙丘曲线法”，用于解决接收点静校正边界问题。应用有炮点区初至层析获得的接收点静校正量，统计出一条“拟沙丘曲线”函数，用以计算工区四周无炮点区的接收点静校正量，快速有效地解决沙漠勘探工区初至层析接收点静校正量的边界问题。
根据本发明的一个方面，提供一种解决沙漠探区层析静校正边界问题的方法，该方法包括：
将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点的沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；
最后输出接收点静校正量。
进一步地，输入数据与参数，包括：工区炮点平面坐标、接收点平面坐标、地表高程、沙丘底部面高程、接收点层析静校正量、静校正替换速度和基准面高程。
进一步地，依据炮点和接收点的坐标数据确定有炮点区范围，将属于有炮点区的接收点打上标记1，属于无炮点区的接收点打上标记0。
进一步地，接收点的沙丘低速带校正量(tsd+tsdb)为接收点沙丘时间tsd与沙丘底部面下低速带静校正量tsdb之和，由下式计算：
tsd+tsdb＝-tsta+(Edatum-Esdb)/vR (2)
其中，Esdb为沙丘底部面高程，Edatum为最终基准面高程，vR为替换速度，tsd为接收点沙丘时间，tsdb为沙丘底部面下低速带静校正量，tsta为层析静校正量。
进一步地，接收点的沙丘厚度hsd等于接收点的地表高程减去沙丘底部面高程：
hsd＝Esurf-Esdb。
其中，Esurf为接收点的地表高程，Esdb为沙丘底部面高程。
进一步地，计算接收点的沙丘底部面下低速带静校正量包括：
将属于有炮点区接收点的坐标范围按给定网格边长分割成若干正方形网格单元；
计算每个网格单元中心点的沙丘底部面下低速带静校正量tsdb；
判断各网格单元tsdb的有效性；
利用有效网格单元中心点的tsdb，内插外推计算所有接收点的tsdb。
进一步地，利用与该网格单元中心点距离小于等于给定半径范围内的所有属于有炮点区的接收点中沙丘厚度hsd小于给定的拟合半长度值的沙丘低速带校正量(tsd+tsdb)与沙丘厚度hsd数据，拟合出一条直线，该直线截距为该网格单元中心点的tsdb；
统计所有网格单元参与tsdb计算的平均接收点数，如果某网格单元参与tsdb计算的接收点数小于平均接收点数的预定比例，该网格单元的tsdb被判定为无效。
进一步地，沙丘低速带校正量(tsd+tsdb)减去沙丘底部面下低速带静校正量tsdb得沙丘时间tsd。
进一步地，利用所有属于有炮点区接收点的沙丘时间tsd、沙丘厚度hsd数据，采用逐样点拟合的方法计算拟沙丘曲线函数tsd＝f(hsd)。
进一步地，利用拟沙丘曲线函数tsd＝f(hsd)和各接收点沙丘底部面下低速带静校正量tsdb，由公式(3)计算无炮点区接收点的静校正量：
tsta＝-f(hsd)-tsdb+(Edatum-Esdb)/vR (3)
其中，tsta为静校正量，Esdb为沙丘底部面高程，Edatum为最终基准面高程，vR为替换速度，tsdb为沙丘底部面下低速带静校正量，hsd为沙丘厚度。
进一步地，计算的无炮点区接收点的静校正量替换原先输入的无炮点区接收点的层析静校正量，和原先输入的有炮点区的接收点层析静校正量一起构成解决了边界问题的接收点静校正量。
根据本发明的另一方面，提供一种解决沙漠探区层析静校正边界问题的系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
依据炮点和接收点的坐标数据将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；
最后输出接收点静校正量。
本发明的方法应用于新疆沙漠地区多块三维地震勘探工区的资料处理，都很好地解决了层析静校正存在的边界问题，效果非常明显。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本发明实施例的沙漠区近地表速度模型的简化表达图。
图2示出了本发明实施例的工区炮点(深色)和接收点(浅色)位置图。
图3示出了本发明实施例的工区接收点高程图。
图4示出了本发明实施例的工区沙丘底部面高程图。
图5示出了本发明实施例的工区接收点层析静校正量图。
图6示出了本发明实施例的属于有炮点区(深色)和属于无炮点区的接收点(浅色)。
图7示出了本发明实施例的工区沙丘厚度图。
图8示出了本发明实施例的属于有炮点区接收点的沙丘低速带校正量图。
图9示出了本发明实施例的接收点沙丘底部面下低速带静校正量图。
图10示出了本发明实施例的属于有炮点区接收点沙丘时间图。
图11示出了本发明实施例的属于有炮点区接收点的沙丘时间与沙丘厚度的关系图。
图12示出了本发明实施例的工区统一的拟沙丘曲线函数表。
图13示出了本发明实施例的属于无炮点区的接收点原静校正量被新静校正量替换后的工区接收点静校正量。
图14示出了本发明实施例的涉及无炮点区的炮道集应用经本发明方法处理前(上)后(下)的静校正量的效果对比图。
图15示出了本发明实施例的应用经本发明方法处理前(左)后(右)的静校正量的涉及无炮点区的叠加效果对比图。
图16示出了本发明实施例的方法流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
在新疆塔里木等沙漠地区，地表高程除沙丘起伏引起的变化比较明显外，沙丘底部面变化相对平缓，沙丘底部面到高速层顶界之间厚度和速度变化相对也比较平稳，沙丘时间和沙丘厚度之间的关系符合一定的规律，可以用一条沙丘曲线大致代表工区沙丘时间与沙丘厚度之间的关系。这些近地表特点使得沙丘曲线法成为解决沙漠区静校正问题的一种有效方法。
沙丘曲线来自对野外近地表调查资料的分析与统计。在剩余静校正技术的配合下，沙丘曲线法可以获得较好的成像效果。沙丘曲线法可能隐含的风险，是没有考虑沙丘底部面到高速顶界面之间的低速带速度和厚度的变化，对低幅度构造的识别存在风险。
由于沙丘曲线法可以有效解决沙丘部分的时间校正问题，因此可以应用有炮点区初至层析计算的接收点静校正量，统计出一条“拟沙丘曲线”函数和平面上相对稳定的沙丘底部面到高速顶界间的剥离与填充时间，工区四周无炮点区接收点静校正量用“拟沙丘曲线”函数和沙丘底部面到高速顶界间的剥离与填充时间来计算。
本发明披露了一种“拟沙丘曲线法”，采用初至层析静校正技术获得的属于有炮点区接收点的静校正量，统计出一条“拟沙丘曲线”函数，用于计算属于无炮点区接收点的静校正量，解决初至层析静校正技术存在的接收点静校正边界问题。
沙漠区近地表速度模型可用图1简化表示，应用式(1)计算静校正量：
Figure BDA0001422895150000061
式中，Esurf、Esdb和Ebol分别为计算点处的地表、沙丘底部面高程和低速带底界的高程，
Figure BDA0001422895150000062
与
Figure BDA0001422895150000063
分别为计算点处的沙丘平均速度(沙丘厚度的函数)和沙丘底部面到高速顶界的平均速度(随平面位置缓慢变化)，Edatum与vR分别为最终基准面高程与替换速度(常数)。右侧第一项tsd是沙丘时间；第2项是沙丘底部面到高速顶界间的剥离与填充时间tsdb，后文称其为“沙丘底部面下低速带静校正量”，是随平面位置平缓变化的量；第3项是沙丘底部面到最终基准面的填充时间。
在与式(1)相关的参数中，层析静校正量tsta、地表高程、沙丘底部面高程、最终基准面高程与替换速度已知，其它参数均未知。将tsd+tsdb移到等式左边，
tsd+tsdb＝-tsta+(Edatum-Esdb)/vR (2)
式(2)的右边均已知。(tsd+tsdb)是沙丘时间与沙丘底部面下低速带静校正量之和，代表沙丘与低速带的校正量，后文称其为“沙丘低速带校正量”。由式(2)从有炮点区接收点静校正量中计算沙丘低速带校正量(tsd+tsdb)，进一步从中分解计算出沙丘底部面下低速带静校正量并内插外推至全区所有接收点，进一步得到各接收点的沙丘时间。用有炮点区所有接收点沙丘时间和沙丘厚度统计出全区统一的以列表方式表达的拟沙丘曲线函数tsd＝f(hsd)(沙丘时间tsd与沙丘厚度hsd间的关系，hsd＝Esurf-Esdb)，无炮点区接收点静校正量用该拟沙丘曲线函数和接收点的沙丘底部面下低速带静校正量计算：
tsta＝-f(hsd)-tsdb+(Edatum-Esdb)/vR (3)
依据炮点和接收点坐标，确定有炮点区范围，将属于有炮点区的接收点与属于无炮点区的接收点分别打上各自的标记，比如属于有炮点区接收点标记1，属于无炮点区接收点标记0。
用(2)式计算所有接收点的沙丘低速带校正量(tsd+tsdb)，它是接收点沙丘时间tsd与沙丘底部面下低速带静校正量tsdb之和。
从沙丘低速带校正量(tsd+tsdb)中分解出沙丘底部面下低速带静校正量tsdb与沙丘时间tsd。先计算各接收点沙丘底部面下低速带静校正量tsdb，再从(tsd+tsdb)中减去tsdb即得沙丘时间。
从沙丘低速带校正量(tsd+tsdb)中分解出沙丘底部面下低速带静校正量包含下列4步：
第1步，将有炮点区接收点坐标范围按给定单元边长(比如1000m)分割成若干正方形网格单元；
第2步，计算每个网格单元中心点的tsdb，利用与该网格单元中心点距离小于给定半径(比如9000m，越大平均效应越强)的所有属于有炮点区的接收点中沙丘厚度hsd小于给定值(拟合半长度(厚度)，比如10m)的接收点的沙丘低速带校正量(tsd+tsdb)与沙丘厚度hsd数据(沙丘低速带校正量(tsd+tsdb)为纵坐标，沙丘厚度为横坐标)，拟合出一条直线，该直线截距(沙丘厚度0处时间)即是该网格单元中心点的tsdb；
第3步，判断各网格单元tsdb的有效性，具体方法是先统计所有网格单元参与tsdb计算的平均接收点数，如果某网格单元参与tsdb计算的接收点数小于平均接收点数的某一个比例(比如2/3)，该网格单元的tsdb被判定为无效；
第4步，利用有效网格单元中心点的tsdb，内插外推计算所有接收点的tsdb。
利用有炮点区内所有接收点沙丘时间tsd与沙丘厚度hsd数据计算一条拟沙丘曲线函数，即工区统一的拟沙丘曲线tsd＝f(hsd)。
计算拟沙丘曲线函数的过程，实际上是从由接收点的沙丘时间tsd、沙丘厚度hsd构成的散点数据中找出它们之间用列表方式表达的函数关系。尽管可以采用多项式拟合的方法计算拟沙丘曲线函数，实际资料试验结果表明采用逐样点拟合的方法实现更加合理可靠。假设有N个接收点的沙丘时间、沙丘厚度数据对可用，即有N个接收点数据，其中最大厚度hmax，计算拟沙丘曲线函数的采样间隔dz。那么整个拟沙丘曲线的采样点数为M＝hmax/dz+1.5取整，第i个采样点的厚度为i*dz，i＝0，1，2，…，M-1。又假设拟合半长度(厚度)参数为FL2，那么参与厚度位于z的拟沙丘曲线函数上的采样点计算所用到的接收点沙丘数据的厚度范围为[z-FL2，z+FL2]，采用线性拟合的方法计算参与该采样点计算的接收点的沙丘时间、沙丘厚度数据对的线性表达式，由该线性表达式计算的该样点沙丘厚度z处的沙丘时间即该采样点对应的沙丘时间，该采样点的沙丘时间、沙丘厚度对就是要计算的拟沙丘曲线函数的一对数据。M个采样点的沙丘时间、沙丘厚度对构成了以列表方式表达的拟沙丘曲线函数。
利用工区统一的拟沙丘曲线函数tsd＝f(hsd)和各接收点沙丘底部面下低速带静校正量tsdb，由(3)式计算无炮点区的接收点的静校正量，并替换原先的层析静校正量。由于拟沙丘曲线函数是以列表方式表达，在应用(3)式计算接收点静校正量时，对沙丘厚度不刚好在沙丘曲线函数采样点上的情况使用内插外推的方法。
计算的无炮点区接收点的静校正量和原先输入的有炮点区的接收点层析静校正量一起构成解决了边界问题的接收点静校正量。
上述所述线性拟合，是指针对一组成对的数据(x1,y1,x2,y2,......,xN,yN)，N为总点数，找出一个y＝ax+b的线性方程，使
Figure BDA0001422895150000091
达到最小。针对本发明技术，沙丘厚度是该成对数据中的x，沙丘时间是该成对数据中的y。
根据本发明的一个方面，提供一种解决沙漠探区层析静校正边界问题的方法，该方法包括：
将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点的沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；
最后输出接收点静校正量。
具体地，参照图16，本发明的方法实现步骤如下：
⑴数据与参数输入
首先输入工区炮点平面坐标，接收点平面坐标和地表高程、沙丘底部面高程、接收点层析静校正量、静校正替换速度和基准面高程等。
⑵将接收点分成属于有炮点区和属于无炮点区两部分
依据炮点和接收点坐标，确定有炮点区范围，将属于有炮点区的接收点打上标记1，属于无炮点区的接收点打上标记0。
⑶计算沙丘低速带校正量和沙丘厚度
用式(2)计算各接收点的沙丘低速带校正量(tsd+tsdb)，它是接收点沙丘时间tsd与沙丘底部面下低速带静校正量tsdb之和。
tsd+tsdb＝-tsta+(Edatum-Esdb)/vR (2)
接收点的沙丘厚度hsd等于接收点地表高程减去沙丘底部面高程hsd＝Esurf-Esdb。
⑷计算接收点沙丘底部面下低速带静校正量
计算各接收点沙丘底部面下低速带静校正量tsdb包含下列4步。
第1步，将有炮点区接收点坐标范围按给定网格边长(该参数命名为GL，比如1000m)分割成若干正方形网格单元；
第2步，计算每个网格单元中心点的tsdb，具体方法是利用与该网格单元中心点距离小于等于给定半径(该参数命名为RD，比如9000m)范围内的所有属于有炮点区的接收点中沙丘厚度hsd小于给定的拟合半长度(厚度)值(该参数命名为FL2，比如10m)的沙丘低速带校正量(tsd+tsdb)与沙丘厚度hsd数据，拟合出一条直线，该直线截距(沙丘厚度0处时间)即是该网格单元中心点的tsdb；
第3步，判断各网格单元tsdb的有效性，具体方法是先统计所有网格单元参与tsdb计算的平均接收点数，如果某网格单元参与tsdb计算的接收点数小于平均接收点数的某一个比例(该参数命名为VR，比如67％)，该网格单元的tsdb被判定为无效；
第4步，利用有效网格单元中心点的tsdb，内插外推计算所有接收点的tsdb。
⑸计算有炮点区接收点的沙丘时间
有炮点区内各接收点的沙丘低速带校正量(tsd+tsdb)减去各自接收点的沙丘底部面下低速带静校正量tsdb，得各自接收点的沙丘时间tsd。
⑹计算拟沙丘曲线函数
利用有炮点区内所有接收点沙丘时间tsd与沙丘厚度hsd数据计算出一个以列表方式表达的拟沙丘曲线函数，即工区统一的拟沙丘曲线函数tsd＝f(hsd)。
利用所有属于有炮点区接收点的沙丘时间tsd、沙丘厚度hsd数据，采用逐样点拟合的方法计算拟沙丘曲线函数。假设有N个接收点的沙丘时间、沙丘厚度数据对可用，即有N个接收点数据，其中最大厚度hmax，计算拟沙丘曲线函数的采样间隔dz。那么整个拟沙丘曲线的采样点数为M＝hmax/dz+1.5取整，第i个采样点的厚度为i*dz，i＝0，1，2，…，M-1。拟合半长度(厚度)参数取与步骤⑷中使用的相同的拟合半长度参数(即FL2)，那么参与厚度位于z的拟沙丘曲线函数上的采样点计算所用到的接收点沙丘数据的厚度范围为[z-FL2，z+FL2]，采用线性拟合的方法计算参与该采样点计算的接收点的沙丘时间、沙丘厚度数据对的线性表达式，由该线性表达式计算的该样点沙丘厚度z处的沙丘时间即该采样点对应的沙丘时间，该采样点的沙丘时间、沙丘厚度对就是要计算的拟沙丘曲线函数的一对数据。M个采样点的沙丘时间、沙丘厚度对构成了以列表方式表达的拟沙丘曲线函数。
⑺计算无炮点区接收点静校正量
利用工区统一的拟沙丘曲线函数tsd＝f(hsd)和各接收点沙丘底部面下低速带静校正量tsdb，由(3)式计算无炮点区接收点静校正量，并替换原先的层析静校正量。由于拟沙丘曲线函数是以列表方式表达，在应用(3)式计算接收点静校正量时，对沙丘厚度不刚好在沙丘曲线函数采样点上的情况采用内插外推的方法。
tsta＝-f(hsd)-tsdb+(Edatum-Esdb)/vR (3)
根据本发明的另一实施方式，提供一种解决沙漠探区层析静校正边界问题的系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
依据炮点和接收点的坐标数据将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区接收点的静校正量并替换原先的层析静校正量；
输出接收点静校正量。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
本实施例以一个新疆沙漠三维工区实际资料作为例子。
图2为工区炮点和接收点相对位置，接收点均匀地布置在整个工区，而炮点只在工区内部布置。图3为工区接收点高程，可见到沙丘的起伏。图4为工区沙丘底部面高程，相对平缓。图5为一个当前主流的层析静校正软件计算的工区接收点静校正量，结合图2，可见属于无炮点区的接收点静校正量存在明显的边界问题。因为静校正量与地表高程之间应该呈现明显的负相关性，图5中只在属于有炮点区的接收点符合这种负相关性。工区静校正计算时使用基准面高程1100m，替换速度2000m/s。
图6是执行“将接收点分成属于有炮点区和属于无炮点区两部分”步骤时，将接收点分成属于有炮点区和属于无炮点区的结果。
图7和图8分别是执行“计算沙丘低速带校正量和沙丘厚度”，得到的工区沙丘厚度和属于有炮点区接收点的沙丘低速带校正量。
图9是执行“计算接收点沙丘底部面下低速带静校正量”，得到的接收点沙丘底部面下低速带静校正量，该工区该静校正量是平面上变化非常平缓的量。计算过程GL参数取1000m，RD参数取9000m，FL2参数取10m，VR参数取67％。
图10是执行“计算有炮点区接收点的沙丘时间”，得到的属于有炮点区接收点的沙丘时间平面图。
图11是有炮点区接收点沙丘时间和沙丘厚度的关系图，其中的浅灰色方点是拟合得到的拟沙丘曲线函数上的点(即图12的拟沙丘曲线表数据的散点图)。
图12是执行“计算拟沙丘曲线函数”，得到的工区统一的拟沙丘时间曲线数据表，计算拟沙丘曲线函数时使用的拟合半长度(厚度)参数FL2取10m，拟沙丘曲线采样间隔dz＝1，沙丘厚度最大值39m。
图13是属于无炮点区接收点的原静校正量被替换后的工区接收点静校正量。用于替换的属于无炮点区接收点的新静校正量是执行“计算无炮点区接收点静校正量”得到的结果。与图5相比，图13在工区中部有炮点区没有变化，而在工区四周无炮点区有明显改变，改变后的静校正量与图3的地表高程之间表现出明显的负相关性。
图14是涉及无炮点区的炮道集应用经本发明方法处理前(上)后(下)的静校正量的效果对比。图中‘B’点左侧的接收点属于有炮点区，本方法处理前后没有变化。但在‘B’点右侧的接收点属于无炮点区，层析静校正存在明显的边界问题，应用层析静校正后的道集(上)初至波同相轴存在明显的起伏。经本方法处理后的接收点静校正边界问题得到较好解决，应用后道集(下)初至波同相轴的起伏得以消除。
图15应用经本发明方法处理前(左)后(右)的静校正量的涉及无炮点区的叠加效果对比。应用本方法后明显提高了涉及无炮点区接收点的边界区域的叠加效果。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (9)

1.一种解决沙漠探区层析静校正边界问题的方法，其特征在于，该方法包括：
将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点的沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；
输出接收点静校正量；
其中，接收点的沙丘厚度hsd等于接收点的地表高程减去沙丘底部面高程：
hsd＝Esurf-Esdb
其中，Esurf为接收点的地表高程，Esdb为沙丘底部面高程；
其中，接收点的沙丘时间利用沙丘低速带校正量(tsd+tsdb)减去沙丘底部面下低速带静校正量tsdb来计算；
其中，利用拟沙丘曲线函数tsd＝f(hsd)和各接收点沙丘底部面下低速带静校正量tsdb，由公式(3)计算无炮点区接收点的静校正量：
tsta＝-f(hsd)-tsdb+(Edatum-Esdb)/vR (3)
其中，tsta为计算得到的接收点静校正量，Esdb为沙丘底部面高程，Edatum为最终基准面高程，vR为替换速度，tsdb为沙丘底部面下低速带静校正量，hsd为沙丘厚度。
2.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，输入数据与参数，包括：工区炮点平面坐标、接收点平面坐标、地表高程、沙丘底部面高程、接收点层析静校正量、静校正替换速度和基准面高程。
3.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，依据炮点和接收点的坐标数据确定有炮点区范围，将属于有炮点区的接收点打上标记1，属于无炮点区的接收点打上标记0。
4.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，接收点的沙丘低速带校正量(tsd+tsdb)为接收点沙丘时间tsd与沙丘底部面下低速带静校正量tsdb之和，由下式计算得到：
tsd+tsdb＝-tsta+(Edatum-Esdb)/vR (2)
其中，Esdb为沙丘底部面高程，Edatum为最终基准面高程，vR为替换速度，tsd为接收点沙丘时间，tsdb为沙丘底部面下低速带静校正量，tsta为层析静校正量。
5.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，计算接收点的沙丘底部面下低速带静校正量包括：
将属于有炮点区接收点的坐标范围按给定网格边长分割成若干正方形网格单元；
计算每个网格单元中心点的沙丘底部面下低速带静校正量tsdb；
判断各网格单元tsdb的有效性；
利用有效网格单元中心点的tsdb，内插外推计算所有接收点的tsdb。
6.根据权利要求5所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，利用与该网格单元中心点距离小于等于给定半径范围内的所有属于有炮点区的接收点中沙丘厚度hsd小于给定的拟合半长度值的沙丘低速带校正量(tsd+tsdb)与沙丘厚度hsd数据，拟合出一条直线，该直线截距为该网格单元中心点的tsdb；
统计所有网格单元参与tsdb计算的平均接收点数，如果某网格单元参与tsdb计算的接收点数小于平均接收点数的预定比例，该网格单元的tsdb被判定为无效。
7.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，利用所有属于有炮点区接收点的沙丘时间tsd、沙丘厚度hsd数据，采用逐样点拟合的方法计算拟沙丘曲线函数tsd＝f(hsd)。
8.根据权利要求1所述的解决沙漠探区层析静校正边界问题的方法，其特征在于，计算的无炮点区接收点的静校正量替换原先输入的无炮点区接收点的层析静校正量，和原先输入的有炮点区的接收点层析静校正量一起构成解决了边界问题的接收点静校正量。
9.一种解决沙漠探区层析静校正边界问题的系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
依据炮点和接收点的坐标数据将接收点分成属于有炮点区和属于无炮点区两部分；
计算接收点的沙丘低速带校正量和沙丘厚度；
计算接收点沙丘底部面下低速带静校正量和沙丘时间；
利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；
利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；
输出接收点静校正量；
其中，接收点的沙丘厚度hsd等于接收点的地表高程减去沙丘底部面高程：
hsd＝Esurf-Esdb
其中，Esurf为接收点的地表高程，Esdb为沙丘底部面高程；
其中，沙丘低速带校正量(tsd+tsdb)减去沙丘底部面下低速带静校正量tsdb得沙丘时间tsd；
其中，利用拟沙丘曲线函数tsd＝f(hsd)和各接收点沙丘底部面下低速带静校正量tsdb，由公式(3)计算无炮点区接收点的静校正量：
tsta＝-f(hsd)-tsdb+(Edatum-Esdb)/vR (3)
其中，tsta为计算得到的接收点静校正量，Esdb为沙丘底部面高程，Edatum为最终基准面高程，vR为替换速度，tsdb为沙丘底部面下低速带静校正量，hsd为沙丘厚度。
CN201710899715.9A 2017-09-28 2017-09-28 一种解决沙漠探区层析静校正边界问题的方法及系统 Active CN109581503B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710899715.9A CN109581503B (zh) 2017-09-28 2017-09-28 一种解决沙漠探区层析静校正边界问题的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710899715.9A CN109581503B (zh) 2017-09-28 2017-09-28 一种解决沙漠探区层析静校正边界问题的方法及系统

Publications (2)

Publication Number Publication Date
CN109581503A CN109581503A (zh) 2019-04-05
CN109581503B true CN109581503B (zh) 2020-08-25

Family

ID=65914500

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710899715.9A Active CN109581503B (zh) 2017-09-28 2017-09-28 一种解决沙漠探区层析静校正边界问题的方法及系统

Country Status (1)

Country Link
CN (1) CN109581503B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112099091B (zh) * 2019-06-17 2024-03-01 中国石油天然气集团有限公司 地震勘探点数据静校正方法及装置
CN112764107B (zh) * 2019-11-04 2024-05-28 中国石油天然气集团有限公司 地表高速覆盖地区的地震资料静校正方法及装置
CN111090118A (zh) * 2020-01-06 2020-05-01 中国石油化工股份有限公司 基于边界炮控制的层析静校正方法
CN112379434B (zh) * 2020-10-30 2024-04-30 中国石油天然气集团有限公司 适用于沙漠区的层析反演静校正的方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393270A (zh) * 2007-09-21 2009-03-25 中国石油集团东方地球物理勘探有限责任公司 一种连续介质地区的表层高精度静校正方法
CN107179553A (zh) * 2017-07-27 2017-09-19 中国石油化工股份有限公司 基于双压实规律的沙漠区表层静校正方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5587968A (en) * 1995-08-25 1996-12-24 Western Atlas International, Inc. Method for measuring the near-surface shear wave velocity for use in determining 3-component 3-D statics
RU2008133896A (ru) * 2008-08-18 2010-02-27 Государственное образовательное учреждение высшего профессионального образования "Саратовский государственный университет им. Н.Г. Ч Способ определения относительных статических поправок в сейсморазведке
US8416644B2 (en) * 2010-06-29 2013-04-09 King Abdul Aziz City for Science and Technology (KACST) System for acquiring seismic data over sand dune surfaces
CN102590864A (zh) * 2011-12-31 2012-07-18 中国石油集团西北地质研究所 两步法层析反演近地表建模方法
CN103293554B (zh) * 2013-05-07 2014-12-24 东北石油大学 一种基于地震数据的静校正方法
CN104656139B (zh) * 2013-11-21 2018-04-06 中国石油天然气集团公司 一种转换波的静校正方法及设备
US11635539B2 (en) * 2015-04-30 2023-04-25 Saudi Arabian Oil Company Imaging shallow heterogeneities based on near-surface scattered elastic waves

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393270A (zh) * 2007-09-21 2009-03-25 中国石油集团东方地球物理勘探有限责任公司 一种连续介质地区的表层高精度静校正方法
CN107179553A (zh) * 2017-07-27 2017-09-19 中国石油化工股份有限公司 基于双压实规律的沙漠区表层静校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"沙丘曲线静校正法在沙漠地区的应用及效果分析";武伟刚;《中国石油和化工标准与质量》;20121231(第1期);第260页 *

Also Published As

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

Similar Documents

Publication Publication Date Title
CN109581503B (zh) 一种解决沙漠探区层析静校正边界问题的方法及系统
CN106227923B (zh) 一种岩体结构面三维粗糙度评价方法
CN104200507B (zh) 一种三维点云点法向量估计方法
CN101393497A (zh) 基于双目立体视觉的多点触摸方法
CN102997875B (zh) 一种复杂螺旋曲面高精度测量不确定度的分析方法
CN104122587B (zh) 一种基于炮检距域的异常初至识别方法以及系统
CN103411533A (zh) 结构光自适应多次曝光法
CN102944896A (zh) 表层调查数据的模型法静校正方法
CN104330827B (zh) 一种地表模型静校正的处理方法
CN104089599B (zh) 用于接触式测头测量中提取二维轮廓的准形态学滤波方法
CN110533603A (zh) 一种基于熵与卷积神经网络的点云降噪方法
CN103315770A (zh) 基于三维超声图像的目标参数处理方法及系统
CN109163674B (zh) 一种面结构光自动化三维测量中传感器测量视点规划方法
CN103324784A (zh) 一种基于局部约束的网格模型碰撞处理方法
CN103728666A (zh) 一种确定静校正中间参考面替换速度的方法
CN106199705B (zh) 地震勘探采集数据的检测方法
CN106842302A (zh) 一种批量编辑初至的方法及装置
CN110153582A (zh) 焊接方案生成方法、装置以及焊接系统
CN105093280A (zh) 表层模型对地震数据影响的低频与高频成分的分解方法
CN107452059B (zh) 一种三维插值建模方法、设备及其存储设备
CN103473812B (zh) 一种基于tps的等深线追踪算法
CN108053481A (zh) 三维点云法向量的生成方法、装置和存储介质
CN106887003A (zh) 基于八邻域深度差的点云边缘提取方法
CN105572741B (zh) 一种计算3d高频静校正量的方法
CN111624655B (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