共查询到20条相似文献,搜索用时 46 毫秒
1.
2.
3.
The optimization inversion method based on derivatives is an important inversion technique in seismic data processing, where
the key problem is how to compute the Jacobian matrix. The computational precision of the Jacobian matrix directly influences
the success of the optimization inversion method. Currently, most of the AVO (amplitude versus offset) inversions are based
on approximate expressions for the Zoeppritz equations to obtain the derivatives of the seismic wave reflection coefficients
(SWRCs) with respect to the stratum parameters. As a result, the computational precision and range of applications of these
AVO inversions are restricted undesirably. In order to improve the computational precision and to extend the range of applications
of AVO inversions, the partial derivative equations of the Zoeppritz equations are established, with respect to the ratios
of wave velocities and medium densities. By solving the partial derivative equations of the Zoeppritz equations accurately,
we obtained the partial derivative of SWRCs with respect to the ratios of seismic wave velocities and medium densities. With
the help of the chain rule for derivatives, the gradient of the SWRCs can be accurately computed. To better understand the
behavior of the gradient of the SWRCs, we plotted the partial derivative curves of the SWRCs, analyzed the characteristics
of these curves, and gained some new insight into the derivatives. Because only a linear system of equations is solved in
our method without adding any new restrictions, the new computational method has both high precision and a quick running speed;
it is not only suitable for small incident angles and weak reflection seismic waves but also for large incident angles and
strong reflection seismic waves. With the theoretical foundations established in the article, we can further study inversion
problems for layered stratum structures and we can further improve the computational speed and precision of the inversions. 相似文献
4.
5.
6.
Ray path of head waves with irregular interfaces 总被引:1,自引:0,他引:1
Head waves are usually considered to be the refracted waves propagating along flat interfaces with an underlying higher velocity.
However, the path that the rays travel along in media with irregular interfaces is not clear. Here we study the problem by
simulation using a new approach of the spectral-element method with some overlapped elements (SEMO) that can accurately evaluate
waves traveling along an irregular interface. Consequently, the head waves are separated from interface waves by a time window.
Thus, their energy and arrival time changes can be analyzed independently. These analyses demonstrate that, contrary to the
case for head waves propagating along a flat interface, there are two mechanisms for head waves traveling along an irregular
interface: a refraction mechanism and transmission mechanism. That is, the head waves may be refracted waves propagating along
the interface or transmitted waves induced by the waves propagating in the higher-velocity media. Such knowledge will be helpful
in constructing a more accurate inversion method, such as head wave travel-time tomography, and in obtaining a more accurate
model of subsurface structure which is very important for understanding the formation mechanism of some special areas, such
as the Tibetan Plateau. 相似文献
7.
8.
讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试. 相似文献
9.
10.
实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法. 相似文献
11.
过套管电阻率测井是通过测量套管壁电势实现测量地层的视电阻率,基于传输线方程理论,针对层状地层,给出了套管壁电势、电流对地层横向电阻导数的微分方程(称Jacobi矩阵微分方程)及边界条件;利用Jacobi矩阵微分方程边值问题导出了过套管电阻率测井反演地层参数的Jacobi矩阵系数的解析表示,利用Marquardt方法实现了过套管测井的地层电阻率反演;通过计算对Jacobi矩阵的特性进行了探讨,并获得了较快的计算速度(因为Jacobi矩阵是用解析解表示的),反演结果与地层模型取得了较好的逼近.本文实现了过套管电阻率测井地层参数的Jacobi系数矩阵的快速计算及地层电阻率反演,为进一步开展电阻率测井数据处理提供了理论依据和快速反演算法. 相似文献
12.
G. NEUMANN 《Geophysical Prospecting》1981,29(2):161-177
Lateral inhomogeneities generate fluctuations in the traveltime of seismic waves. By evaluation of these traveltime fluctuations from different source and receiver positions, lateral inhomogeneities can be located using a pseudo inverse matrix method (Aki, Christoffersson and Husebye 1977). The formulation of the problem is possible for transmitted waves as well as for reflected and refracted waves. In reflection seismics this method is of importance, if no reflections from the inhomogeneities themselves, but only reflections from lower boundaries can be observed. The basic assumptions for the mathematical formulation are (1) the average velocities and depths of the reflecting horizons are known already from standard processing methods, and (2) the traveltime residuals are due to lateral velocity changes between different reflectors or between reflectors and the surface. The area of the earth to be considered is divided into layers and the layers into rectangular blocks. The parallel displacement of a ray after passing a disturbed block is neglected, only the traveltime residual is taken into account. In this paper the method and its application to data obtained with two-dimensional models are described. 相似文献
13.
A review of the 6 × 6 anisotropic interface ray propagator matrix in Cartesian coordinates and within the framework of the Hamiltonian formalism shows that there is one unique propagator satisfying the symplectic property. This is essential, since the symplecticity furnishes an exact inverse, while an eigenvalue analysis indicates that the propagator may be arbitrarily ill-conditioned. As such, the symplectic interface propagator naturally connects to symplectic ray integration algorithms for smooth media, designed to maintain accuracy. Moreover, several ray invariants for smooth media remain invariant across interfaces. It is straightforward to derive expressions for the interface propagator, both explicit and implicit. Symplecticity is equivalent to the condition that the propagator preserves the eikonal constraint across the interface. The symplectic interface propagator complies with phase matching of the incident and reflected/transmitted ray field, and is therefore in accordance with the earlier derived 4 × 4 matrix in ray-centred coordinates. The symplectic property is related to the symmetry of the second derivative matrix of the reflected/transmitted traveltime field. Thanks to the analytic expression of the symplectic interface propagator, relating interface curvature directly to second derivatives of traveltimes observed at a datum level, numerous applications are available in the area of processing and inversion. 相似文献
14.
《中国地震研究》2021,(1)
Accurate calculations of travel times and raypaths of reflection waves are important for reflection travel time tomography. The multistage shortest path method(MSPM) and multistage fast marching method(MFMM) have been widely used in reflection wave raytracing, and both of them are characterized by high efficiency and accuracy. However, the MSPM does not strictly follow Snell's law at the interface because it treats the interface point as a sub-source, resulting in a decrease in accuracy. The MFMM achieves high accuracy by solving the Eikonal equation in local triangular mesh. However, the implementation process is complex. Here we propose a new method which uses linear interpolation to compute the incident travel time of interface points and then using Snell's law to compute the reflection travel time of grid points just above the interface. Our new method is much simpler than the MFMM; furthermore, numerical simulations show that the accuracy of the MFMM and our new method are basically the same, thus the reflection tomography algorithms which use our new method are easier to implement without decreasing accuracy. Besides, our new method can be extended easily to other grid-based raytracing methods. 相似文献
15.
—Within the "Integrated Seismics Oberpfalz 1989 (ISO89)" a three-component Moving Source Profiling (MSP) experiment, also named walk-away VSP, was carried out at the drilling site of the "Kontinentales Tiefbohrprogramm der Bundesrepublik Deutschland (KTB)" in Germany. Analysis of transmitted waves traveling from the source locations at the surface down to the receiver array in the borehole reveals velocity information about the illuminated part of the subsurface. Complementary to the widely used evaluation of travel-time perturbations to locate velocity inhomogeneities we suggest the use of the directivity of transmitted wave types down in the borehole. To determine the wave-field directivity we focus on transmitted arrivals by employing principles of "Controlled Directional Reception (CDR)." We calculate local slant-stacks for three different depth positions as a function of the source offset, thus obtaining the variation of the vertical slowness (vertical ray parameter) of incident waves along the horizontal source profile and the vertical receiver array. The slowness data combined with travel times are interpreted by forward modeling taking into account geological information of the survey area. Our findings confirm results from gravity measurements which suggest the existence of large amphibolite/metabasite complexes in the vicinity of the borehole. The described method is also used to identify P-to-S converted energy originating from fracture zones above the receiver array and to locate the region in which conversion occurs. 相似文献
16.
采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法. 相似文献
17.
在研究地壳结构的人工源宽角反射地震资料解释中,常规宽角反射波走时和射线路径计算大都假定地壳模型为层状块状均匀介质.为了逼近实际地壳结构模型,要求模型尺度较大,为了提高地震资料解释的可靠性,须减小模型离散单元的尺寸,但同时计算量大大增加,使资料解释的效率较低.为此,本文尝试同时提高宽角反射地震资料解释效率和可靠性的方法,即使用双重网格计算宽角反射地震波走时和射线路径的最小走时树方法.双重网格法在均匀介质内部仅计算大网格节点,在速度变化点、震源点和检波点区域,同时计算小网格节点;在界面边界点使用比介质内部节点更大的子波传播区域.模型计算结果表明,对于大尺度的层状块状均匀介质模型,在保证精度的条件下,本文所提出的双重网格射线追踪方法的计算效率比单网格方法显著提高. 相似文献
18.
RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA 下载免费PDF全文
Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection. 相似文献
19.
20.
To conduct forward and simultaneous inversion in a complex geological model, including an irregular topography (or irregular reflector or velocity anomaly), we in this paper combined our previous multiphase arrival tracking method (referred as triangular shortest-path method, TSPM) in triangular (2D) or tetrahedral (3D) cell model and a linearized inversion solver (referred to as damped minimum norms and constrained least squares problem solved using the conjugate gradient method, DMNCLS-CG) to formulate a simultaneous travel time inversion method for updating both velocity and reflector geometry by using multiphase arrival times. In the triangular/tetrahedral cells, we deduced the partial derivative of velocity variation with respective to the depth change of reflector. The numerical simulation results show that the computational accuracy can be tuned to a high precision in forward modeling and the irregular velocity anomaly and reflector geometry can be accurately captured in the simultaneous inversion, because the triangular/tetrahedral cell can be easily used to stitch the irregular topography or subsurface interface. 相似文献