[发明专利]利用沙漠地表高程数据计算沙丘底部面高程的方法及系统

说明书

本发明提出了一种利用沙漠地表高程数据计算沙丘底部面高程的方法及系统，该方法包括：输入工区测点坐标与参数；计算地形趋势面；计算所述地形趋势面与沙丘谷底间的高程差函数；确定属于沙丘谷底测点的高程分界面；计算沙丘谷底趋势面以及沙丘底部面。与常规沙丘底部面的确定方法相比，该方法显著提高了效率，不受人为因素的影响，结果稳定可靠，得到的沙丘底部面就在沙丘谷底附近的下方，整体趋势平滑。该发明已经在中国西部沙漠区实际资料处理中得到广泛应用，取得了良好的效果。

技术领域

本发明属于地震勘探资料数据处理等技术领域，尤其涉及沙漠工区地震勘探静校正计算技术中的从沙漠地表高程数据中计算沙丘底部面高程的技术。

背景技术

沙漠工区地震勘探资料处理中，应用沙丘曲线法相关技术计算静校正量时需要沙丘底部面的高程数据。地震勘探测量阶段只能得到沙漠工区随沙丘起伏的地表高程，无法直接得到沙丘底部面的高程数据。沙丘底部面高程都是通过对地表高程数据的分析处理得到。现有技术是由技术人员人工选择一系列被认为是沙丘谷底的测量样点，利用内插外推的方式计算整个工区的沙丘底部面，效率较低，结果受人为因素及所选沙丘谷底的测量样点的影响。

发明内容

本发明提出一种依据沙漠工区测量得到的地表高程数据，快速准确地计算沙丘底部面高程的方法，达到提高沙丘底部面计算效率和避免人为因素对计算结果影响的目的。

根据本发明的一个方面，提供一种利用沙漠地表高程数据计算沙丘底部面高程的方法，该方法包括：

输入工区测点坐标与参数；

计算地形趋势面；

计算所述地形趋势面与沙丘谷底间的高程差函数；

确定属于沙丘谷底测点的高程分界面；

计算沙丘谷底趋势面以及沙丘底部面；

输出沙丘底部面。

进一步地，输入工区测点的平面与高程坐标(x_i,y_i,z_i)，i＝1,2,...,N；N为测点总数。

输入的参数包括：平滑半径参数R1、R2与R3，测点的最大高程大于谷底测点分界面高程的幅度t，以及高程差最大值剔除百分比α。

进一步地，以高程属性数据(x,y,z)为输入，采用逐点平面拟合的方法，计算地形趋势面(x,y,u)：

以测点(x_i,y_i)为中心，半径R1范围内的所有测点，以平面拟合方法构建u＝ax+by+c的平面方程，由该平面方程计算的该点高程即所述地形趋势面在该点的高程u_i；所有测点的地形趋势面高程构成地形趋势面；其中下标i是测点序号。

进一步地，计算所述地形趋势面与沙丘谷底间的高程差函数，包括：

构建地形趋势面与地表之间的高程差属性数据(x,y,u-z)，将其中z_i≥u_i即地表高程大于等于地形趋势面高程的测点设置成无效高程差，z_iu_i即地表高程小于地形趋势面高程的测点设置为有效高程差；

以高程差属性数据(x,y,u-z)为输入，移动平均计算平均高程差分布函数(x,y,v)：以测点(x_i,y_i)为中心，半径R2范围内的除被设置为无效高程差的点以外的所有测点都参与计算平均值，得到该点的平均高程差v_i；所有测点的平均高程差构成平均高程差分布函数；其中下标i是测点序号。

进一步地，确定属于沙丘谷底测点的高程分界面，包括：