Seismic profiles across the middle Tan-Lu fault zone in Laizhou Bay, Bohai Sea, eastern China

The stratigraphic architecture, structure and Cenozoic tectonic evolution of the Tan-Lu fault zone in Laizhou Bay, eastern China, are analyzed based on interpretations of 31 new 2D seismic lines across Laizhou Bay. Cenozoic strata in the study area are divided into two layers separated by a prominent and widespread unconformity. The upper sedimentary layer is made up of Neogene and Quaternary fluvial and marine sediments, while the lower layer consists of Paleogene lacustrine and fluvial facies. In terms of tectonics, the sediments beneath the unconformity can be divided into four main structural units: the west depression, central uplift, east depression and Ludong uplift. The two branches of the middle Tan-Lu fault zone differ in their geometry and offset: the east branch fault is a steeply dipping S-shaped strike-slip fault that cuts acoustic basement at depths greater than 8 km, whereas the west branch fault is a relatively shallow normal fault. The Tan-Lu fault zone is the key fault in the study area, having controlled its Cenozoic evolution. Based on balanced cross-sections constructed along transverse seismic line 99.8 and longitudinal seismic line 699.0, the Cenozoic evolution of the middle Tan-Lu fault zone is divided into three stages: Paleocene-Eocene transtension, Oligocene-Early Miocene transpression and Middle Miocene to present-day stable subsidence. The reasons for the contrasting tectonic features of the two branch faults and the timing of the change from transtension to transpression are discussed. Crown Copyright (C) 2008 Published by Elsevier Ltd. All rights reserved.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data

The stress release model, a stochastic version of the elastic-rebound theory, is applied to the historical earthquake data from three strong earthquake-prone regions of China, including North China, Southwest China, and the Taiwan seismic regions. The results show that the seismicity along a plate boundary (Taiwan) is more active than in intraplate regions (North and Southwest China). The degree of predictability or regularity of seismic events in these seismic regions, based on both the Akaike information criterion (AIC) and fitted sensitivity parameters, follows the order Taiwan, Southwest China, and North China, which is further identified by numerical simulations. (c) 2004 Elsevier Ltd. All rights reserved.

碳酸盐岩缝洞型储层地球物理响应特征及预测方法研究

Proven by the petroleum exploration activities, the karsts-fissure reservoir in carbonate rocks is significant to find out the large scale oil & gas field. They are made up of the four reservoir types: karsts-cave, karsts-crack, crack-cave and fracture-pore-cave. Each reservoir space and each reservoir bed has different features of reservoir heterogeneity and small scale of pore-crack-cave. The fracture-cave reservoir in carbonate rocks is characteristic by multi-types and long oiliness well. The reservoir shape is controlled by the irregular pore-crack-cave. The development level of fracture and karst-cave is the key element of hydrocarbon enriching, high productivity and stable production. However, most of Carbonate formation are buried deeply and the signal-ration-noise of seismic reflection are very low. It is reason why the fracture-cave reservoir are difficult to be predicted effectively. In terms of surveyed and studied lots of the former research outcome, The author applied the methods of synthetical reservoir geophysical prediction from two ways including macrosopic and microcomic technics in terms of the reservoir-cap condition, geophysics and geology feature and difficulty of prediction in carbonate rocks. It is guiden by the new ideas of stratigraphy, sedimentology, sedimentography, reservoir geology and karst geology. The geophysics technology is key technics. In aspects of macroscopic studies, starting off the three efficiencies of controlling the reservoir distribution including sedimental facies, karst and fracture, by means of comprehensive utilization of geology, geophysics, boring well and well log, the study of reservoir features and karst inside story are developed in terms of data of individual well and multiple well. Through establishing the carbonate deposition model, karstic model and fracture model, the macro-distribution laws of carbonatite are carried out by the study of coherence analysis, seismic reflection feature analysis and palaeotectonics analysis. In aspects of microcosmic studies, starting off analysis in reservoir geophysical response feature of fracture and karst-cave model according to guidance of the macroscopic geological model in carbonate reservoir, the methods of the carbonate reservoir prediction are developed by comprehensively utilization of seismic multi-attribution intersection analysis, seismic inversion restricted by log, seismic discontinuity analysis, seimic spectrum attenuation gradient, moniliform reflection feature analysis and multiparameter karst reservoir appraisement.Through application of carbonate reservoir synthetical geophysics prediction, the author r successfully develops the beneficial reservoir distribution province in Ordovician of Katake block 1in middle Tarim basin. The fracture-cave reservoir distributions are delineated. The prospect direction and favorable aims are demonstrated. There are a set of carbonate reservoir prediction methods in middle Tarim basin. It is the favorable basic technique in predicting reservoir of the Ordovician carbonate in middle Tarim. Proven by exploration drilling, the favorable region of moniliform reflection fracture and pore-cave and cave-fracture in lower-middle Ordovician are coincidence with the region of hydrocarbon show. It’s indicated that the reservoir prediction methods described in the study of Ordovician carbonate formation are feasible practicably.

全信息精细储层预测方法研究与实践-以周清庄油田沙三段油藏为例

With the continuously proceeding of petroleum exploratory development in China, exploratory development becomes more and more difficult. For increasing reserve volume and production, lithologic hydrocarbon reservoir has been the most workable, potential and universality exploration targets. In the past, Dagang Oil Field use the complicated fault reservoir theory as the guide, develop and form a suit of matching construction and instrument in prospecting complicated fault reservoir that reach top of exploration industry in China. But the research of lithologic hydrocarbon reservoir is not much, which affects the exploitation progress of lithologic hydrocarbon reservoir. In this thesis, is object, through the depth study of lithologic deposition in Shasan segment of Zhouqingzhuang Oil Field, a suit of holographic fine reservoir bed forecasting techniques is built up and finally gets following main results: 1. Applying geology, seism, drilling, logging and other information to sensitivity preferences, geological model, inversion and integrated stratum evaluation, realizing the method and flow of refined multi-information stratum forecast. 2. Built up a full three dimensional fine structural interpretation method: in view of r problem of accurately demarcating 90% inclined well, propose a inclined well air space demarcating method, make bed demarcating more exactly; in view of problem of faults demarcating and combination in seismic interpretation, propose a computational method of seismic interference based on wavelet translation, make identify the fault in different level more dependable and reasonable; for exactly identifying structural attitude, propose a velocity modeling method under multi-well restriction, make structural attitude closer to the facts. 3. Built up a high accuracy reservoir bed inversion method: in view of problem in exactly identifying reservoir and nonreservoir with conventional wave impedance inversion method in this place, propose a reservoir log response characteristic analysis and sensible log parameter inversion method. ①analysis log response of reservoir and nonreservoir in region of interest, make definite the most sensible log parameter in identifying reservoir and nonreservoir in this region; ②make sensible log parameter inversion based on wave impedance inversion, to improve inversion accuracy, the thickness of recognizable reservoir bed reach 4-5m. 4. Built up a 4-D reservoir forcasting circuit: in view of difficulty that in lithologic hydrocarbon reservoir making reservoir space characteristic clear by using structural map and reservoir forecasting techniques once only, propose a 4-D reservoir forcasting circuit. In other words, based on development conceptual design, forcast reservoir of different time, namely multiple 3D reservoir forcasting in time queue, each time the accuracy degree of reservoir forcasting is improved since apply the new well material, thereby achieve high quality and highly efficient in exploratory development. During exploratory development lithologic depositin in Shasan segment of Zhouqingzhuang Oil Field, there are thirteen wells get 100% success rate, which sufficiently proves that this suit of method is scientific and effective.

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

天然气藏地震识别技术研究

The Research of Seismic Recognition Techniques for Gas Reservoir Shang Yong_sheng(Geophysics) Directed by Yang Chang-chun Abstract Gas reservior is one of the most important nature resources. Someone forecast that the output will exceed crude oil in 2015 and become the largest energy source. Recently,more and more gas reservior are discovered as the oil field and gas filed exploration go deep into development. Although the gas proved reserves rise greatly the explorative degree of natural gas resource in our country is still very low. The potential of gas exploration is very great and our task is so hard. How to recognise and discover new gas reservoir is the first task based on the great gas reservior resources foreground. the gas reservior in different oil and gas field have its special gas generation, reservoiring, physical property conditions. However,it may have the same geophysical characters. So,it is very important to analyse, research, summarizing the geophysical characters of the gas reservior and make use of the characters to identify the gas layer effectively. This paper start with modeling,and it discuss the geophysical characters of the gas reservior response. It analyse the seismical wave characters of the gas reservoir. Furthermore, it summarize the method of using the seismica profile to identify the gas reservior directly. The paper discuss the research of extracting diffraction wave for mass diffraction wave grow at the edge of the gas reservoir at the seismic section. Making use of the technique of extracting diffraction wave to identify the gas reservior is the first experiment of the gas reservoir prediction technique. The avo technology is a new geophysical method. From the pre-stack analysis, this paper discuss the technique of using the rich information to identify the gas reservoir. Based on the case study of the Qidam basin and the Hailaer basin it discuss the method of predicating gas reservoir using pre-stack information. It include pre-stack amplitude preserve process, AVO modeling, fluid replacement technique, AVO analysis and interpretation technique. The paper summarize a gas reservoir prediction procedure focusing on the pre-stack information. The seismic wave will cause great attenuation when it pass through the gas layer and the high frequency component loss. This paper discuss the technique of extracting seismic attributes to represent the attenuation degree of seismica wave. Based on the attenuation attributes,it does the research of the gas reservor identification and prediction. At last,the paper discuss the method of calculating the azimuthal anisotropy to predict the fracture reservoir. Keyword: gas reservoir, diffraction wave, AVO, attenuation attribute,fracture prediction

博兴洼陷新生代构造-沉积演化及其对鲁西隆起的响应

The West Shandong Uplift and its adjacent basins, with same evolutional history before Mesozoic, are an important basin-orogenic systems in North China. After late Mesozoic, tectonic differentiation between basin and orogenic belt gradually displayed in the study area. The Boxing sag is a part of Jiyang Depression near to West Shandong Uplift, in which the whole Mesozoic and Cenozoic strata are preserved. Based on the analysis of sedimentary records in the Boxing sag, the Cenozoic structural and sedimentary evolutions in Boxing Sag and its response to Western Shandong uplift are discussed in this dissertation. The main conclusions in this research are presented as follows. Based on Seismic and well logging profile interpretation, fault growth index, thickness difference between bottom wall and top wall and fault activity rate from Eocene to Pliocene are studied. Boxing sag had three main faults, NE, NW and NEE trending faults. Research shows that the activity of the NW trending fault in the Boxing sag became weaken from E1-2S4 to N2m gradually. The evolution of NE and the NEE trending fault can be divided into three episodes, from E1-2k to E2s4, from E2s3 to E3s1, from N2m to E3d. The analysis of Paleogene samples of heavy mineral assemblages shows that metamorphic rocks represented by garnet, intermediate-acid igneous rocks represented by the assemblage of apatite, zircon and tourmaline became less from E1-2k to N2g, and sedimentary rocks represented by the assemblage of pyrite, barite and limonite also became less. Intermediate-basic igneous rocks represented by the assemblage of leucoxene, rutile and ilmenite and metamorphic rocks represented by epidote became more and more. Electronic microprobe analysis shows that glaucophane and barroisite are existed in Kongdian Formation and the 4th member of Shahejie Formation, and they demonstrate that Western Shandong and Eastern Shandong are all the source regions of the Boxing Sag, and they also indicate that oceanic crust existed before the collision between the Yangtze and North China continent. The fact that Eastern Shandong is the source region of Boxing Sag also indicates that Western Shandong was not high enough to prevent sediment from Eastern Shandong at E1-2k and E2s4. The results of the dating of five detrital zircons of Boxing Sag show Kongdian Formation and the 4th member of Shahejie Formation have the age peaks of 2800Ma and 700-800. It means that Eastern Shandong is the source region of Boxing Sag at early Paleogene and Western Shandong is not high enough to prevent the sediment from Eastern Shandong. The ages of 160-180 and 220-260 Ma, which exist in the Guantao Formation and Paleogene, are common in Eastern Shandong and rare in Western Shandong,and it implied that Western Shandong is a low uplift at 24Ma. The Paleogene strata have almost same age groups, while the Guantao Formation has significant variations of age groups, and this indicates that Boxing Sag and Western Shandong uplift had taken place tremendous changes. The results of apatite fission track in Boxing sag show that three times uplifts happened at the source region at 60 Ma, 45Ma and 15Ma respectively, and the Boxing sag experienced two subsidences at 60Ma, 45Ma and one uplift at 20Ma.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

地震P-S波时差耦合作用的斜坡崩滑效应研究

China locates between the circum-Pacific and the Mediterranean-Himalayan seismic belt. The seismic activities in our country are very frequent and so are the collapses and slides of slope triggered by earthquakes. Many collapses and slides of slope take place mainly in the west of China with many earthquakes and mountains, especially in Sichuan and Yunnan Provinces. When a strong earthquake happening, the damage especially in mountains area caused by geological hazards it triggered such as rock collapses, landslides and debris flows is heavier than that it caused directly. A conclusion which the number of lives lost caused by geological hazards triggered by a strong earthquake in mountains area often accounts for a half even more of the total one induced by the strong earthquake can be made by consulting the statistical loss of several representative earthquakes. As a result, geological hazards such as collapses and slides of slope triggered by strong earthquakes attract wide attention for their great costs. Based on field geological investigation, engineering geological exploration and material data analysis, chief conclusions have been drawn after systematic research on formation mechanism, key inducing factors, dynamic characteristics of geological hazards such as collapses and slides of slope triggered by strong earthquakes by means of engineering geomechanics comprehensive analysis, finite difference numerical simulation test, in-lab dynamic triaxial shear test of rock, discrete element numerical simulation. Based on research on a great number of collapses and landslides triggered by Wenchuan and Xiaonanhai Earthquake, two-set methods, i.e. the method for original topography recovering based on factors such as lithology and elevation comparing and the method for reconstructing collapsing and sliding process of slope based on characteristics of seism tectonic zone, structural fissure, diameter spatial distribution of slope debris mass, propagation direction and mechanical property of seismic wave, have been gotten. What is more, types, formation mechanism and dynamic characteristics of collapses and slides of slope induced by strong earthquakes are discussed comprehensively. Firstly, collapsed and slided accumulative mass is in a state of heavily even more broken. Secondly, dynamic process of slope collapsing and sliding consists of almost four stages, i.e. broken, thrown, crushed and river blocked. Thirdly, classified according to failure forms, there are usually four types which are made up of collapsing, land sliding, land sliding-debris flowing and vibrating liquefaction. Finally, as for key inducing factors in slope collapsing and sliding, they often include characteristics of seism tectonic belts, structure and construction of rock mass, terrain and physiognomy, weathering degree of rock mass and mechanical functions of seismic waves. Based on microscopic study on initial fracturing of slope caused by seismic effect, combined with two change trends which include ratio of vertical vs. horizontal peak ground acceleration corresponding to epicentral distance and enlarging effect of peak ground acceleration along slope, key inducing factor of initial slope fracturing in various area with different epicentral distance is obtained. In near-field area, i.e. epicentral distance being less than 30 km, tensile strength of rock mass is a key intrinsic factor inducing initial fracturing of slope undergoing seismic effect whereas shear strength of rock mass is the one when epicentral distance is more than 30 km. In the latter circumstance, research by means of finite difference numerical simulation test and in-lab dynamic triaxial shear test of rock shows that initial fracture begins always in the place of slope shoulder. The fact that fracture strain and shear strength which are proportional to buried depth of rock mass in the place of slope shoulder are less than other place and peak ground acceleration is enlarged in the place causes prior failure at slope shoulder. Key extrinsic factors inducing dynamic fracture of slope at different distances to epicenter have been obtained through discrete element numerical simulation on the total process of collapsing and sliding of slope triggered by Wenchuan Earthquake. Research shows that combined action of P and S seismic waves is the key factor inducing collapsing and sliding of slope at a distance less than 64 km to initial epicenter along earthquake-triggering structure. What is more, vertical tensile action of P seismic wave plays a leading role near epicenter, whereas vertical shear action of S seismic wave plays a leading role gradually with epicentral distance increasing in this range. On the other hand, single action of P seismic wave becomes the key factor inducing collapsing and sliding of slope at a distance between 64 km and 216 km to initial epicenter. Horizontal tensile action of P seismic wave becomes the key factor gradually from combined action between vertical and horizontal tensile action of P seismic wave with epicentral distance increasing in this distance range. In addition, initial failure triggered by strong earthquakes begins almost in the place of slope shoulder. However, initial failure beginning from toe of slope relates probably with gradient and rock occurrence. Finally, starting time of initial failure in slope increases usually with epicentral distance. It is perhaps that the starting time increasing is a result of attenuating of seismic wave from epicenter along earthquake-triggering structure. It is of great theoretical and practical significance for us to construct towns and infrastructure in fragile geological environment along seism tectonic belts and conduct risk management on earthquake-triggered geological hazards by referring to above conclusions.

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

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.