共查询到20条相似文献，搜索用时 62 毫秒
1.
地震正演模拟是逆时偏移和全波形反演中的核心问题之一,因为它们都需要高效、高精度地模拟波场正向和反向传播。为了提高数值模拟的精度,人们广泛采用高阶有限差分方法,但是大多数方法仅在空间上具有更高的精度,在时间上只有二阶精度。首先系统介绍时空域高精度交错网格有限差分方法的基本原理,然后利用模型验证方法的有效性,结果表明:时空域高精度交错网格有限差分方法拥有比常规交错网格有限差分方法更低的数值频散。 相似文献
2.
通常工业界实现逆时偏移算法时采用有限差分数值方法模拟地震波场,波场模拟常常受稳定性条件限制,且易产生数值频散,成像精度降低.本文引入了一步法波场延拓方法,首先构建声波传播算子,借助Chebyshev多项式和Jacobi-Anger展开式近似传播算子中的e指数项,进而实现波场递推,该方法时间步长的选取不受稳定性条件限制而且不存在空间频散现象.本文将一步法波场延拓方法用于逆时偏移成像的波场模拟,并提出双缓冲区存储策略,在不增加计算量的前提下,大幅降低了逆时偏移方法的波场存储量.波场模拟和逆时偏移成像测试表明,本文提出的一步法波场延拓方法模拟地震波场精度高,消除了频散影响,可在较大时间步长的情况下实现高精度波场模拟;提出的基于一步法波场延拓的逆时偏移方法成像质量好;基于双缓冲区存储策略的逆时偏移成像方法存储成本低. 相似文献
3.
声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.为精确求解空间偏导数,相继发展了高阶差分格式优化方法和伪谱方法.近期,为更好地缓解数值频散,提出了时间-空间域有限差分方法,该方法采用了泰勒展开近似方法来确定有限差分格式系数,因而只能保证在一定的小范围内很好的拟合波场传播规律.为进一步压制数值频散效应,本文引入了时间-空间域特定波数点满足频散关系的方法,根据震源、波速和网格间距确定波数范围,同时考虑了多个传播角度,然后建立方程确定了相应的有限差分格式系数,使得差分系数能在更大范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性. 相似文献
4.
时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质. 相似文献
5.
有限差分法广泛应用于地震波数值模拟、成像和波形反演中,差分数值解的精度直接影响着地震成像和反演的效果.因为有限差分算子可以通过截断伪谱法的空间褶积序列得到,而截断窗函数的属性影响有限差分算子逼近微分算子的精度.具体地讲,窗函数的幅值响应的主瓣和旁瓣决定了有限差分算子逼近的精度,主瓣越窄,旁瓣衰减越大,则有限差分算子逼近微分算子的精度越高,更好地压制数值频散.基于此认识,本文提出了一种基于Chebyshev自褶积组合窗截断逼近的有限差分算子优化方法.Chebyshev自褶积组合窗的主瓣较窄,且旁瓣衰减大,其可通过只调节三个参数,更直观和可视化地控制主瓣和旁瓣的形状,改变有限差分算子逼近微分算子的精度;该窗函数截断逼近的有限差分算子不仅有较大的谱覆盖范围,而且精度误差波动较小,这表明低阶的差分算子可以达到高阶算子的精度,且逼近误差更稳定;从经济上来讲,将有效地减少模拟计算花费,提高计算效率. 相似文献
6.
压制数值频散,提高正演模拟精度,一直是有限差分正演模拟研究的重要内容.基于时空域频散关系的有限差分法,比基于空间域频散关系的传统有限差分法,模拟精度更高.时空域声波方程数值模拟,普遍采用常规十字交叉型高阶有限差分格式.而在频率-空间域,普遍采用旋转网格和常规网格混合的有限差分格式,有效提高了模拟精度和计算效率.本文将频率-空间域混合网格有限差分的思想引入到时空域,提出了时空域混合网格2 M+N型声波方程有限差分方法.首先推导出基于时空域频散关系的混合网格差分系数计算方法,然后进行频散分析、稳定性分析,并和传统高阶、时空域高阶有限差分法对比,结果表明:计算量相同时,新方法能有效压制数值频散,显著提高模拟精度;新方法相比传统2 M阶有限差分法,稳定性增强,与时空域2 M阶有限差分法稳定性基本相当.最后利用新方法进行均匀介质、层状介质、盐丘模型的数值模拟和盐丘模型的逆时偏移,模拟效果和成像质量进一步证实了该方法的有效性和普遍适用性. 相似文献
7.
叠前逆时偏移是对地下介质进行精确成像的方法之一。由于实际地下介质具有粘滞性,研究粘滞声波叠前逆时偏移具有一定的现实意义。逆时偏移的步骤之一是求解波动方程,对地震波场进行正向和反向外推,因此,精确、高效地求解波动方程对逆时偏移的成像效果和计算效率具有重要影响。本文中,我们利用基于优化时空域频散关系的高阶有限差分方法求解粘滞声波方程。频散分析和数值模拟的结果证明了优化时空域有限差分方法具有较高的精度,可以很好地压制数值频散。利用混合吸收边界条件处理边界反射,然后利用震源归一化互相关成像条件进行成像,并利用拉普拉斯滤波方法去除低频噪音。数值模型的测试结果显示,在考虑地下介质的粘滞性时,粘滞声波方程逆时偏移比声波方程逆时偏移具有更高的成像分辨率。另外,在进行波场外推的时候,采用自适应变长度的有限差分算子计算空间导数,在不影响求解精度的情况下,有效地提高了计算效率。 相似文献
8.
叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率. 相似文献
9.
声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性. 相似文献
10.
有限差分时间域正演是弹性波逆时偏移和全波形反演的基础,正演的计算精度也控制着偏移结果的准确性,若精度不高,则在偏移、反演后会带来假象.为了有效提高正演精度,本文结合窗函数优化方法,在窗函数截断伪谱法空间褶积序列以逼近有限差分算子的基础上,提出了一种基于Chebyshev窗的余弦调制模型,在原始Chebyshev窗的基础上引入了调制次数和调制范围,通过调节这两个参数可以人工可视化的调节截断误差,新的窗函数继承了Chebyshev窗的特点,在不明显降低截断谱范围的基础上明显降低了截断误差.本文针对不同正演阶数N,给出了一组经验调制系数,并通过数值模拟方法,对比了新方法、改进二项式窗和基于最小二乘优化方法的正演效果.结果表明,基于余弦调制的Chebyshev窗控制数值频散的能力更强,在大网格下可以得到更精确的正演结果.从经济角度分析,该方法减小了计算花费,提高了计算效率. 相似文献
11.
《Advances in water resources》2005,28(7):701-710
The patterns of temporal variations of precipitation (P), streamflow (SF) and baseflow (BF) as well as their nitrate-nitrogen (nitrate) concentrations (C) and loads (L) from a long-term record (28 years) in the Raccoon River, Iowa, were analyzed using variogram and spectral analyses. The daily P is random but scaling may exist in the daily SF and BF with a possible break point in the scaling at about 18 days and 45 days, respectively. The nitrate concentrations and loads are shown to have a half-year cycle while daily P, SF, and BF have a one-year cycle. Furthermore, there may be a low-frequency cycle of 6–8 years in C. The power spectra of C and L in both SF and BF exhibit fractal 1/f scaling with two characteristic frequencies of half-year and one-year, and are fitted well with the spectrum of the gamma distribution. The nitrate input to SF and BF at the Raccoon watershed seems likely to be a white noise process superimposed on another process with a half-year and one-year cycle. 相似文献
12.
13.
14.
燕郊等测点迁移优化与地磁观测研究 总被引：1，自引：2，他引：1
为优化地磁观测条件，开展了燕郊、夏垫等测点迁移工作；按照测眯迁移原则与实施技术方案，完整地收集并整理了地质构造与地球物理等方面的基本资料；进行了野外实地勘察，磁场梯度的测量，确定了新测点，在新老测点上进行了较长时间的地磁场对比观测；应用多种方法分析研究了地磁对比观测资料，结果表明，新老测点与有关测点的地磁变化具有良好的一致性，并得到了新老测点之间的地磁数据的按点差。 相似文献
15.
16.
17.
We rederive and generalize hyperbolic moveout formulae for the common-midpoint (CMP) gather and for the zero-offset (ZO) section that can be efficiently used for macro-model-independent reflection imaging in two-dimensional media. The hyperbolic moveout formulae for the common-midpoint gather are obtained from different Taylor series expansions of a particular parametric moveout surface defined in the multicoverage data space. Such a moveout surface involves three kinematic wave-field attributes of two hypothetical waves, which have to be determined by a coherency analysis. By using hyperbolic moveout curves in the CMP gather and in the ZO section one can determine these attributes in two steps. The relationships between the local shapes of the interfaces and the attributes of the hypothetical wave-fields attributes are considered by means of geometrical optics. The determination of these attributes allows to perform a macro-model-independent ZO simulation and a subsequent inversion. 相似文献
18.
Summary This paper discusses the influence of system stiffness on the dynamic instability of fault surfaces under laboratory conditions for a number of test modes. In conjunction with shear load stiffness, the normal load stiffness, often neglected, is shown to have a considerable effect on the stick-slip process —its presence or absence and its characteristics. Also appropriate stiffnesses are suggested for an earthquake sequence modeled as a growing dislocation. 相似文献
19.
Vindell Hsu Charles E. Helsley Eduard Berg David A. Novelo-Casanova 《Pure and Applied Geophysics》1984,122(6):878-893
A close correlation in spatial distribution of local seismic activity and energy release patterns before and after the 1979 Petatlan, Mexico earthquake suggests heterogeneity within the fault plane of this major low-angle thrust event associated with subduction along the Middle America Trench. A simple two-asperity model is proposed to account for the complexity. Foreshocks and aftershocks of the neighboring 1981 Playa Azul earthquake showed a similar pattern. As both events occurred at the junction of the Orozco Fracture Zone and the Middle America Trench, we speculate that the observed complex fault plane is caused by subduction of the rugged ocean floor of the Orozco Fracture Zone. Short-term precursory seismicity prior to the Petatlan earthquake can be explained by using the asperity model and migration of a slip front from the south-east to the north-west across the main shock source region. 相似文献
20.