地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

湖南、广西泥盆纪硅岩地球化学与沉积构造背景

During the Devonian, a complicated carbonate platform-basin configuration was created through transtensional rifting in the context of opening of Devonian South China Sea; extensive bedded chert, commonly interbedded with tuffaceous beds, occurred in the narrow, elongate interplatform basins (or troughs) in South China, where they occurred earlier (Early Devonian) in southern Guangxi and later (early Late Devonian) in northern Guangxi-south central Hunan. In order to unravel the origin and distribution of the bedded chert successions, and their relationships to basement faulting activities during the opening of the Devonian South China Sea, studies of element (major, minor and REE) geochemistry and Rb-Sr, Sm-Nd isotopic systematics are carried out upon the chert deposits. These chert deposits commonly have high SiO2 contents and (average 94.01%) and low TFe2O3 (average 0.55%), together with other geochemical parameters, suggestive of both biogenic and hydrothermal origins. However, Fe/Ti ratio are high along the elongate interplatform basins(troughs) to the northwest along Wuxiangling-Zhaisha-Chengbu, and to the southeast along Xiaodong-Mugui-Xinpu, suggesting relatively intense hydrothermal activities there. They generally contain very low total REE contents (∑REE average 31.21ug/g) with mediate negative Ce anomalies (mean Ce/Ce*=0.83) and low Lan/Cen values (average 1.64), indicating an overall continental margin basin where they precipitated. The northward increases in Ce/Ce* values, particularly along the elongate troughs bounded both to the east and west of the Guangxi-Huanan rift basin, suggest a northward enhancement of terrigenous influences, thereby reflecting a gradual northward propagation of open marine setting. Generally low positive Eu anomalies in the chert, except for the apparently high Eu anomalies in the chert from Chengbu (Eu/Eu* up to 4.6), suggest mild hydrothermal venting activities in general, except for those at Chengbu. The initial 87Sr/86Sr (0) ratios of chert generally vary from 0.712000 to 0.73000 , suggesting influences both from terrigenous influx and seawater. The Nd isotopic model ages (tDM or t2DM) and initial εNd (0) values of chert vary mostly from 1.5 to 2.1 Ga, and from –16 to –21, respectively, implying that the silica sources were derived from the provenances of the Palaeoproterozoic crust relics at depth. The high εNd (0) values of chert (-0.22 to 14.7) in some localities, mostly along the elongate troughs, suggest that silica sources may have been derived from deeper-seated mantle, being channeled through the interplate boundary fault zones extending downwards to the mantle. At Wuxiangling, Nanning, chert occurs extensively from the Emsian through the Frasnian strata, both U/Th ratios and tDM ages of chert reached up to a maximum in the early Frasnian corresponding to the extensive development of chert in South China, pointing to a maximum extensional stage of Devonian South China basin, which is supported by the Ce/Ce* values as is opposed to the previous datasets as the coeval minimum values.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

纵波方位属性裂缝检测技术综合研究

Fractured oil and gas reservoir is an important type of oil and gas reservoir, which is taking a growing part of current oil and gas production in the whole world. Thus these technologies targeted at exploration of fractured oil and gas reservoirs are drawing vast attentions. It is difficult to accurately predict the fracture development orientation and intensity in oil and gas exploration. Focused on this problem, this paper systematically conducted series study of seismic data processing and P-wave attributes fracture detection based on the structure of ZX buried mountain, and obtained good results. This paper firstly stimulated the propagation of P-wave in weak anisotropic media caused by vertical aligned cracks, and analyzed the rule of P-wave attributes’ variation associated with observed azimuth, such as travel-time, amplitude and AVO gradient and so on, and quantitatively described the sensitive degree of these attributes to anisotropy of fracture medium. In order to further study the sensitive degree of these attributes to anisotropy of fractures, meanwhile, this paper stimulated P-wave propagation through different types and different intensity anisotropic medium respectively and summarized the rule of these attributes’ variation associated with observed azimuth in different anisotropic medium. The results of these studies provided reliable references for predicting orientation, extensity and size of actual complicated cracked medium by P-wave azimuth attributes responses. In the paper, amounts of seismic data processing methods are used to keep and recover all kinds of attributes applied for fracture detection, which guarantee the high accurate of these attributes, thus then improve the accurate of fracture detection. During seismic data processing, the paper adopted the three dimensional F-Kx-Ky field cone filter technique to attenuate ground roll waves and multiple waves, then enhances the S/N ratio of pre-stack seismic data; comprehensively applying geometrical spread compensation, surface consistent amplitude compensation, residual amplitude compensation to recover amplitude; common azimuth processing method effectively preserves the azimuthal characteristics of P-wave attributes; the technique of bend ray adaptive aperture pre-stack time migration insures to obtain the best image in each azimuth. Application of these processing methods guaranteed these attributes’ accuracy, and then improved the accuracy of fracture detection. After comparing and analyzing a variety of attributes, relative wave impedance (relative amplitude) attribute is selected to inverse the orientation of fracture medium; attenuation gradient and corresponding frequency of 85% energy are selected to inverse the intensity of fracture medium; then obtained the fracture distribution characteristics of lower Paleozoic and Precambrian in ZX ancient buried mountains. The results are good accord with the characteristics of faults system and well information in this area.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.

水-岩土反应的力学与环境效应研究

The research on mechanical effects of water-rock and soil interaction on deformation and failure of rocks and soils involves three aspects of mechanics, physics and chemistry. It is the cross between geochemistry and rock mechanics and soil mechanics. To sum up, the mechanical effects of water-rock and soil interaction is related to many complex processes. Research in this respect has been being an important forward field and has broad prospects. In connection with the mechanism of the effects of the chemical action of water-rock on deformation and failure of rocks and soils, the research significance, the present state, the developments in this research domain are summarized. Author prospects the future of this research. The research of the subject should be possessed of important position in studying engineering geology and will lead directly to a new understand on geological hazard and control research. In order to investigation the macroscopic mechanics effects of chemical kinetics of water-rock interaction on the deformation and failure, calcic rock, red sandstone and grey granite reacting chemically with different aqueous solution at atmospheric temperature and atmospheric pressure are uniaxially compressed. The quantitative results concerning the changes of uniaxially compressive strength and elastic modulus under different conditions are obtained. It is found that the mechanical effects of water on rock is closely related to the chemical action of water-rock or the chemical damage in rock, and the intensity of chemical damage is direct ratio to the intensity of chemical action in water-rock system. It is also found that the hydrochemical action on rock is time-dependent through the test. The mechanism of permeation and hydrochemical action resulting in failure of loaded rock mass or propagation of fractures in rocks would be a key question in rock fracture mechanics. In this paper, the fracture mechanical effects of chemical action of water-rock and their time- and chemical environment-dependent behavior in grey granite, green granite, grey sandstone and red sandstone are analyzed by testing K_(IC) and COD of rock under different conditions. It is found that: ①the fracture mechanical effect of chemical action of water-rock is outstanding and time-dependent, and high differences exist in the influence of different aqueous solution, different rocks, different immersion ways and different velocity of cycle flow on the fracture mechanical effects in rock. ②the mechanical effects of water-rock interaction on propagation of fractures is consistent with the mechanical effects on the peak strength of rock. ③the intensity of the mechanical fracture effects increases as the intensity of chemical action of water-rock increases. ④iron and calcium ion bearing mineral or cement in rock are some key ion or chemical composition, and especially iron ion-bearing mineral resulting in chemical action of water-rock to be provided with both positive and negative mechanical effects on rock. Through the above two tests, we suggest that primary factors influencing chemical damage in rock consist of the chemical property of rock and aqueous solution, the structure or homogeneity of rocks, the flow velocity of aqueous solution passing through rock, and cause of formation or evolution of rock. The paper explores the mechanism on the mechanical effects of water-rock interaction on rock by using the theory of chemistry and rock fracture mechanics with chemical damage proposed by author, the modeling method and the energy point of view. In this paper, the concept of absorbed suction between soil grains caused by capillary response is given and expounded, and the relation and basic distinction among this absorbed suction, surface tension and capillary pressure of the soil are analyzed and established. The law of absorbed suction change and the primary factors affecting it are approached. We hold that the structure suction are changeable along with the change of the saturation state in unsaturated soils. In view of this, the concept of intrinsic structure suction and variable structure suction are given and expounded, and this paper points out: What we should study is variable structure suction when studying the effective stress. By IIIy κHH's theory of structure strength of soils, the computer method for variable structure suction is analyzed, the measure method for variable structure suction is discussed, and it reach the conclusions: ①Besides saturation state, variable structure suction is affected by grain composition and packing patter of grains. ②The internal relations are present between structure parameter N in computing structure suction and structure parameter D in computing absorbed suction. We think that some problems exit in available principle of effective stress and shear strength theory for unsaturated soil. Based on the variable structure suction and absorbed suction, the classification of saturation in soil and a principle of narrow sense effective stress are proposed for unsaturated soils. Based on generalized suction, the generalized effective stress formula and a principle of generalized effective stress are proposed for unsaturated soils. The experience parameter χ in Bishop's effective stress formula is defined, and the principal factors influencing effective stress or χ. The primary factor affecting the effective stress in unsaturated soils, and the principle classifying unsaturated soils and its mechanics methods analyzing unsaturated soils are discussed, and this paper points out: The theory on studying unsaturated soil mechanics should adopt the micromechanics method, then raise it to macromechanics and to applying. Researching the mechanical effects of chemical action of water-soil on soil is of great importance to geoenvironmental hazard control. The texture of soil and the fabric of soil mass are set forth. The tests on physical and mechanical property are performed to investigate the mechanism of the positive and negative mechanical effects of different chemical property of aqueous solution. The test results make clear that the plastic limit, liquid limit and plasticity index are changed, and there exists both positive and negative effects on specimens in this test. Based on analyzing the mechanism of the mechanical effects of water-soil interaction on soil, author thinks that hydrochemical actions being provided with mechanical effects on soil comprise three kinds of dissolution, sedimentation or crystallization. The significance of these tests lie in which it is recognized for us that we may improve, adjust and control the quality of soils, and may achieve the goal geological hazard control and prevention.The present and the significance of the research on environmental effects of water-rock and soil interaction. Various living example on geoenvironmental hazard in this field are enumerated. Following above thinking, we have approached such ideals that: ①changing the intensity and distribution of source and sink in groundwater flow system can be used to control the water-rock and soil interaction. ②the chemical action of water-rock and soil can be used to ameliorate the physical and mechanical property of rocks and soils. Lastly, the research thinking and the research methods on mechanical effects and environmental effects of water-rock and soil interaction are put forward and detailed.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

中-新元古代海洋的硫和碳同位素的演化

Data on seawater carbon isotope in the Mesoproterozoic and Neoproterozoic is abundant. However, the sulfur isotopic age curve of seawater sulfates determined through the analysis of sulfur isotopic composition of marine evaporite is uncertain in the Mesoproterozoic and Neoproterozoic since evaporites are generally rare in Precambrian. The Mesoproterozoic and Neoproterozoic Carbonate Formations preserve not only the carbon isotopic records, but also the sulfur isotopic records of coeval seawater in the Huabei Platform and the Yangtze Platform, China. Sulfur isotopic composition can be determined by the extraction of trace sulfate from carbonate samples. Successive measurements of sulfur and carbon isotopic compositions of carbonate samples from the Mesoproterozoic and Neoproterozoic strata in the Huabei Platform and the Yangtze Platform was accomplished through the extracting of trace sulfate from carbonates. Sulfur and carbon isotopic compositions of coeval seawater were obtained from analytical results of sulfur and carbon isotopes of the same sample without diagenetic alteration. The high-resolution age curve of sulfur isotope given in this paper may reflect the trend of variations in sulfur isotope composition of seawater sulfates during the Mesoproterozoic and Neoproterozoic. It can be correlated with the characteristics of variation in age curve of carbon isotope of coeval seawater carbonates. The δ34S values of seawater varied from +10.3-37.0‰ during the Mesoproterozoic, which took on oscillated variation on the whole. The δ34S values took on high values in the Mesoproterozoic Chuanlinggou stage, Tuanshanzi stage Tieling stage and in Neoproterozoic Jing'eryu stage. The average of those was about +30‰. The sulfates have low δ34S values in the Mesoproterozoic Yangzhuang stage and Hongshuizhuang stage, The average of those was all lower than +20‰. There occured large-amplitude changs in δ34S values of seawater during the Mesoproterozoic. Large-amplitude oscillate of 534S values occured in the intervals of 1600~1400Ma and 1300~1200Ma. The δ13C values of seawater are mostly negative in Changcheng stage of late Paleoproterozoic, -0 ± 1‰ range in Jixian stage of Mesoproterozoic , and the positive 2±2‰ commonly in early Neoproterozoic Jing'eryu stage. From 1000 Ma to 900 Ma, about 108 years interval of oceanic 513C record is shortage. At the end of Paleoproterozoic (1700 - 1600 Ma), the oceanic 813C values change from -3‰ to 0‰, but strongly oscillate near 1600 Ma. Two larger variations of seawater 513C values occur in the Mesoproterozoic: one is a cycle of about 4%o happens at ca. 1400 Ma; another is rise from >2‰ to>5‰ at ca. 1250 Ma and then become stable at the near 1000 Ma. There appears a large positive excursion over +20‰ in 534S value of ancient seawater sulfates in the early Doushantuo stage. Simultaneously, 8 C values of ancient seawater occur a positive excursion reaching 10‰. These allow δ4S values and 513C values to reach high values of+51.7‰ and +6.9‰, respectively. The range of variation in 834S values of seawater is relatively narrow and 513C values are quite high in the middle Doushantuo stage. Then, δ34S values of seawater become oscillating, the same happens in δ13C values. Negative excursions in 834S values and 813C values occur simultaneously at the end of the Doushantuo stage, and the minimum of δ34S values and δ13C values dropped to -11.3‰ and -5.7‰, respectively. The ancient seawater in the Dengying stage has high δS values and δ13C values. Most of the δ34S values of the trace sulfate samples varied between +23.6‰ and +37.9‰ except two boundaries of the Dengying Formation, and the S13C values of the carbonate samples of the Dengying Formation varied between +0.5‰ and +5.0‰. There appeared large negative excursion in 834S values and δ13C values of ancient seawater at the bounder of Precambrian-Cambrian. The isotopic characteristics of sulfur and carbon implicated that the organic productivity and isotopic fractionation caused by biology were low and the palaeoceanic environment was quite unstable during the Mesoproterozoic. The increase and subsequent oscillation of seawater δ13C value occurred from 1700 to 1600 Ma and near 1300 Ma may be responsible to the two global tectonic events happened at coeval time. The characteristics of variation in sulfur and carbon isotopes of ancient seawater imply strong changes in oceanic environment, which became beneficial to inhabitation and propagation of organism. The organic production and the burial rate of organic carbon once reached a quite high level during the Doushantuo stage. However, the state of environment became unstable that means the global climate and the environment possibly were fluctuating and reiterating after the global glaciation. The negative excursions of S34S values and δ13C values occurring at the end of the Doushantuo stage represent a global event, which might be relative to the oxidation of deep seawater. The isotopic characteristics of sulfur and carbon implicated that there were a high organic productivity and a high burial rate of organic carbon in the Dengying stage. It is obvious that the palaeoceanic environment in Dengying stage was stable corresponding and beneficial for biology to inhabit and propagate except for the two boundaries. The tendency of sulfur and carbon isotopic variations maybe resulted from the gradual oxygenation of ocean environment during the Dengying stage. It has been reported that the secular variations of the sulfur isotopic compositions in seawater was negative correlated with that of carbon isotopic compositions. However, our results show that it is not the case. They were negatively correlated in some intervals and positively in some other intervals of the Mesoproterozoic and Neoproterozoic. The difference in correlation may be associated with the changes in conditions of redox in oceanic environment, e.g. sharp change of the oxidation-reduction interface. The strong changes in global environment may induce the abnormality to occur in the biogeo chemical S and C cycles in the ocean and accordingly sharp Variations in isotopic composition of seawater sulfur and carbon during the Mesoproterozoic and Neoproterozoic. Simultaneously, the global tectonism caused large changes of 87Sr/86Sr ratios. The leading factor that causes the variation in isotopic composition is different in the different intervals of the Mesoproterozoic and Neoproterozoic. Thus, there may exist different models of the biogeochemical S and C cycles in the ocean during the Mesoproterozoic and Neoproterozoic.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

地震偏移宽角吸收边界条件和Toeplitz矩阵反演中的边界效应研究

The content of this paper is based on the research work while the author took part in the key project of NSFC and the key project of Knowledge Innovation of CAS. The whole paper is expanded by introduction of the inevitable boundary problem during seismic migration and inversion. Boundary problem is a popular issue in seismic data processing. At the presence of artificial boundary, reflected wave which does not exist in reality comes to presence when the incident seismic wave arrives at the artificial boundary. That will interfere the propagation of seismic wave and cause alias information on the processed profile. Furthermore, the quality of the whole seismic profile will decrease and the subsequent work will fail.This paper has also made a review on the development of seismic migration, expatiated temporary seismic migration status and predicted the possible break through. Aiming at the absorbing boundary problem in migration, we have deduced the wide angle absorbing boundary condition and made a compare with the boundary effect of Toepiitz matrix fast approximate computation.During the process of fast approximate inversion computation of Toepiitz system, we have introduced the pre-conditioned conjugate gradient method employing co circulant extension to construct pre-conditioned matrix. Especially, employment of combined preconditioner will reduce the boundary effect during computation.Comparing the boundary problem in seismic migration with that in Toepiitz matrix inversion we find that the change of boundary condition will lead to the change of coefficient matrix eigenvalues and the change of coefficient matrix eigenvalues will cause boundary effect. In this paper, the author has made an qualitative analysis of the relationship between the coefficient matrix eigenvalues and the boundary effect. Quantitative analysis is worthy of further research.