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

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

A Method and System for Solving the Boundary Problem of Tomographic Static Correction in Desert Exploration Area

The invention provides a method and system for solving the boundary problem of tomographic static correction in desert exploration area. The method includes: dividing the receiving point into two parts: the receiving point belongs to the area with gunpoint and the area without gunpoint; calculating the correction amount of low-speed band of sand dune at receiving point and the thickness of sand dune; calculating the static correction amount and dune time of low-speed band under the bottom surface of sand dune at receiving point; utilizing the receiving point belongs to the area with gunpoint. Dune time and dune thickness data, statistical work area unified quasi-dune curve function; using quasi-dune curve function to calculate static correction of receiving points in non-shot area. The method of the present invention is applied to data processing of multi-block three-dimensional seismic exploration areas in Xinjiang desert area, and the boundary problem of tomographic static correction is well solved, and the effect is very obvious.

【技术实现步骤摘要】
一种解决沙漠探区层析静校正边界问题的方法及系统
本专利技术属于地震勘探资料数据处理领域，尤其涉及沙漠工区地震勘探静校正计算技术。
技术介绍
基于初至时间的层析静校正技术是解决沙漠工区静校正问题的有效手段。初至层析近地表速度模型反演的精度取决于近道初至信息，在有近道初至信息的有炮点区反演的近地表速度模型可靠，相应地计算的静校正量也可靠。而在工区四周无炮点区，由于缺乏近道初至信息，反演的近地表速度模型精度明显下降，计算的接收点静校正量精度随之下降。所有基于初至层析的静校正技术都有这样的边界问题，只是程度上有差异。现有的解决方法多采用把有炮点区近地表速度模型或静校正量进行外推，外推的精度决定了静校正精度，工区四周无炮点区接收点静校正量的精度难以令人满意。
技术实现思路
本专利技术提出一种“拟沙丘曲线法”，用于解决接收点静校正边界问题。应用有炮点区初至层析获得的接收点静校正量，统计出一条“拟沙丘曲线”函数，用以计算工区四周无炮点区的接收点静校正量，快速有效地解决沙漠勘探工区初至层析接收点静校正量的边界问题。根据本专利技术的一个方面，提供一种解决沙漠探区层析静校正边界问题的方法，该方法包括：将接收点分成属于有炮点区和属于无炮点区两部分；计算接收点的沙丘低速带校正量和沙丘厚度；计算接收点的沙丘底部面下低速带静校正量和沙丘时间；利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；最后输出接收点静校正量。进一步地，输入数据与参数，包括：工区炮点平面坐标、接收点平面坐标、地表高程、沙丘底部面高程、接收点层析静校正量、静校正替换速度和基准面高程。进一步地，依据炮点和接收点的坐标数据确定有炮点区范围，将属于有炮点区的接收点打上标记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.一种解决沙漠探区层析静校正边界问题的方法，其特征在于，该方法包括：将接收点分成属于有炮点区和属于无炮点区两部分；计算接收点的沙丘低速带校正量和沙丘厚度；计算接收点的沙丘底部面下低速带静校正量和沙丘时间；利用属于有炮点区接收点的沙丘时间和沙丘厚度数据，统计工区统一的拟沙丘曲线函数；利用拟沙丘曲线函数计算属于无炮点区的接收点的静校正量，并替换原先的层析静校正量；输出接收点静校正量。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小于给定的拟合半长度值的沙丘...

【专利技术属性】
技术研发人员：林伯香，徐颖，吕秋玲，邵文潮，马方正，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11