离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

均质弹性条件下多洞平面问题解析算法及其应用研究

The scholars in the world have been trying to find an effective analytic algorithm of multiple hole problems usually meet in engineering designs. Though some studies on circular or elliptic holes had been achieved under specific conditions, no efforts were made to any multiple hole problems that is most significant for engineering designs. The author has made further studies on any multiple hole problems, using complex variable function method and Schwarz alternating method. After solving a series of technological difficulties, the author obtains an effective analytic algorithm, and acquires stress field and displacement field with high accuracy, which can be conducted for arbitrary many iterations according to practical accuracy requirements. In addition, th solution of stress and displacement fields, even for multiple holes of complex shapes and smaller distances. Further, the author made preliminary studies on viscoelastic displacement solution for any double holes. In terms of the obtained displacement solution of any multiple holes, this paper studies displacement back-analysis for the excavations of two tunnels, and find that the back-analysis method is accurate. Additionally, the author presents the mathematical prove of inversion uniqueness for ground stresses, elastic modulus and Poisson ratio. The author believes that the accurate analytic algorithm provided in this paper will presents an effective way to stress and displacement analysis for any multiple hole problems, optimal arrangement of multiple holes, hole shape optimization of multiple holes, etc..