具体实施方式
以下,将参照附图来详细说明本发明的实施例。
可通过计算机软件模块实现本发明的表层调查数据的模型法静校正方法。在本申请中描述的各个步骤不限制为上述步骤,其中的一些步骤可被进一步拆分为更多的步骤,并且一些步骤可合并为较少的步骤。
本发明还提供一种实现表层调查数据的模型法静校正方法的装置。
图1是示出根据本发明的示例性实施例的表层调查数据的模型法静校正方法的流程图。在所述实施例中,假设在通用计算机系统上通过实现各步骤的功能的软件模块执行所述模型法静校正方法。
参照图1,在步骤S110,计算机系统加载二维或三维的表层调查数据。根据本发明的示例性实施例,计算机系统从存储有所述地震调查数据的信息存储介质或数据库加载二维或三维的表层调查数据。计算机系统可选择性地加载使用不同的表层资料采集方法获得的表层调查数据。采集中表层控制点的控制程度必须反映表层岩性的变化,才可反映出不同岩性中的速度变化趋势,从而建立准确的表层速度模型。
图2示出根据用户的选择大炮初至数据、小折射表层调查数据、微测井表层调查数据等的图形界面。图3示例性地示出加载的表层调查数据。
在步骤S120,计算机系统对加载的表层调查数据执行统一处理,其中,按照时深变化趋势进行统一分层,并计算每个控制点在各层的速度和厚度参数。
例如,对于微测井方法取得的表层调查数据,可使用以下方法进行分层并计算每个控制点在各层的速度和厚度参数:
对各层的点做直线回归拟合,求出回归系数a、b和平均标准差s。在点数大于3的情况下,可剔除个别与拟合线偏差大于平均标准差2.5倍的点,然后再作最后的直线回归拟合,得到最终的a和b值,b为截距时间T,l/a为速度值。假设有n层,则求出所述n个层的截距时间Ti和层速度值Vi。即,第i层的时深曲线为:Ti(h)=Ti+h/v,本层的厚度=本层的深度-上一层的深度。采用解析法求解下式,求取层深度Hi:
Ti+h/vi=Ti+1+h/Vi+1 ......式1
......式2
又例如,对于小折射资料采集方法取得的表层调查数据,对各层的点做直线回归拟合,求出回归系数a、b和平均标准差s。在点数大于3的情况下,可剔除个别与拟合线偏差大于平均标准差2.5倍的点,然后再作最后的直线回归拟合,得到最终的a、b值,其中,b为截距时间T,l/a×1000为速度值。如有n层,则求出所述n个层的截距时间和视速度值。
当处理单支的小折射表层调查数据时,按照式3仅对各层数据执行校正:
TiL*=TiL-T0L ......式3
其中,L代表小折射左支数据。
当处理双支的小折射表层调查数据时,如果左右支的层数一致,则由式4和式5求出该点的平均速度Vi:
Vi=2ViL×ViR/(ViL+ViR) ......式4
V0=(V0L+V0R)/2 ......式5
同时,仍然存在TiL*=TiL-T0L,TiR*=TiR-T0R。此时将T0作为0点进行校正。
最终,最终获得的平均截距时间Ti为
Ti=(TiL*-TiR*)/2 ......式6
在式4、式和式6中,L代表小折射左支数据,R代表小折射右支数据。
当双支层数不一致时(例:左支有2层,右支有3层)时,仍然按照式5求取直达波V0的速度;其它层通常的做法是,先求得高速层(最后一层)的速度(即取左、右支得最后一层速度值),此后用式4分别求得左支、右支第i层的Vi。对于左、右支剩下的层,比较左右各层的速度,对速度相对接近的层求平均速度和截距时间值,再剩下的其中一支的速度层单独作为一层。
可以按下式7和式8计算每层的厚度:
......式7
......式8
在步骤S130,计算机系统基于统一处理的表层调查数据建立二维或三维近地表模型。
如果在步骤S110加载了二维表层调查数据,则根据表层控制点,利用相关系数原理进行内插得到炮检点的表层结构。这里,可以使用的控制点类型包括:小折射控制点、微测井控制点、大折射控制点、相交测线处交点的表层成果。还可以通过人为插入控制点来解决大范围缺失控制点的情况。根据式9建立二维近地表模型:
WH0j=F H0j+DH0j×K0 ......式9
WHij=F Hij+DHij×Ki
其中,K为层间相关系数,其大小与层间相关度成正比,在致密表层K值较大,在疏松表层K值较小,K值范围是0~1。当K=0时,确定不相关;当K=1时,确定强相关。
图4示例性地示出根据本发明的示例性实施例建立的二维近地表模型。
如果在步骤S110加载了三维表层调查数据,则根据表层控制点、采用层间相关三角网格插值原理进行内插得到炮检点的表层结构,以构建三维近地表模型。图5示例性地示出根据本发明的示例性实施例建立的三维近地表模型。
三维近地表模型以式10的数学模型来表示:
......式10
其中,At是小折射控制点构成的大三角形的面积,An是被内插点与小折射控制点构成的小三角形的面积,P是被内插点的位置,V为内插的数值项。
在三维模型建立中,对炮检点相关系,还可以根据地表高程的变化情况对其进行编辑,目的是衰减地形波峰处的层间相关性,增强地形波谷处的层间相关性,其实现方法如下:
定义地形波峰处衰减系数Kr和地形波谷处增强系数Kv(二者一般在0-0.5之间)、最小、最大相关系数值(一般最小可为0,最大可为1)以及给定控制点密度。此后,可通过如下过程对炮检点相关系数进行编辑和控制:
首先对炮检点相关系数做最大和最小控制,再以P点为中心,在略大于控制点密度的区域内,根据下式统计其N个地表高程E
i的平均值
......式12
再根据式13计算实际高程E
i与平均高程
的差的绝对值的平均差
......式13
该点高程Ei与平均值E的差DE0为
......式14
如果 那么比例系数K为:
......式15
否则,K=0。
如果DE0>0,则有:
R[P]=R[P]*R[P]*(1-K*Kr) ......式16
R[P]为原相关系数,在波峰处衰减,否则
R[P]=R[P]*R[P]*(1-K*Kv) ......式17
R[P]为原相关系数,在波谷处增强。
在步骤S140,计算机系统对建立的近地表模型进行优化处理。
其中,对于在步骤S130建立的二维近地表模型,计算机系统采用双向均值内插方法对建立的二维近地表模型执行交点模型闭合。具体为,采用双向均值内插方法在交点上求出一个控制点C点。假设L1测线和L2测线相交于C点,L1测线的控制点为R1和R2,L2测线的控制点为R3和R4。根据层间相关系数法,由R1点、R2点求出L1线上C点的成果C1,由R3点、R4点求出L2线上C点的成果C2,然后在将双向的成果均值得到C点的成果作为L1、L2线共用的控制点,以此达到交点模型的闭合。
另一方面,对于在步骤S130建立的三维近地表模型,计算机系统采用中心区域平滑方法消除三维近地表模型的模型畸变,以使近地表模型的速度模型符合实际地表状况。
首先,根据平面速度模型中的速度异常区域,确定出需要调整的一个或多个速度区域。所述每个区域的形状由一组参与计算的表层调查点确定,这些确定调整区域形状的表层调查点称为边界点。将相邻两个边界点连接所构成的所有边称为调整区域的边界。如果一个控制点恰好在边界上,则称其为边界上的控制点。在平面内由所有边界点所组成的区域形状必须构成一个凸多边形。如果调整区域不是凸多边形,则会造成调整所得到的相关系数曲面产生突变。
然后,在每个确定需要调整的区域内确定调整点,并且将区域内其余的表层调查点确定为被调整点。以调整点为中心,通过样条插值的方法计算得到其它被调整点层间相关系数的大小,通过层间相关系数的改变,调整速度异常区域的速度,以此方法对每个需要调整的区域执行区域平滑处理,消除速度突变。
此后,在步骤S150,计算机系统通过时深曲线关联操作建立二维或三维近地表各层的速度、厚度模型,并计算基于表层调查数据的模型法静校正量。
由于不同的近地表建模技术都有其适应范围,不适应工区近地表条件的建模技术所建立的近地表模型是不可靠的。如果根据对大炮初至特征、微测井资料等信息的分析,认为工区近地表不存在稳定的折射层,或者存在高速折射屏蔽层、速度倒转或薄层构成的隐蔽层,又没有足够的微测井资料约束,那么利用折射静校正技术建立的近地表速度模型就不一定可信。同样的,在没有足够微测井资料约束的情况下,对浅层地表存在高速初至屏蔽层(如高速砾岩层、永冻层等)、速度倒转或薄层构成的隐蔽层时,层析静校正技术往往也无法建立可靠的近地表速度模型。
因此,本发明基于表层调查资料的模型法静校正技术充分利用了不同类型的表层调查资料,在统一的数据平台上,对表层调查数据进行统一解释,再对从数据源中解析出表层调查点数据,采用合适的层间相关系数值,并在这些表层调查点的控制下,计算出各层的速度和厚度,再通过交互方式反复求证,建立准确的近地表模型,克服了在大炮初至拾取困难的地区建立近地表模型困难的难题,解决了复杂地表的静校正问题。
尽管已参照优选实施例表示和描述了本发明,但本领域技术人员应该理解,在不脱离由权利要求限定的本发明的精神和范围的情况下,可以对这些实施例进行各种修改和变换。