共查询到20条相似文献，搜索用时 187 毫秒
1.
2.
3.
4.
5.
长86km、南北向横跨乌鲁木齐坳陷的深地震反射剖面，揭示了北天山山前地壳的薄皮构造特征．共深度点叠加剖面的石河子以南部分显示了天山北缘平行山体的第一和第二排背斜构造．与双程时间分别为2.5～3.0s和5.5～6.0ｓ的反射事件对应的滑脱构造，将地壳深部构造与地表逆断裂 褶皱构造联系在一起．玛纳斯断裂以铲形方式向下延伸，在2.5ｓ左右深度归并于第一滑脱面，向南与清水河断裂汇合．在5.5～6.0ｓ深度上为与玛纳斯下背斜相连的主滑脱面．它们最终汇集到准噶尔南缘断裂．石河子以北的坳陷沉积深度达12～14km．沿剖面的莫霍界面深度在准噶尔盆地为45ｋｍ 左右，往南加深至50ｋｍ．对该区域内的深地震测深剖面和布格重力异常资料的分析结果，与深反射剖面的地壳结构图象具有一致性．深地震反射剖面通过1906年玛纳斯7.7级地震宏观震中区，共深度点叠加剖面用于推断玛纳斯7.7级地震与北天山山前地壳构造之间的关系:玛纳斯地震属于一类褶皱地震，其发震构造是由准噶尔南缘断裂、清水河逆冲断裂、滑脱面和玛纳斯浅部断坡组成的断层系． 相似文献
6.
7.
8.
复杂近地表散射衰减所有的地面观测波场,形成半随机半相干的近地表强散射噪音背景,弥漫整个炮集,淹没深层反射信号,是导致地震资料极低信噪比的主要原因.如何研究和评价近地表散射强度一直是石油勘探未解决的问题,这与起伏地表的粗糙度、近地表速度横向变化和结构倾角分布密切相关.基于前期复杂近地表边界元法波动方程数值模拟研究,本文提出一种复杂近地表散射振幅矩阵方法来分析近地表散射强度.首先对复杂近地表结构进行边界元配置方法离散,根据边界积分方程生成矩阵方程.我们不求解该矩阵方程(涉及海量计算),只是利用矩阵分析技术来解析矩阵方程中的散射振幅系数矩阵,研究复杂近地表结构对不同频率波场的散射强度.该方法利用边界元对近地表结构几何特征的精确表征,研究起伏地表和非规则地质分界面对地震波传播的影响,由基本解及其在边界上的法向导数经过高斯数值积分计算得到的散射振幅系数矩阵,不仅描述了任意两点之间的相互影响,同时还刻画了边界形状特征的影响,为评价不同地质结构的散射强度提供了可能性.作为初步评价手段,我们采用矩阵元素总和与矩阵维数之比作为表征散射振幅系数矩阵散射特征的标量复杂系数,通过理论和实际模型测试,形成了一套行之有效、计算快速的近地表复杂性分析方法. 相似文献
9.
核－幔边界的动力学背景 总被引：1，自引：1，他引：1
根据传播矩阵方法，并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷，以求取核－幔边界的动力学背景．计算结果显示：１．所求的核－幔边界起伏图像与Ｈａｇｅｒ等根据格林函数方法所求得的核－幔边界的起伏在全球范围内基本相符．２．核－幔边界处的环型场仅在数量上降为地表处的环型场的１／８左右，而极型场较地表处的极型场的流动图像有显著变化，数值也增为地表处极型场的３倍左右． 相似文献
10.
位于北京断陷西北边界的黄庄—高丽营断层是区域新构造运动中的一条重要的隐伏断层,本研究利用伪随机可控震源及其高精度反射地震方法,调查了该断层在北京奥林匹克公园段的特征.研究工作首先从理论上分析了可控震源解码原理和抗噪声能力,然后利用三分量波动方程正演指导了黄庄—高丽营断层的野外地震观测,在城市强噪声环境下采集了高精度的地震数据,经过数据处理在该地区第一次得到地震偏移深度剖面.获得的地震偏移深度剖面清晰地显示了黄庄—高丽营断层的位置,地质解释结果为进一步明确剖面上断层两侧新生代地层的厚度提供了依据. 相似文献
11.
从可控源电磁法的基本原理出发,推导了基于电场矢量波动方程的三维边值问题,利用广义变分原理,把边值问题转换为变分问题,并引入散度条件,避免了伪解的出现,使有限元计算在理论上更加完备.在准静态近似条件下,把水平电偶极子在空中和大地的远区电场闭合表达式作为有限元计算中的区域外边界条件,解决了边界条件加载的困难;把应用于地震模拟中的伪delta函数引入到可控源电磁法中的三维有限元模拟中,消除了源点的奇异性,提高了方程组的稳定性.通过对均匀大地和层状介质模型的模拟,检验了程序的正确性,并对典型的地质体模型进行了数值模拟,分析了其变化规律. 相似文献
12.
Viscoacoustic prestack reverse time migration based on the optimal time-space domain high-order finite-difference method 总被引：3，自引：1，他引：2
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution. 相似文献
13.
ABSTRACT Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE. 相似文献
14.
以二维声波方程为模型,在时间域深入研究了全波形速度反演.全波形反演要解一个非线性的最小二乘问题,是一个极小化模拟数据与已知数据之间残量的过程.针对全波形反演易陷入局部极值的困难,本文提出了基于不同尺度的频率数据的"逐级反演"策略,即先基于低频尺度的波场信息进行反演,得出一个合理的初始模型,然后再利用其他不同尺度频率的波场进行反演,并且用前一尺度的迭代反演结果作为下一尺度反演的初始模型,这样逐级进行反演.文中详细阐述和推导了理论方法及公式,包括有限差分正演模拟、速度模型修正、梯度计算和算法描述,并以Marmousi复杂构造模型为例,进行了MPI并行全波形反演数值计算,得到了较好的反演结果,验证了方法的有效性和稳健性. 相似文献
15.
为研究瑞雷面波的形成机理及传播规律，促进瑞雷面波资料处理方法的发展，本文根据弹性波方程，采用交错网格有限差分数值求解算法，对浅层各向同性弹性介质进行了包括瑞雷面波和体波在内的全波场模拟. 提出了变系数吸收边界条件并将之应用于正演模拟，使边界条件的处理简单而高效，同时给出了角点的处理方法. 对工程勘察中常见的连续和层状介质模型进行了模拟，获得了更加接近实际情况的地震记录. 结合模拟记录，探讨了瑞雷面波的形成条件，同时讨论了震源埋深对面波能量的影响. 相似文献
16.
目前,对于可控源电磁法各向异性介质2.5维问题,主要采用一次场、二次场分离的方法消除场源奇异性并降低截断边界对计算区域的影响.该方法数值计算精度高,但是很难适用于复杂地形条件下的数值模拟.针对复杂地形问题,基于总场的有限元方法表现出一定的优越性,然而,这种方法存在场源奇异性问题和截断边界问题.本文采用基于总场计算的方法对带地形的可控源电磁法2.5维各向异性介质进行模拟研究,推导了考虑电导率和介电常数各向异性的2.5维控制方程;引入网格加密-收缩算法降低场源奇异性的影响范围,提升数值计算效率;引入行波分解边界条件降低截断边界的影响;提出任意采样反傅里叶变换方法,快速、高精度地计算出空间域电磁场分量.理论模型数值算例中:首先,验证了本文算法的有效性;其次,对任意各向异性倾角产生的可控源电磁响应规律进行研究;最后,采用山丘模型对各向异性介质电磁场的响应规律进行了模拟和分析. 相似文献
17.
地震自发破裂模拟是震源动力学研究的重要内容,了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用,大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征,为了更为灵活地刻画断层几何复杂性,我们建立断层模型的三角形网格离散方案,通过精确的解析解形式来计算断层各个单元之间的应力格林函数,联立滑动弱化摩擦准则和非奇异边界积分方程,对断层的自发破裂过程进行了模拟.在简单的平面断层模型下,将计算结果与前人的结果进行了对比,验证了方法的正确性与有效性.对于几种常见的复杂断层模型,例如弯折、阶跃、含障碍体断层等,我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明,非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层,同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景. 相似文献
18.
本文研究了三层介质中膨胀球模拟震源的应力、应变场以及由这个应力、应变场产生的电阻率变化量和地表视电阻率变化的关系。结果表明,对于一定强度的“震源”,在地电装置探测范围内的应力场随具有不同弹性参数的层位有显著的变化;而应变场却看不出有明显的突变;由地表视电阻率变化量的计算结果揭示了地电异常受台站介质电性结构与力学结构的综合影响。因此作为预报地震的地电阻率法一方面要寻找具有高电阻率——应变灵敏度岩层的台址,另一方面台址的视电阻率变化对各层真电阻率变化要具有好的响应。 相似文献
19.
本文提出了一种新的基于平均导数优化方法(average-derivative optimal method,简称ADM)的二维VTI介质qP波波动方程频率空间域二阶9点格式,这种新算法将二维VTI介质qP波波动方程中中心空间导数项的差分近似表示为正交方向上3个网格点的加权平均形式.通过最小二乘优化方法求取空间导数项和加速度项的加权优化系数从而使数值频散达到极小化,每个波长所需要的网格点数在1%的误差范围内仅为3.57个网格点数,而VTI介质常规9点差分格式在相同的误差范围内则需要约12个网格点数,新方法的计算精度明显提高.复杂BP2007 2D VTI海洋标准模型数值模拟结果也验证了本文VTI介质9点ADM算法的有效性和准确性. 相似文献
20.
Double-square-root one-way wave equation prestack tau migration in heterogeneous media 总被引：1，自引：0，他引：1
In this paper, source‐receiver migration based on the double‐square‐root one‐way wave equation is modified to operate in the two‐way vertical traveltime (τ) domain. This tau migration method includes reasonable treatment for media with lateral inhomogeneity. It is implemented by recursive wavefield extrapolation with a frequency‐wavenumber domain phase shift in a constant background medium, followed by a phase correction in the frequency‐space domain, which accommodates moderate lateral velocity variations. More advanced τ‐domain double‐square‐root wave propagators have been conceptually discussed in this paper for migration in media with stronger lateral velocity variations. To address the problems that the full 3D double‐square‐root equation prestack tau migration could meet in practical applications, we present a method for downward continuing common‐azimuth data, which is based on a stationary‐phase approximation of the full 3D migration operator in the theoretical frame of prestack tau migration of cross‐line constant offset data. Migrations of synthetic data sets show that our tau migration approach has good performance in strong contrast media. The real data example demonstrates that common‐azimuth prestack tau migration has improved the delineation of the geological structures and stratigraphic configurations in a complex fault area. Prestack tau migration has some inherent robust characteristics usually associated with prestack time migration. It follows a velocity‐independent anti‐aliasing criterion that generally leads to reduction of the computation cost for typical vertical velocity variations. Moreover, this τ‐domain source‐receiver migration method has features that could be of help to speed up the convergence of the velocity estimation. 相似文献