【技术实现步骤摘要】
本专利技术涉及地质勘探，具体为一种面向构造控制型花岗岩体高波数阻抗地震成像方法。
技术介绍
1、逆时偏移(rtm)与其他波方程成像方法相同，可分为两个步骤来获取地下反射体的图像。第一步是从地面地震资料中外推两个波场：一个发在震源端位置，另一个发在检波点位置，第二步骤是将这些外推的波场做时间上交叉相关，并形成一个零滞后相关的图像。
2、尽管逆时偏移技术具有高精度和高保真的成像能力，但在当前普遍采用的互相相关成像条件下，将会出现低波数伪影，从而影响到成像图像的质量，尤其是对于速度模型中存在较强的速度梯度异常值或者较强的波阻抗反射界面，这是因为正演模拟的震源波场和逆时延拓的检波点波场都会在这些波阻抗反射界面上产生强烈的逆散射能量，这些逆散射能量满足互相关成像条件，因此，成像过程会产生较强的低波数伪影，从而干扰成像质量。
3、因此，为了在得到高精度，清晰成像结果的基础上，需要消除这些低波数伪影，基于此，本专利技术提出一种面向构造控制型花岗岩体高波数阻抗地震成像方法。
技术实现思路
1、(一)解决的技术问题
2、针对现有技术的不足，本专利技术提供了一种面向构造控制型花岗岩体高波数阻抗地震成像方法，解决了上述
技术介绍
提到的“在当前普遍采用的互相相关成像条件下，将会出现低波数伪影，从而影响到成像图像的质量以及成像过程会产生较强的低波数伪影扰成像质量”的问题。
3、(二)技术方案
4、为实现以上目的，本专利技术通过以下技术方案予以实现：一
5、步骤s1：获取输入参数；
6、步骤s2：根据步骤s1中获取的参数并且通过成像表达式生成炮点xs所对应的成像结果；
7、步骤s3：叠加各炮点所对应的成像结果，生成并输出偏移成像结果。
8、优选的，所述步骤s2中还包括：
9、步骤s201：在笛卡尔坐标系中，考虑声波介质，基于速度场和地震子波，并且通过带有震源的二维声波波动方程来表示构建的震源波场；
10、步骤s202：通过根据观测的地震记录和速度场，并且通过波动方程来表示构建的检波点波场；
11、步骤s203：通过数值离散化的波动方程来表示更新震源波场的策略以及震源的正向波场；
12、步骤s204：利用波场延拓算子求取进行逆时延拓，并且通过公式来表示检波点波场逆时延拓；
13、步骤s205：对于体积模量和密度参数化的各向同性声介质，由于模型参数的变化而引起的misfitfunctionχ的变化，可以得到高波数阻抗地震成像条件；
14、步骤s206：记录每一计算递推时刻的震源波场和检波点波场，在时间窗内，正、反向波场子波相位相同、二者互相关输出达到最大值时，然后选取高波数阻抗地震成像条件并且对多炮数据叠加减少噪音影响，可以得到面向构造控制型花岗岩体高波数阻抗地震成像方法公式。
15、优选的，所述步骤s1还包括
16、步骤s101：计划制定
17、根据研究区域和目标制定相关的计划；
18、步骤s102：设备布置
19、部署相关的地震仪器以及传感器；
20、步骤s103：数据采集
21、通过步骤s102中部署的设备对相关数据进行采集并记录；
22、步骤s104：数据处理
23、对步骤s103中获取的数据进行处理，其中包括去噪、时域分析和频域分析。
24、优选的，所述s1中所需要获取的参数包括p波速度参数场、观测炮记录。
25、优选的，所述步骤s201中具体的带有震源的二维声波波动方程公式为
26、
27、优选的，所述步骤s202中波动方程具体公式为
28、
29、优选的，所述步骤s203中数值离散化的波动方程具体公式为
30、
31、优选的，所述步骤s204中表示检波点波场逆时延拓的公式为
32、
33、优选的，所述步骤s206中得到面向构造控制型花岗岩体高波数阻抗地震成像方法公式为
34、iz(x)＝kz(x)＝kκ(x)+kρ(x)
35、所述步骤s205中高波数阻抗地震成像条件表达公式为
36、δχ＝∫kκ(x)δlogκ(x)+∫kρ(x)δlogρ(x)dx
37、优选的，所述步骤s2中成像表达式具体为如下
38、
39、(三)有益效果
40、本专利技术提供了一种面向构造控制型花岗岩体高波数阻抗地震成像方法。具备以下有益效果：
41、(1)该面向构造控制型花岗岩体高波数阻抗地震成像方法在使用时，能够解决传统互相关逆时偏移成像方法产生的低波数伪影问题，进而产生更加清晰的成像结果，开发面向构造控制型深部高温花岗岩体高波数阻抗地震成像技术，为复杂构造区的地震资料处理和后续地震解释工作提供高精度成像的理论依据和技术支持。
42、(2)该面向构造控制型花岗岩体高波数阻抗地震成像方法在使用时，通过采用高波数阻抗地震成像条件，所有的低波数伪影都被去除干净，使得高温岩体、所有断层等构造反射界面都能清晰成像，对比传统互相关逆时偏移成像结果，本专利技术能够解决低波数伪影问题，进而产生更加清晰的成像结果。
本文档来自技高网...【技术保护点】
1.一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：包括以下步骤：
2.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S2中还包括：
3.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S1还包括
4.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：S1中所需要获取的参数包括P波速度参数场、观测炮记录。
5.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S201中具体的带有震源的二维声波波动方程公式为
6.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S202中波动方程具体公式为
7.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S203中数值离散化的波动方程具体公式为
8.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成
9.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S206中得到面向构造控制型花岗岩体高波数阻抗地震成像方法公式为
10.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤S2中成像表达式具体为如下
...【技术特征摘要】
1.一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：包括以下步骤：
2.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤s2中还包括：
3.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤s1还包括
4.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：s1中所需要获取的参数包括p波速度参数场、观测炮记录。
5.根据权利要求1所述的一种面向构造控制型花岗岩体高波数阻抗地震成像方法，其特征在于：所述步骤s201中具体的带有震源的二维声波波动方程公式为
6.根据权利要求1所述的一种面向构造控制型...
【专利技术属性】
技术研发人员：李伟，王鹏飞，隋建才，张鑫，王扬州，张辉，刘朋，高建明，郭威，杨继东，黄建平，
申请(专利权)人：山东能源集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。