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

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

中国东部上第三系硬粘土的工程地质特性及边坡开挖工程环境效应研究

The engineering geological properties of Neogene hard clays and related engineering problems are frontiers in the fields of Engineering Geology, Soil Mechanics and Rock Mechanics. Recently, it has been recognized that Neogene hard clay is the intermediate type of material between the soil and the rock. Many aspects of them, such as sampling, testing, calculating and engineering process, are special, which could not be researched by means of theories and methods of traditional Soil Mechanics of Rock Mechanics. In order to get real knowledge and instruct the engineering practice, intersect studying of multiple disciplines, including Engineering Geology, Soil Mechanics and Rock Mechanics, etc., is necessary. Neogene hard clay is one of the important study objects of regional problem rocks & soils in our country, which extensively distributed in China, especially in Eastern China. Taking the related areas along the middle line of the Project of Transferring Water from the South to the North (e.g. Nanyang basin, Fangcheng-Baofeng area and Handan-Yongnian area), South-west of Shandong, Xu-Huai area and Beijing area, etc. as main study areas, the paper divided Neogene hard clays into reduction environment dominated origin and oxidation environment dominated origin, which distributed on areas western and eastern to Mount Taihangshan respectively. Intermediate types are also existed in some areas, which mainly distribute near the edges of depositional basins; they are usually of transitions between diluvial and lacustrine deposits. As to Neogene hard clays from Eastern China, the clay particle content is high, and montmorillonite or illite/montmorillonite turbostratic mineral is the dominating clay mineral. The content of effective montmorillonite is very high in each area, which is the basis for the undesirable engineering properties of Neogene hard clays. For hard clays from the same area, the content of effective montmorillonite in gray-greenish hard clay is much higher than that in purple-brownish or brown-yellowish hard clay, which is the reason why the gray-greenish hard clay usually has outstanding expansive property. On the other hand, purple-brownish or brown-yellowish hard clay has relatively less montmorillonite, so its property is better. All of these prove that the composition (clay mineral) of Neogene hard clay is the control factor for the engineering properties. Neogene hard clays have obvious properties such as fissured, overconsolidated and expansive, which are the main reasons that many engineering problems and geological harzards usually occur in Neogene hard clays. The paper systematically elaborates the engineering properties of Neogene hard clays from Eastern China, analyses the relationships between engineering properties and basic indexes. The author introduces the ANN method into the prediction of engineering property indexes of hard clays, which provides a new way for quantitatively assessment and prediction of engineering property indexes. During investigation in the field, the author found that there exists obvious seam-sheared zone between different hard clays in Miocene Xiacaowan formation in Xu-Huai area. Similar phenomenon also exists near the borderline between Neogene hard clays and underlying coal measures in the Southwest of Shandong province, which could be observed in the cores. The discovery of seam-sheard zone has important theoretical and practical significance for engineering stability analysis and revealing the origin of fissures in Neogene hard clays. The macrostructure, medium structure and microstructure together control the engineering properties of hard clays. The author analyses and summarizes the structural effects on hard clays in detail. The complex of the strength property of hard clays is mostly related to the characteristics of fissures, which is one of the main factors that affect the choice of shear strength parameters. So structure-control theory must be inseparably combined with composition-control theory during the engineering geological and rock/soil mechanics research of hard clays. The engineering properties, such as fissured, overconsolidated and expansive, control the instability of engineering behaviors of Neogene hard clays under the condition of excavation, i.e. very sensitive to the change of existence environment. Based on test data analysis, the author elaborates the effects of engineering environment change on the engineering properties. Taking Nanyang basin as example, the author utilizes FEM to study the effects of various factors on stability of cutting canal slopes, than sets forth the characteristics, development laws and formation mechanism of the deformation and failure of hard clay canal slopes, summarizes the protection and reinforcement principles, as well as the protection and remedy steps. On the basis of comparison of engineering properties of domestic and foreign Neogene muddy deposits, in the view of whole globe and associated with the geological characteristics of China, the paper demonstrates that the intermediate type of the material between the soil and the rock, named "hard clay/soft rock", which can not be separated abruptly, really exists in China. The author has given a preliminary classification based on its geological origin and distribution law, which is very significant for promoting the mixture of Engineering Geology, Soil Mechanics and Rock Mechanics. In the course of large scales engineering construction in China, many engineering experiences and testing data are gained, summarizing these testing results and automatically managing them with computer technology are very necessary. The author develops a software named "Hard Clay-Soft Rock Engineering Geological Information Management and Analysis System (HRGIMS)", realizes the automatic and visual management of geo-engineering information, on the basis of information management, the functions of test data analysis and engineering property prediction are strengthened. This system has well merits for practice and popularization.

各向同性与各向异性介质中地震波旅行时正演模拟与反演成像

By seismic tomography, interesting results have been achieved not only in the research of the geosphere with a large scale but also in the exploration of resources and projects with a small scale since 80'. Compared with traditional inversion methods, seismic tomography can offer more and detailed information about subsurface and has been being paid attention by more and more geophysicists. Since inversion based on forward modeling, we have studied and improved the methods to calculate seismic traveltimes and raypaths in isotropic and anisotropic media, and applied the improved forward methods to traveltime tomography. There are three main kinds of methods to calculate seismic traveltime field and its ray path distribution, which are ray-tracing theory, eikonal equation by the finite-difference and minimum traveltime tree algorithm. In ray tracing, five methods are introduced in the paper, including analytic ray tracing, ray shooting, ray bending, grid ray tracing and rectangle grid ray perturbation with three points. Finite-difference solution of eikonal equation is very efficient in calculation of seismic first-break, but is awkward in calculation of reflection traveltimes. We have put forward a idea to calculate traveltimes of reflected waves using a combining way of eikonal equation method and other one in order to improve its capability of dealing with reflection waves. The minimum traveltime tree algorithm has been studied with emphases. Three improved algorithms are put forward on the basis of basic algorithm of the minimum traveltime tree. The first improved algorithm is called raypath tracing backward minimum traveltime algorithm, in which not only wavelets from the current source but also wavelets from upper source points are all calculated. The algorithm can obviously improve the speed of calculating traveltimes and raypaths in layered or blocked homogeneous media and keep good accuracy. The second improved algorithm is raypath key point minimum traveltime algorithm in which traveltimes and raypaths are calculated with a view of key points of raypaths (key points of raypths mean the pivotal points which determine raypaths). The raypath key point method is developed on the basis of the first improved algorithm, and has better applicability. For example, it is very efficient even for inhomogeneous media. Another improved algorithm, double grid minimum traveltime tree algorithm, bases upon raypath key point scheme, in which a model is divided with two kinds of grids so that the unnecessary calculation can be left out. Violent undulation of curved interface often results in the phenomenon that there are no reflection points on some parts of interfaces where there should be. One efficacious scheme that curved interfaces are divided into segments, and these segments are treated respectively is presented to solve the problem. In addition, the approximation to interfaces with discrete grids leads to large errors in calculation of traveltimes and raypaths. Noting the point, we have thought a new method to remove the negative effect of mesh and to improve calculation accuracy by correcting the traveltimes with a little of additional calculation, and obtained better results.

保重心脉冲压缩滤波及其在速度自动拾取中的应用

To pick velocity automatically is not only helpful to improve the efficiency of seismic data process, but also to provide quickly the initial velocity for prestack depth migration. In this thesis, we use the Viterbi algorithm to do automatic picking, but the velocity picked usually is immoderate. By thorough study and analysis, we think that the Viterbi algorithm has the function to do quickly and effectually automatic picking, but the data provided for picking maybe not continuous on derivative of its curved surface, viz., the curved face on velocity spectrum is not slick. Therefore, the velocity picked may include irrational velocity information. To solve the problem above, we develop a new method to filter signal by performing nonlinear transformation of coordinate and filter of function. Here, we call it as Gravity Center Preserved Pulse Compressed Filter (GCPPCF). The main idea to perform the GCPPCF as follows: separating a curve, such as a pulse, to several subsection, calculating the gravity center (coordinate displacement), and then assign the value (density) on the subsection to gravity center. When gravity center departure away from center of its subsection, the value assigned to gravity center is smaller than the actual one, but non other than gravity center anastomoses fully with its subsection center, the assigned value equal to the actual one. By doing so, the curve shape under new coordinate breadthwise narrows down compare to its original one. It is a process of nonlinear transformation of coordinate, due to gravity center changing with the shape of subsection. Furthermore, the gravity function is filter one, because it is a cause of filtering that the value assigned from subsection center to gravity center is obtained by calculating its weight mean of subsetion function. In addition, the filter has the properties of the adaptive time delay changed filter, owing to the weight coefficient used for weight mean also changes with the shape of subsection. In this thesis, the Viterbi algorithm inducted, being applied to auto pick the stack velocity, makes the rule to integral the max velocity spectrum ("energy group") forward and to get the optimal solution in recursion backward. It is a convenient tool to pick automatically velocity. The GCPPCF above not only can be used to preserve the position of peak value and compress the velocity spectrum, but also can be used as adaptive time delay changed filter to smooth object curved line or curved face. We apply it to smooth variable of sequence observed to get a favourable source data ta provide for achieving the final exact resolution. If there is no the adaptive time delay-changed filter to perform optimization, we can't get a finer source data and also can't valid velocity information, moreover, if there is no the Viterbi algorithm to do shortcut searching, we can't pick velocity automatically. Accordingly, combination of both of algorithm is to make an effective method to do automatic picking. We apply the method of automatic picking velocity to do velocity analysis of the wavefield extrapolated. The results calculated show that the imaging effect of deep layer with the wavefield extrapolated was improved dominantly. The GCPPCF above has achieved a good effect in application. It not only can be used to optimize and smooth velocity spectrum, but also can be used to perform a correlated process for other type of signal. The method of automatic picking velocity developed in this thesis has obtained favorable result by applying it to calculate single model, complicated model (Marmousi model) and also the practical data. The results show that it not only has feasibility, but also practicability.

基于结构分析地震资料去噪方法研究

With the development of oil and gas exploration, the exploration of the continental oil and gas turns into the exploration of the subtle oil and gas reservoirs from the structural oil and gas reservoirs in China. The reserves of the found subtle oil and gas reservoirs account for more than 60 percent of the in the discovered oil and gas reserves. Exploration of the subtle oil and gas reservoirs is becoming more and more important and can be taken as the main orientation for the increase of the oil and gas reserves. The characteristics of the continental sedimentary facies determine the complexities of the lithological exploration. Most of the continental rift basins in East China have entered exploration stages of medium and high maturity. Although the quality of the seismic data is relatively good, this areas have the characteristics of the thin sand thickness, small faults, small range of the stratum. It requests that the seismic data have high resolution. It is a important task how to improve the signal/noise ratio of the high frequency of seismic data. In West China, there are the complex landforms, the deep embedding the targets of the prospecting, the complex geological constructs, many ruptures, small range of the traps, the low rock properties, many high pressure stratums and difficulties of boring well. Those represent low signal/noise ratio and complex kinds of noise in the seismic records. This needs to develop the method and technique of the noise attenuation in the data acquisition and processing. So that, oil and gas explorations need the high resolution technique of the geophysics in order to solve the implementation of the oil resources strategy for keep oil production and reserves stable in Ease China and developing the crude production and reserves in West China. High signal/noise ratio of seismic data is the basis. It is impossible to realize for the high resolution and high fidelity without the high signal/noise ratio. We play emphasis on many researches based on the structure analysis for improving signal/noise ratio of the complex areas. Several methods are put forward for noise attenuation to truly reflect the geological features. Those can reflect the geological structures, keep the edges of geological construction and improve the identifications of the oil and gas traps. The ideas of emphasize the foundation, give prominence to innovate, and pay attention to application runs through the paper. The dip-scanning method as the center of the scanned point inevitably blurs the edges of geological features, such as fault and fractures. We develop the new dip scanning method in the shap of end with two sides scanning to solve this problem. We bring forward the methods of signal estimation with the coherence, seismic wave characteristc with coherence, the most homogeneous dip-sanning for the noise attenuation using the new dip-scanning method. They can keep the geological characters, suppress the random noise and improve the s/n ratio and resolution. The rutine dip-scanning is in the time-space domain. Anew method of dip-scanning in the frequency-wavenumber domain for the noise attenuation is put forward. It use the quality of distinguishing between different dip events of the reflection in f-k domain. It can reduce the noise and gain the dip information. We describe a methodology for studying and developing filtering methods based on differential equations. It transforms the filtering equations in the frequency domain or the f-k domain into time or time-space domains, and uses a finite-difference algorithm to solve these equations. This method does not require that seismic data be stationary, so their parameters can vary at every temporal and spatial point. That enhances the adaptability of the filter. It is computationally efficient. We put forward a method of matching pursuits for the noise suppression. This method decomposes any signal into a linear expansion of waveforms that are selected from a redundant dictionary of functions. These waveforms are chosen in order to best match the signal structures. It can extract the effective signal from the noisy signal and reduce the noise. We introduce the beamforming filtering method for the noise elimination. Real seismic data processing shows that it is effective in attenuating multiples and internal multiples. The s/n ratio and resolution are improved. The effective signals have the high fidelity. Through calculating in the theoretic model and applying it to the real seismic data processing, it is proved that the methods in this paper can effectively suppress the random noise, eliminate the cohence noise, and improve the resolution of the seismic data. Their practicability is very better. And the effect is very obvious.

东南极格罗夫山土壤研究与古环境分析

Three soil spots were found in Grove Mountains, east Antarctica during 1999-2000, when the Chinare 16th Antarctic expedition teams entered the inland Antarctica. The characteristics of soils in Grove Mountains are desert pavement coating the surface, abundant water soluble salt, negligible organ matter, and severe rubification and salinization, scarces of liquid water, partly with dry permafrost, corresponding with the soils of McMurdo, Transantarctic. The soils age of Grove Mountains is 0.5-3.5Ma. Podzolization and redoximorphism are the main features in coastal Wilks region, in addition, there is strong enrichment of organic matter in many soils of this region. The main soil processes of Fildes Peninsula of King George Island include the intense physical weathering, decalcification and weakly biochemical processes. Peat accumulation is the main processes in Arctic because of humid and cold environment.Based on synthesis of heavy minerals, particle size, quartz grain surface textures, as well as pollen in soils, the soils parent materials of Grove Mountains derived from alluvial sediment of the weathering bedrocks around soils, and formed during the warm period of Pliocene. The detailed information is followed .l)The results of heavy minerals particle size showed the parent minerals derived form the weathering bedrocks around soils. 2)The quartz sand surface textures include glacial crushing and abrasion such as abrasive conchoidal fractures and grain edges, abrasive subparallel linear fractures and angularity, subaqueous environments produce V-shaped and irregular impact pits, polished surface, and chemical textures, such as beehive solution pits, which showed the water is the main force during the sediment of the soil parent minerals. 3)The pollen consist of 40 plant species, of which at least 5 species including Ranunculaceae, Chenopodiaceae, Artemisia, Gramineae, Podocarpus belong to the Neogene vegetation except the species from the old continent. Compared with Neogene vegetation of Transantarctic Mountains, Antarctic, we concluded that they grow in warm Pliocene.