地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

北天山新生代构造演化与古环境演变

The Tianshan Mountains is located about 1000-2000 km north of the India-Asia suture and is the most outstanding topography in central Asia, with transmeridional length of nearly 2500 km, north-southern wideness of ~ 300-500 km, peaks exceeding 7000 m above sea level (asl.), and average altitude of over 4000 m asl. Much of the modern relief of the Tianshan Range is a result of contraction driven by the collision of the India subcontinent with the southern margin of Asia, which began in early Tertiary and continues today. Understanding where, when and how the deformation of the Tianshan Mountains occurred is essential to decipher the mechanism of intracontinental tectonics, the process of foreland basin evolution and mountain building, and the history of climate change in central Asia. In order to better constrain the Cenozoic building history of the Tianshan Mountains and the climate change in the southern margin of the Junggar Basin, we carried out multiple studies of magnetostratigraphy, sedimentology, and stable isotopes of paleosol carbonate at the Jingou River section, which is located at the Huoerguosi anticline, the westernest one of the second folds and thrust faults zone in the northern piedmont of the Tianshan Mountains. The Jingou River section with a thickness of about 4160 m is continuous in deposits according to the observed gradual change in sedimentary environments and can be divided into five formations: Anjihaihe, Shawan, Taxihe, Dushanzi and Xiyu in upward sequence. Characteristic remamences were isolated by progressive thermal demagnetization, generally between 300 and 680℃. A total of 1133 out of 1607 samples yielded well-defined ChRMs and were used to establish the magnetostratigraphic column of a 3270-m-thick section from the exposed base of the Anjihaihe Formation to the middle of the Xiyu Formation. Two vertebrate fossil sites and a good correlation with the CK95 geomagnetic polarity time scale suggest that the section was deposited from ~30.5 to ~4.6 Ma and the age of the top of the Xiyu formation is ~2.6 Ma based on an extrapolation of the sedimentation rates. A plot of magnetostratigraphic age vs. height at the Jingou River section shows that significant increases in sedimentation rates as well as notable changes in depositional environments occurred at ~26-22.5 Ma, ~13-11 Ma and ~7 Ma, which represent the initial uplift of the Tianshan Mountains and two subsequent rapid uplift events. In addition, changes in sedimentation rates display characteristic alternations between increases and decreases, which probably indicate that the uplift of the Tianshan Mountains was episodic. We discussed the history of C4 biomass and climatic conditions in the southern margin of the Junggur Basin using the stable carbon and oxygen isotope composition of paleosol carbonates from the Jingou River section during ~17.5-6.5 Ma. The δ13C values indicate that the proportion of C4 biomass was uniform and moderate (15-20 %) during the interval of ~17.5-6.5 Ma. We proposed three hypotheses for this pattern of C4 biomass: (1) counteraction of two opposed factors (global cooling since ~15 Ma and thereafter increased dry and seasonality in central Asia) controlling the growth of C4 grasses, (2) variability in abundance of C3 grasses relative to C3 trees and shrubs if vegetation had ever changed in ecosystems, and (3) the higher latitude of the studied region. The δ18O values show a stepwise negative trend since ~13 Ma which may be attributed to three factors: (1) the temperature decreasing gradually after the middle Miocene (~15 Ma), (2) the increasing contribution of the moistures carried by the polar air masses from the Arctic Ocean to precipitation, and (3) the gradual retreat westward and disappearance of the Paratethys Ocean. Among them, which one played a more important role will need further study of the paleoclimate in central Asia.

典型人类活动对边坡变形及稳定性的影响研究

China is a mountainous country in which geological hazards occurred frequently, especially in the east of China. Except the geology, topography and extreme climate, the large scale human activities have become a major factor to landslides. Typical human activities which induced landslides are fill, cut and underground mining. On the topic of the deformation mechanism and slope stability, taking three different man-made slopes as examples, deformation mechanism and slope stability were studied by several methods, such as field work, numerical modeling and monitor. The details are as following: (1) The numerical modeling approach advantages over other conventional methods such as limit methods, so the numerical modeling is the major tool in this thesis. So far, there is no uniform failure criterion for numerical simulation. The failure criterion were summarized and analyzed firstly, subsequently the appropriate criterion was determinated. (2) Taking 220kV Yanjin transformation substation fill slope as example, the deformable characteristic, unstable mode and laboratory tests were studied systematically. The results show: the slope deformation was probably caused by a combination effect of unfavorable topographic, geological and hydro geological conditions, and external loading due to filling. It was concluded that the creep deformation of the slope was triggered by external loading applied at the back of the slope. In order to define the calculating parameters, a set of consolidated drained (CD) tests, consolidated undrained (CU) tests, repeated direct shear tests and UCS tests were carried out. The stability of the slope before and after reinforcement was assessed using 3D numerical modeling and shear strength reduction technique. The numerical modeling results showed: the factor of safety (FOS) of the slope was 1.10 in the natural state, and reduced to 1.03 after fill, which was close to the critical state and it caused creeping slip or deformation under rainfall. The failure surface in the slope is in active shear failure, whereas tensile failure occurs at the slope crest. After the site was reinforced with piles, the FOS was 1.27. Therefore, the slope is stable after reinforcement measures were taken. (3) The cut slope stability is a complex problem. Taking the left cut slope of Xiangjiaba as example in this thesis, the deformation and slope stability were studied systematically by numerical modeling and monitor methods. The numerical results show: the displacement is gradually increasing along with the cutting, and the largest displacement is 27.5mm which located at the bench between the elevation 340 and 380. Some failure state units distribute near the undermining part and there is no linked failure state occurred from crest to bottom during cutting. After cutting, some failure units appeared at the ground surface between elevation 340 and 360. The increasing tense stress made the disturbed rock failed. The slope is stable after cutting by the monitor method, such as surface monitor, multipoint displacement meter, inclinometer and anchor cable tensometer. (4) The interaction between underground mining and slope stability is a common situation in mountainous. The slope deformation mechanism induced by underground mining may contributed significantly to slope destabilization. The Mabukan slope in xiangjiaba was analyzed to illustrate this. Failure mechanism and the slope stability were presented by numerical modeling and residual deformation monitor. The results show: the roof deformed to the free face and the floor uplift lightly to the free face. The subsidence basin is formed, but the subsidence and the horizontal movement is small, and there is no failure zone occurred. When the underground mining is going on, the roof deformation, subsidence and the horizontal movements begin increasing. The rock deformation near the free face is larger than the ground surface, and the interaction between these coal seams appeared. There are some tensile failures and shear failures occurred on the roof and floor, and a majority of failure is tensile failure. The roof deformation, subsidence and the horizontal movements increased obviously along with the underground mining. The failure characteristic is shear failure which means the tensile stress transformed to the compressive stress. So the underground mining will induced tensile stress first which lead to structure crack, subsequently the compressive stress appeared which result in slippage. The crest was subjected to horizontal tension which made the rock crack along with the joint. The long term residual deformation monitor demonstrates that the slope is stable after the underground mining stopped.

地震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.

波动方程基准面校正方法研究

Datuming which has not been well solved in complex areas is a long-existing problem in seismic processing and imaging. Theoretically, Wave-equation datuming(WED) works well in the areas with substantial surface topography and areas of complex velocity structure. However, many difficulties still exist in practice. There are three main reasons: (1) It’s difficult to obtain the velocity model. (2) The computational cost is high and the efficiency is low. (3) Reflection waveform distortions are introduced by low S/N ratio in seismic data. The second and third problems are involved in the paper. To improve computational efficiency, DP1 proposed by Fu Li-Yun is applied in WED. Some quantitative and semi-quantitative conclusions of assessing the computational accuracy and efficiency have been obtained by comparing the adaptation of three operators( PS, SSF, DP1) to the surface topography and the lateral velocity variation. Moreover, the impacts of near surface scattering associated with complex surface topography on WED is analyzed theoretically. According to the analysis results, the following conclusions have been obtained. WED is stable and effective when the field data has high S/N ratio and velocity model is accurate. However ,it doesn’t work well when S/N ratio of field data is low. So denoising techniques in process of WED is important for low S/N data. The paper presents the theoretical analysis for the issues facing WED, which is expected to provide a useful reference to the further development of this technology.

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

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.

西南极南设得兰群岛中新生代火山岩基底沉积学及其大地构造制约

Livingston Island, the second island of South Shetland Island, constains Mesozoic-Cenozoic basement, Mesozoic-Cenozoic volcanic sequences, plutonic intrusions and post-subduction volcanic rocks, which document the history and evolution of an important part of the South Shetland Islands magmatic arc. The sedimentary sequence is named the Miers Bluff Formation (MBF) and is interpreted as turbidite since the first geological study on South Shetland Islands, and is interpreted as turbidite. It base and top are not exposed, but a thickness of more than 3000m has been suggested and seems plausible. The turbidite is overlain by Mid - Cretaceous volcanic rocks and intruded by Eocene tonalites. The age of the Miers Bluff Formation is poorly constrained Late Carboniferous -Early Triassic. Sedimentary Environment, tectonic setting and forming age of sedimentary rocks of the Miers Bluff Formation were discussed by means of the methods of sedimentology, petrography and geochemistry, combinig with the study of trace fossils and microfossil plants. The following conclusions are obstained. A sedimentary geological section of Johnsons Dock is made by outside measuring and watching, and then according the section, the geological map near the Spanish Antarctic station was mapped. Four pebbly mudstone layers are first distinguished, which thickness is about 10m. The pebbly mudstone is the typical rock of debris flow, and the depostional environment of pebbly mudstone may be the channel of mid fan of submarine fan. The sedimentsry structural characteristics and size analysis of sandstones show the typical sedimentary feature of turbidity flow and the Miers Bluff Formation is a deep-water turbidite (include some gravity-flow sediments). The materials of palaeocurrents suggest the continental slope dip to southeast, and indicate the provenance of turbidity sediment in the northwest area. By facies analysis, six main facies which include seven subfacies were recognized, which are formed in mid-fan and lower-fan of submarine, meanwhile, the sedimentary features of each facies and subfacies are summarized. The study of clastic composition, major elements, trace elements and rare earth elements indicates the forming setting of the Miers Bluff Formaton is active continental margin and continental island arc and the provenance is dissected magmatic arc which main composition is felsic gneiss. Many trace fossils of the whole succession were found in the turbidites of the Miers Bluff Formation. All these trace fossils are deep sea ichnofossils. There are fifteen ichnogenus, sixteen ichnospecies. Moreover, a new trace fossil was found and a new ichnogenus and new ichnospecies was proposed - Paleaichnus antarctics ichnogen, et ichnosp, nov.. Except the new ichnogenus and ichnospecies, others had been found in deep-sea flysch turbidites. Some are in mudstone and are preserved in the cast convex of overlying sandstone sole, they formed before turbidity flows occurred and belong to the high-different Graphoglyptida of fiysch mudstone. Others as Fucusopsis and Neonereites are preserved in sandstones and stand for trace assemblages after turbidity sedimentation. These trace fossils are typical members of abyssal "Nereites" ichnofacies, and provide for the depositional environment of the Miers Bluff Formation. Fairly diverse microfossil plants have been recovered from the Miers Bluff Formation, Livingston Island, including spores, pollen, acritarchs, wood fragments and cuticles. Containing a total of about 45 species (forms) of miospores, the palynofiora is quantitatively characterized by the dominance of non-striate bisaccate pollen, but spores of pteridophytes and pollen of gymnosperms are proportionate in diversity. It is somewhat comparable to the subzone C+D of the Alisporites zone of Antarctica, and the upper Craterisporites rotundus zone and the lower Polycingulatisporites crenulatus zone of Australia, suggesting a Late Triassic (possibly Norian-Rhaetian) age, as also evidenced by the sporadic occurrence of Aratrisporites and probable Classopollis as well as the complete absence of bisaccate Striatiti. The parent vegetation and paleoclimate are preliminarily deduced. At last, the paper prooses the provenance of sedimentary rocks of the Miers Bluff Formation locates in the east part to the southern Chile(or Southern South American). In the Triassic period, contrasting with New Zealand, Australia and South American of the Pacific margin of Gondwanaland, the Miers Bluff Formation is deposited in the fore-arc basin or back-arc basin of magmatic arc.

松辽徐家围子断陷层序地层充填模式及含油气性特征

Based on the principle and method of sequence, the author describes the sequence-filling model of the rifting basin of Xujiaweizi and its gas exploration potential. The object of this paper belongs to the area around Shengping-Wangjiatun anticline. Its srtatigraphy includes Huoshiling Formation (neutral and basic volcanic rocks), Shahezi Formation (coal bedding and mud and some sandstone) and Yingcheng Formation from bottom to top. These stratigraphy units are defined by author as mesosequences respectively. The author emphasizes that the main control factors of sequence change with the types of basin and stage of basin. So the sequence is researched according to the types of basin. This viewpoint is very new, and it is consistent with the principle of sequence. Volcanic action is very frequent and acute, topography difference is obvious. Between the volcanic events, Shahezi Formation is formed, which mainly consists-of sedimentary rocks. Based on the datum from seismic section and drilling core and well-logging, the author analyzes the single unit and unit set and system tract and sedimentary fancies, then, according to the accommodation space change and marking of sequence boundary, Shahezi Formation is divided into two Third-scale sequences. The sedimentary fancies and depth distribution are described. The author also pointed out that the volcanic rocks consume the accommodation space, so volcanic rocks can influence the development of sequence. Based on the concept of accommodation space, the author put volcanic rocks into sequence frame, which normally consists of sedimentary rocks. The topography of volcanic is controlled by lithology of volcanic rocks, the pattern of volcanic eruption and the topography before volcanic eruption. The topography of volcanic can influence sedimentation and the filling pattern of sedimentary rocks. The author describes the composition and lithology fancies and depth distribution of volcanic rocks. The volcanic rocks and Volcanic building, volcanic structure is recognized on seismic section. The author paid a special attention to the relationship between sedimentation and volcanism. Finally, the author analyses the combination of source-reservoir-cover unit in sequence frame. The mudstone of Shahezi Formation has a great depth, the Kerogene in it belongs to type II and III, which tends to produce gas. The Yingcheng Formation lies between Shahezi Formation and Denglouku Formation, belonging to good reservoir. The volcanic rocks of Huoshiling Formation often formed high building, which can capture the gas produced from Shahezi Formation. The stratigraphy of rifting basin of Xujiaweizi has the great potential of gas exploration. This paper claims the following creative points: 1. The author applied the principle and method of sequence to rifting basin, greatly extending its research area and topic issues. 2. The author pointed out that basin of different type and of different stage has a different type of sequence. This is caused by the different main control factors of sequence. 3. Put volcanic rocks into the sequence frame, discussing the probability of regarding the volcanic rocks as the component of sequence, dealing with the relationship between sedimentation and volcanism and its influence to the source-reservoir-cover system. 4. The author pointed out that the filling pattern of rifting basin are determined by the filling pattern of megasequence, whose filling pattern is determined by the filling pattern of system tract and the change of accommodation space.

大别山岩石圈有效弹性厚度初步研究

As powerful tools to study the lithosphere dynamics, the effective elastic thickness (Te) as well as the envelope of yielding stress of lithosphere have been attracted great attention of geoscientists in the past thirty years. The oceanic lithosphere, contrary to the continental lithosphere, has more fruits for its simple structures and evolution process. In continent, the lithosphere commonly is complex and variable in the rheological, thermal structures, and has a complicated history. Therefore, the application of the effective elastic thickness in continent is still a subject to learn in a long time. Te, with the definition of the thickness of an elastic plate in theory flexured by the equal benging of the real stress in the lithosphere plate (Turcotte, 1982), marks the depth of transition between elastic and fluid behaviors of rocks subjected to stress exceeding 100 MPa over the geological timescales (McNutt, 1990). There are three methods often adapted: admittance or isostatic response function, coherence and forwarding. In principle, the models of Te consist of thermal-rheological, non-linear Maxwell, non-linear work hardening and rheological layered models. There is a tentative knowledge of Te that it is affected by the following factors: crustal thickness, crust-mantle decoupling, plate bending, boundary conditions of plate (end forces and bending moments), stress state, sedimentary layer, faulting effect, variation in the mountain belts' strike, foreland basin, inheritance of tectonic evolution, convection of mantle, seismic depth and lithosphere strength. In this thesis, the author introduces the geological sketch of the Dabie collisional orogenic belt and the Hefei Basin. The Dabie Mts. is famous for the ultra-high pressure metamorphism. The crustal materials subducted down to the depth of at least 100 km and exhumed. So that the front subjects arise such as the deeply subduction of continent, and the post-collisional crust-mantle interaction. In a geological journey at June of 1999, the author found the rarely variolitic basaltic andesite in the Dabie Mts. It occurs in Susong Group, near Zhifenghe Countryside, Susong County, Anhui Province. It is just to the south of the boundary between the high-grade Susong melange and the ultra-high grade South Dabie melange. It has a noticeable knobby or pitted appearance in the surface. The size of the varioles is about 1-4 mm. In hand-specimen and under microscope, there are distinct contacts between the varioles and the matrice. The mineralogy of the varioles is primarily radiate plagioclase, with little pyroxene, hornblende and quartz. The pyroxene, hornblende and quartz are in the interstices between plagioclase. The matrix is consisted of glass, and micro-crystals of chlorite, epidote and zoisite. It is clearly subjected and extensive alteration. The andesite has an uncommon chemical composition. The SiO_2 content is about 56.8%, TiO_2 = 0.9%, MgO = 6.4%, (Fe_2O_3)_(Total) = 6.7% ～ 7.6%, 100 Mg/(Mg+Fe) = 64.1 ～ 66.2. Mg# is significantly high. The andesite has higher abundances of large-lithophile trace elements (e.g. K, Ba, Sr, LREE), e.g. La/Nd = 5.56-6.07, low abundances of high-strength-field elements (HFSE, e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subducted-related magmas (Pearcce, 1982; Sun and McDonaugh, 1989). In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a characteristic of the continental margins (Pearce, 1982). There has a strongly enrichment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)_N is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)_N = 28.63-36.74, (La/Y)_N = 70.33 - 82.84. The elements Y and Yb depleted greatly, Y < 20 ppm, Y_N = 2.74-2.84, Yb_N = 2.18 - 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is -18.7 ～ -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the basaltic andesite was resulted from the post-collisional crust-mantle interaction between the subducted Yangze carton and the mantle of Sino-Korea carton. To obtain the Te of the lithosphere beneath the Dabie Mts. and the Hefei Basin, the author applies the coherence method in this thesis. The author makes two topography-gravity profiles (profiles 7~(th) and 9~(th)) across the Dabie Mts. and the Hefei Basin, and calculates the auto-coherence, across coherence, power spectrum, across power spectrum of the topography and gravity of the two profiles. From the relationships between the coherence and the wave-number of profiles. From the relationships between the coherence and the wave-number of profiles 7~(th) and 9~(th), it is obtained that the characteristic wavelengths respectively are 157 km and 126 km. Consequently the values of effective elastic thickness are 6.5 km and 4.8 km, respectively. However, the Te values merely are the minimum value of the lithosphere because the coherencemethod in a relative small region will generate a systemic underestimation. Why there is a so low Te value? In order to check the strength of the lithosphere beneath the Dabie Mts., the authore tries to outline the yielding-stress envelope of the lithosphere. It is suggested that the elastic layers in the crust and upper mantle are 18 km and 35 km, respectively. Since there exist a low viscosity layer about 3-5 km thickness, so it is reasonable that the decoupling between the crust and mantle occurred. So the effective thickness of the lithosphere can be estimated from the two elastic layers. Te is about 34 km. This is the maximum strength of the lithosphere. We can make an approximately estimation about the strength of the lithosphere beneath the Dabie Mts.: Te is about 20-30 km. The author believes that the following factors should be responsible for the low Te value: (1) the Dabie Mts. has elevated strongly since K_3-J_1. The north part of the Dabie Mts. elevates faster than the south part today; (2) there occur large active striking faults in this area. And in the east, the huge Tan-Lu striking fault anyway tends to decrease the lithosphere strength; (3) the lithosphere beneath the Dabie Mts. is heter-homogeneous in spatio-temporal; (4) the study area just locates in the adjacent region between the eastern China where the lithosphere thickness is significantly reduced and the normal western China. These factors will decrease the lithosphere strength.

流体地球化学：成因、流动和流体-岩石相互作用及塔里木盆地实例研究

Geofluid in sedimentary basins is related to petroleum generation, migration, accumulation and preservation, and is a topic of geological frontier. By integrating the multi-discipline methods of petroleum geochemistry, sedimentology, hydrogeology, petroleum geology and experimental geochemistry, the thesis has carried out experiments of microcline dissolution in solutions with organic acids, crude oil, brines with high total dissolved solids (TDS), and has dealt with Al distribution between the crude oil and the brines after the experiments. Cases for study includes Central Tarim, Hetianhe Gas Field and Kucha forland basin with data containing fluid chemistry and isotopic compositions, thin sections of sandstones and carbonates, homogenization temperatures and salinities of fluid inclusions, isotopic compositions of bulk rock and autigenic minerals. The aims are to elucidate fluid origin and flow in the three areas, effect of hydrocarbon emplacement on diagenesis, and to show occurrence of microbe-mediated, and thermochemical sulfate reduction in the Tarim Basin. Microcline dissolution experiments show that after 100 hour, part of the dissolved Al distributes in the crude oil, and the Al concentrations in the crude oil rise when organic acids are added. The result can be used to explain that most oilfield waters in the Tarim Basin are characterized by less than 3mg/L Al. Crude oil added to the solutions can enhance microcline dissolution, which is also observed in the case - Silurian sandstones with early crude oil emplacement in the Central Tarim. Al and Si have higher concentrations in the experiments of oxalic acid than of acetic acid under the same pH conditions, suggesting that there exist Al-oxalate and Si-oxalate complexes. Presence of acetate can enhance the activity of Ca and Al, but Al concentrations have not been increased significantly due to formation of small Al-acetate complex during the experiments. Relationships between δD and δ~(18)O in conjunction with chemistry of oilfield waters show that the waters are evaporated connate waters, which subsequently mixed with meteoric water, and were influenced by water-rock interactions such as salt dissolution, dolomitization of calcite, albitization of feldspar. In the Hetianhe Gas Field where salt dissolution took place, δD and δ~(18)O values can be used to trace nicely meteoric water recharge area and flow direction, but TDS can not. Part of the waters have high TDS but very light δD and δ~(18)O. When combined with paleo-topography, or fluid potentials, meteoric water is suggested to flow eastward in the Hetianhe Gas Field, which is the same with the Central Tarim. Whist in the Kuche forland basin, meteoric water may have permeated Cambrian-Ordovician strata. Relationship between ~(87)Sr/~(86)Sr and 1/Sr can be used to indicate migration and mixing of brines from carbonate strata (low ~(87)Sr/~(86)Sr ratio but high Sr content), clastic strata (high ~(87)Sr/~(86)Sr ratio but low Sr content) and crystalline basement (high ~(87)Sr/~(86)Sr ratio and heavy δ~(18)O value). Using this approach, it can be found that ~(87)Sr-depleted brine from Ordovician carbonates have migrated up to and mixed with ~(87)Sr-enriched waters from Silurian and Carboniferous sandstones, and that Silurian brines have mixed with meteoric water. In the Kuche forland basin, brines from the Cambrian and Ordovician carbonates have higher ~(87)Sr/~(86)Sr ratios than those from the overlying sandstones, when combined with chemistry, δ~(15)N and ~3He/~4He ratios of the coexisting natural gases, suggesting that the brines were derived from the basement. There exists some debate on the effect of hydrocarbon emplacement on mineral diagenesis. Case-study from Silurian sandstones in the Central Tarim show that quartz has kept overgrowing secondarily when oil saturation was decreased by meteoric water flushing subsequently to hydrocarbon emplacement. Silicon precipitates on the water-wet quartz surface, leading to decreased Si concentration close to the surface. A Si grads can result in Si diffusion, which supplies Si for quartz overgrowth. Hydrocarbon oxidation-sulfate reduction is an important type of organic-inorganic interaction. Not only can it make secondary alteration of hydrocarbons, but generate H_2S and CO_2 gases which can improve reservoir property. Thermochemical sulfate reduction took place at the temperatures more than 125 ℃ to 140 ℃ in the Cambrian-Ordovician carbonates, the products - H_2S and CO_2 gases migrated up to the Silurian, and precipitated as pyrite and calcite, respectively. The pyrite has an average δ~(34)S value close to those of Ordovician seawater and anhydrite, and calcite has δ~(13)C value as low as -21.5‰. In the Hetianhe Gas Field, sulfate reduction bacteria carried by meteoric water flowing eastward may have preferentially depleted ~(12)C of light hydrocarbon gases, and results in heavier δ~(13)C values of the residual hydrocarbon gases and higher molar CO_2 in the natural gases in the west than in the east. Coexisting pyrite has δ~(34)S values as low as -24.9‰.

油气运聚动力学的模拟分析——以准噶尔盆地为例

The petroleum migration, happening in the geologic past, is the very important and complex dynamic processes in the petroleum systems. It plays a linking role among all static factors in a system. The accumulation is in fact the result of the petroleum migration. For the petroleum geology, the dynamics research of the petroleum migration refers to the mechanism and process research, as well as the use of the quantitative methods. In this thesis, combining with the qualitative analysis and quantitative modeling, the author manages to discuss theoretically some key problems dealing with migration processes, which have not been solved yet, and to apply the studied results in petroleum system analysis in actual basins. The basin analysis offers the base of the numerical modeling for geological phenomena occurring in sedimentary basins, that consists of the sedimentary facies analysis, the section reconstructing technique, eroded thickness estimating, etc. The methods to construct the geologic model, which is needed in the research of oil and gas migration and accumulation, are discussed. The basin analysis offers also the possibility for the latter modeling works to get and select the parameters, such as stratum's thickness, age, stratigraphy etc. Modeling works were done by using two basin modeling softwares: Basin_Mod and TPC_Mod. The role of compaction during the secondary migration and the heterogeneity of migrating paths within the clastic carrier are modeled. And the conclusions were applied in the migration studies in the Jungaer Basin, lying on the Northwest part of the China. To construct a reliable migration model, the author studied the characteristics of the sedimentation, the pore fluid pressure evolution, as well as the distribution and the evolution of fluid potential, following the tectonic evolution of the Jungaer Basin. The geochemical prospecting results were used to evidence and to calibrate the migration processes: the oil-source correlation, the distribution of the properties of oil, gas and water. Finally, two important petroleum systems, Permian one and Jurassic one were studied and identified, according, principally, to the studies on the petroleum migration within the Jungaer Basin. Since the oil, as well as the gas, moves mainly in separate phase during the secondary migration, their migrating behaviors would be determined by the dynamics conditions of migration, including the driving forces and pathways. Based on such a consideration, the further understandings may be acquired: the roles played by permeable carriers and low-permeable source rock would be very different in compaction, overpressure generation, petroleum migration, and so on. With the numerical method, the effect of the compaction on the secondary migration was analyzed and the results show that the pressure gradient and the flux resulted from compaction are so small that could be neglected by comparing to the buoyancy of oil. The main secondary migration driving forces are therefore buoyancy and capillary within a hydrostatic system. Modeling with the commercial software-Basin_Mod, the migration pathways of petroleum in clastic carriers seem to be inhomogeneous, controlled by heterogeneity of the driving force, which in turn resulted from the topography of seals, the fabrics and the capillary pressure of the clastic carriers. Furthermore, the direct and indirect methods to study fault-sealing properties in the course of migration were systemically summarized. They may be characterized directly by lithological juxtaposition, clay smear and diagenesis, and indirectly the comparing the pressures and fluid properties in the walls at two apartments of a fault. In Jungaer Basin, the abnormal pressures are found in the formations beneath Badaowan or Baijantan Formation. The occurrence of the overpressure seems controlled by the stratigraphy. The rapid sedimentation, tectonic pressuring, clay sealing, chemical diagensis were considered as the principal pressuring mechanisms. The evolution of fluid pressure is influenced differently at different parts of the basin by the tectonic stresses. So the basin appears different pressure evolution cycles from each part to another during the geological history. By coupling the results of thermal evolution, pressure evolution and organic matter maturation, the area and the period of primary migration were acquired and used to determine the secondary migration time and range. The primary migration in Fengcheng Formation happened from latter Triassic to early Jurassic in the main depressions. The main period of lower-Wuerhe Formation was at latter Jurassic in Changji, Shawan and Pen-1-jing-xi Depression, and at the end of early Cretaceous in Mahu Depression. The primary migration in Badaowan and Sangonghe Formation is at the end of early-Cretaceous in Changji Depression. After then, the fluid potential of oil is calculated at the key time determined from area and time of the primary migration. Generally, fluid potential of oil is high in the depressions and low at the uplifts. Synthetically, it is recognized that the petroleum migration in the Jungaer Basin is very complex, that leads us to classify the evolution of petroleum systems in Northwestern China as a primary stage and a reformed one. The remigration of accumulated petroleum, caused by the reformation of the basin, results in the generation of multiple petroleum systems. The faults and unconformities are usually the linkers among the original petroleum systems. The Permian petroleum system in Jungaer Basin is such a multiple petroleum system. However, the Jurassic petroleum system stays still in its primary stage, thought the strong influences of the new tectonic activities.

塔里木盆地奥陶系不整合面岩溶作用特征与油气储层研究

Based on the principle and methods of carbonate sedimentology and reservoir geology, and guided by the theories of carbonate reservoir geology, the palaeokarst of Ordovician carbonate rocks in Tarim Basin has been comprehensively studied with multiple methods from different branches of geology. It is indicated that the features and distribution of palaeokarstification have developed in Ordovician carbonates. The controlling of karstification to Ordovician carbonate reservoirs has been discussed. Regional distribution of carbonate reservoirs controlled by karstification has been predicted within this basin. The main consents and conclusions of the this dissertation is as follows: Nine key indicators to the recognition of palaeokarst are proposed in terms of careful observation upon the well cores, lithological and geochemical analyses, and drilling and logging responses to the karst caves and fractures. The time and environment of cave filling are documented from careful research of lithofacies, mineralogy, and geochemistry of the physical and chemical fillings within karst caves. The caves in Ordovician carbonates were filled in Early Carboniferous in Lunnan area. The muddy filling in upper caves was deposited under subaerial fresh-water setting, while the muddy filling in lower caves was formed in the mixed water body of fresh-water and dominated sea water. Although most chemical fillings are suggested being precipated in the burial diagenetic environment after karstification but mineralogic and geochemical characteristics of some chemical fillings indicates they formed in meteoric environment during the karstification. It is obvious that the palaeokarst has been zoned in vertical profile. It can be divided into four units from top to bottom: surface karst, vadose karst, phreatic and tranquil flow zones. Between two types of limestone karst and dolostone karst are firstly differentiated in Tarim Basin, based on the comparison of features of each karst zone in limestone and dolostone regions. In Tabei area, the lowest depth of karstification is approximately 300 m below the Upper Ordovician unconformity interface, while the bottom depth of karstification in Tazhong area ranges commonly from 300 to 400 m, in rare cases may be up to 750 m below the upper Ordovician unconformity interface. In Lunnan and Tazhong areas, the palaeokarst morphology and the surface hydrosystem are firstly reconstructed based on the top of carboniferous "Shuangfeng limestone bed (Double-Peaks limestone)" as basal. According to the palaeomorphologic feature, karst topography can be divided into three units: karst upland, karst slope, and karst valley. Vadose zone was well developed in karst upland, and it can be found in a quite depth. Both vadose and phreatic zones were well developed in karst slope and upstream valley. In downstream valley, the karstification is not strong, the vadose and phreatic zones are thin in thickness. In Tazhong and Yingmaili areas, karstification is also developed in relict carbonate palaeo-hills which existed as isolated blocks admits clastic strata.

工程地质三维建模与可视化技术应用开发

Synthetic Geology Information System (SGIS) is an important constituent part of the theory of Engineering Geomechanics Mate-Synthetic (EGMS), and is the information system more suited for the collection, storage, management, analysis and processing to the information coming from engineering geology,' geological engineering and geotechnical engineering. Its contents involve various works and methods of the investigation, design, and construction in different stages of the geological engineering. Engineering geological and three-dimensional modeling and visualization is the fundamental part of the SGIS, and is a theory, method and technique by which, adopting the computer graphics and image processing techniques, the data derived from engineering geological survey and the calculated results obtained from the geomechanical numerical simulation and analysis are converted to the graphics and images displayed on the computer screen and can be processed interactively. In this paper, the significance and realizing approaches of the three-dimensional modeling and visualization for the complex geological mass in the engineering geology are discussed and the methods of taking advantage of the interpolation and fitting for the scattered and field-surveyed data to simulate the geological layers, such as the topography and earth surface, the groundwater table and the stratum boundary, are researched into. At the mean time, in mind the characteristics of the structure of the basic data for three-dimensional modeling, its visual management can be resolved into the engineering surveyed database management module, plot parameter management module and data output module and the requirement for basic data management can be fulfilled. In the paper, the establishment and development of the three-dimensional geological information system are probed tentatively, and an instance of three-dimensional visual Engineering Distribution Information System (EDIS), theConstruction Management Information System for an airport, in which the functions, such as the real-time browse among the three-dimensional virtual-reality landscapes of the airport construction from start to finish, the information query to the airport facility and the building in the housing district and the recording and playback of the animation sets for the browse and the takeoff and landing of the planes, is developed by applying the component-mode three-dimensional virtual-reality geological information system (GIS) software development kits (SDK), so the three-dimensional visual management platform is provided for the airport construction. Moreover, in the gaper, integrated with the three-dimensional topography visualization and its application in the Sichuan-Tibet Highways, the method of the digital elevation model (DEM) data collection from the topographic maps is described, and the three-dimensional visualization and the roaming about the terrain along the highway are achieved through computer language programming. Understanding to the important role played by the varied and unique topographical condition in the gestation and germination of the highly-dense, frequently-arising and severely-endangered geological hazards can be deepened.

利用重力与地形数据探讨我国大陆若干地区岩石圈有效弹性厚度

The lithosphere behaves as an elastic thin plate overlying a weak fluid. When loaded by topography, it deflects to compensate the topography and isostasy is achieved. The response of the plate is often characterized by the flexural rigidity, or equivalently, by the effective elastic thickness(EET). The relationships between gravity and topography have commonly been employed to investigate quantitively the isostatic compensation of the lithosphere. Cross-spectral techniques have been developed to estimate the admittance and coherence between gravity and topography. The observed admittance and coherence functions, when compared with theoretical admittance and coherence computed for various models and parameters, can provide an estimate of the effective elastic thickness. As for coherence, the wavelength at which the coherence drops from approaching unity at long wavelengths to value approaching zero at short wavelengths is a measure of the rigidity of the lithosphere. This research takes advantage of the high-resolution gravity and topography data in Erdos. Using the coherence technique, we have estimated the effective elastic thickness of this region. Subsurface loads are also considered in our calculation. When the admittance between topography and gravity is obtained, we can filter the topography to give an estimate of gravity anomaly. This would be very important when the gravity data is scarce. Several other regions of China have been selected to investigate the effective elastic thickness of lithosphere. We compare those result with lithosphere thickness obtained through seismological technique and heat flow of the region. We find the effective elastic thickness is always smaller than the seismogenic thickness. This is due to what the effective elastic thickness reflects is only the up elastic part of the lithosphere. And we also find there is some degree of correlation between effective elastic thickness and heat flow. This suggests EET is probably controlled by the thermal state of the lithosphere and correlates with their tectonic age. Thus the estimations of the effective elastic thickness of lithosphere can help to investigate geophysical features of active tectonics of continental lithosphere.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN