基于弹性参数分析的碳酸盐岩储层预测方法研究

Prediction of Carbonate Reservoir Based on the Elastic Parameter Analysis Zhang Guangzhi (Solid Geophysics) Directed by Professor Liu Hong Abstract With the exploration and development of Puguang Oilfield, oil-gas exploration of carbonate rock in China has shown good prospects. Research on earthquake prediction methods for carbonate reservoir becomes the key of oil and gas exploration. Starting with analysis of geological characteristics of carbonate rock, prestack AVO inversion method, prestack elastic impedance inversion and parameter calculation method and seismic attribute extraction and optimization method were studied based on the analysis of rock physics in this work. First, variation characteristic and law of carbonate rock reservoir parameters were studied based on experimental data of rock physics, log data, analysis assay data, mud logging data and seismic data, so as to lay a foundation for the further reservoir identification and description. Then, the structure, type and propagation law of seismic wave field were analyzed through seismic forward modeling of the reservoir, and contact between information from log and geology data with elastic parameters, such as compressional wave and shear wave velocity and density were established, so as to provide a standard for reservoir identification and hydrocarbon detection using seismic reflection characteristics of the research area. Starting with the general concept of inverse problem, through analysis of Zoeppritz equation, three kinds of pre-stack inversion methods were derived and analyzed in detail, the AVO 3-parameter inversion based on Bayesian theory, the prestack AVO waveform inversion method and the simultaneous inversion method, based on the statistical hypothesis of inversion parameters and observation data and the Gauss distribution assumption of noise. The three methods were validated by model data and real data. Then, the elastic wave impedance inversion method of carbonate reservoir was investigated and the method of elastic parameter extraction from elastic impedance data was put forward. Based on the analysis of conventional methods of seismic attribute extraction and optimization, the time-frequency attributes and the wavelet attributes with time and amplitude feature were presented, and the prestack seismic attribute calculation method which can characterize the reservoir rock and fluid characteristic was presented. And the optimization of seismic attribute using the nonlinear KPCA method was also put forward. A series of seismic prediction technologies for carbonate reservoir were presented based on analysis of rock physics and seismic forward simulation technology. Practical application of these technologies was implemented in A oil field of Southern China and good effect has been achieved. Key words: carbonate rock; reservoir prediction; rock physics, prestack seismic inversion; seismic attribute

青藏高原横波速度与各向异性结构及其壳幔形变印迹

As the most spectacular and youngest case of continental collision on the Earth, to investigate the crust and mantle of Tibetan plateau, and then to reveal its characters of structure and deformation, are most important to understand its deformation mechanism and deep process. A great number of surface wave data were initially collected from events occurred between 1980 and 2002, which were recorded by 13 broadband digital stations in Eurasia and India. Up to 1,525 source-station Rayleigh waveforms and 1,464 Love wave trains were analysed to obtain group velocity dispersions, accompanying with the detail and quantitative assessment of the fitness of the classic Ray Theory, errors from focal and measurements. Assuming the model region covered by a mesh of 2ox2o-sized grid-cells, we have used the damped least-squares approach and the SVD to carry out tomographic inversion, SV- and SH-wave velocity images of the crust and upper mantle beneath the Tibetan Plateau and surroundings are obtained, and then the radial anisotropy is computed from the Love-Rayleigh discrepancy. The main results demonstrate that follows, a) The Moho beneath the Tibetan Plateau presents an undulating shape that lies between 65 and 74 km, and a clear correlation between the elevations of the plateau and the Moho topography suggests that at least a great part of the highly raised plateau is isostatically compensated. b) The lithospheric root presents a depth that can be substantiated at ~140 km (Qiangtang Block) and exceptionally at ~180 km (Lhasa Block), and exhibits laterally varying fast velocity between 4.6 and 4.7 km/s, even ~4.8 km/s under northern Lhasa Block and Qiangtang Block, which may be correlated with the presence of a shield-like upper mantle beneath the Tibetan Plateau and therefore looked as one of the geophysical tests confirming the underthrusting of India, whose leading edge might have exceeded the Bangong-Nujiang Suture, even the Jinsha Suture. c) The asthenosphere is depicted by a low velocity channel at depths between 140 and 220 km with negative velocity gradient and velocities as low as 4.2 km/s; d) Areas in which transverse radial anisotropy is in excess of ~4% and 6% on the average anisotropy are found in the crust and upper mantle underlying most of the Plateau, and up to 8% in some places. The strength, spatial configuration and sign of radial anisotropy seem to indicate the existence of a regime of horizontal compressive forces in the frame of the convergent orogen at the same time that laterally varying lithospheric rheology and a differential movement as regards the compressive driving forces. e) Slow-velocity anomalies of 12% or more in southern Tibet and the eastern edge of the Plateau support the idea of a mechanically weak middle-to-lower crust and the existence of crustal flow in Tibet.

地下水位上升造成的黄土远程滑坡形成机理与数值模拟研究

Based on field survey, laboratory testing and numerical modeling, engineering characteristics of undisturbed loess and the mechanism of long-runout loess landslides caused by underground water level rise, as well as the formation conditions and spatial distribution of landslides, are systematically studied and analyzed. Loess landslides at south Plateau of Jingyang County are mainly classified as flowslide, slide and fall. Flowslide is the main type characteristic of high velocity, long runout and multi-stages. The steep relief composed of loose structured loess-old aged soil serials and the rise of groundwater table are the predominant conditions for landslides in the study area. To study loess mechanic poperties and loess landslides mechanisims, isotropically and anisotropically consolidated undrained compression(ICU and ACU) tests and constant-deviator-drained compression (CQD) tests were carried out on undisturbed samples. The results of undrained compression tests performed at the in-situ stress level show that the soils are of consistently strain-softening in the stress-strain relations and cause high excess pore pressure. The steady-state line and the potential region of instability are obtained from ICU and ACU test results. A necessary condition for liquefaction is that the soil state initially lies in or is brought into the potential instability region. In addition, a strong strain-softening model is also formed. CQD tests demonstrate that the mobilized friction angle is far less than the steady-state angle and that the soil experiences undrained contractive failure suddenly at very small strains when its stress path during drained loading tries to cross the potential instability region，thus validates the proposed instability region. Based on the location of the region of potential instability and the stress state of slope soil, a method of static liquefaction analysis is proposed for loess landslides caused by rise in groundwater table. Compared with other liquefaction analysis methods, this method overcomes the limitations inherent in conventional slope stability method and undrained brittleness index method. Triaxial tests composed of constant water content (CW) and wetting tests at constant deviator stress are performed on undisturbed unsaturated samples. The stress-strain relation of CW tests takes on strain-hardening behavior; The results of wetting tests at constant deviator stress designed to study the mechanics of failure of unsaturated loess caused by an increase in the degree of saturation (wetting) shows that a contractive failure occurs in the undisturbed samples. On the basis of the above triaxial test results, the initiation of static liquefaction is presented for long-runout loess landslides caused by rise in groundwater table, that is, the loess slope soil gradually transfer from unsaturated to saturated state under the infiltration of irrigation. A contractive failure occurs in the local region at very small strain by increasing the pore-water pressure at constant deviator stresses under drained conditons. It is the contractive failrue resulting from rise of pore pressure that leads to high excess pore pressure in the neighbour soil which reduces shear resistance of soil. The neighbour soils also fail due to the rapid increase in pore-water pressure. Thus a connected failure surface is developed quickly and a flowslide occurs. Based on the saturated-unsaturated seepage theory, transient seepage is computed using the finite element method on loess slope under groundwater table rise. Pore-water pressure distribution for every time step after irrigation are obtained. The phreatic surface in the slope increases with the groundwater table. Pore-water pressure distribution within 8m above the phreatic surface changes very quickly，but the water content and pore water pressure in the region ranging from 8m above the phreatic surface up to ground surface is almost not affected and the matric suction usually is kept at 100～120 kPa. Based on the results of laboratory tests and seepage flow analysis, the development process of loess landslide is modeled considering groundwater table rise. The shearing plastic zone first occurs at the slope toe where the soil is soaked for long term during rise in groundwater table. As irrigation continues, the shearing plastic zone gradually extends to the interior soils, with the results that the tensile plastic zone occurs at the slope crown. As time goes on, both the shearing plastic zone and tensile plastic zone continue to extend. Then a connected plastic zone is formed and fowslide occurs. In comparision to laboratory test results, the results of numerical simulation quite well verify the presented mechanism of static liquefaction of long-runout loess landslides caused by rise in groundwater table.

KEL-TI介质的地震波场响应与成像方法研究

To deal with the problems in multi-component converted seismic wave exploration in coal fields, the wave propagating features and imaging methods of multi-component converted waves in coal measure strata are researched in this thesis firstly. The relations between viscoelasticity and anisotropy in coal measure strata are analyzed to build KEL-TI model, and which seismic wave propagating and attenuating features are researched. The disadvantages of converted wave imaging methods based on common converted point gather are analyzed and constant velocity no NMO converted wave imaging method based on common scattering point gather is put forward, according to Huygens-Fresnel principle, which applicabilities in the elastic isotropic, elastic TI and KEL-TI situations are discussed. To different model simulation data, the common scattering point gathers’ and stacked profiles’ features are analyzed. The results show that the method can image compressional waves and converted waves with high precision. Secondly, the resolution enhancing theories and methods of converted wave are researched by Rayleigh wave suppressing, converted wave static correction and poststack inverse-Q filtering. 1) The polarization filter is designed by the instantaneous polarization information of seismic waves, and the Rayleigh wave suppressing method is researched. From the spectrum analysis before and after filtering, it can be derived that the amplitudes are kept relatively. 2) To constant velocity no NMO converted wave imaging method, the static correction method based on common equivalent offset point gather is put forward and tested to the actual converted waves. 3) The relation between equivalent quality factor of converted wave, compressional wave quality factor and the ratio of compressional to shear wave velocity is derived. The compressional wave quality inversion method by first arrivals of none-offset VSP is researched, and which is then transformed to the equivalent quality factor to perform inverse-Q filtering of actual converted waves. The result has shown that the method can recover the high frequency energy of converted waves. At last, the theories and methods researched in this thesis are practiced to the 3D3C seismic exploration in Guqiao coal mine in Huainan and achieve good results.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

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

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

非饱和膨胀土力学特性及其工程处置研究

The unsaturated expansive soil is a hotspot and difficulty in soil mechanics inland and outland. The expansive soil in our China is one of the widest in distributing and greatest in area, and the disaster of expansive soil happens continually as a result. The soil mechanics test, monitor, numerical simulation and engineering practice are used to research swell and shrinkage characteristic, edge strength characteristic and unsaturated strength characteristic of Mengzi expansive soil. The seep and stability of the slope for expansive soil associated with fissure are analyzed and two kinds of new technique are put forward to be used in expansive soil area, based on disaster mechnics proposed of the slope.The technique of reinforcement in road embankment is optimized also. Associated with engineering geology research of Mengzi expansive soil, mineral composition, chemical composition, specific area and cation content, dissolubility salt and agglutinate, microcosmic fabric characteristic, cause of formation and atmosphere effect depth are analyzed to explain the intrinsic cause and essence of swell and shrinkage for expansive soil. The rule between swell-shrinkage and initial state, namely initial water content, initial dry density and initial pressure, can be used to construction control. Does Response model is fit to simulate the rule, based on ternary regression analysis. It has great meaning to expansive soil engineering in area with salt or alkali. The mechanics under CD, CU and GCU of expansive soil is researched by edge surface theory to explain the remarkable effect of consolidation pressure, initial dry density, initial water content, cut velocity, drainage and reinforcement to the edge strength characteristic. The infirm hardening stress strain curves can be fitted with hyperbola model and the infirm softening curves can be fitted with exponential model. The normalization theory can be used to reveal the intrinsic unity of the otherness which is brought by different methods to the shear strength of the same kinds of samples. The unsaturated strain softening characteristic and strength envelope of remolding samples are researched by triaxial shear test based on suction controlled, the result of which is simulated by exponential function. The strength parameters of the unsaturated samples are obtained to be used in the unsaturated seep associated with rainfall. The elasticity and plasticity characters of expansive soil are researched to attain the model parameters by using modified G-A model. The humidification destroy characteristic of expansive soil is discussed to research the disaster mechanism of the slope with the back pressure increasing and suction decreasing under bias pressure consolidation. The indoor and outdoor SWCCs are measured to research the effect factors and the rule between different stress and filling environment. The moisture absorption curves can express the relationship between suction and water content in locale. The SWCCs of Mengzi expansive soil are measured by GDS stress path trixial system. The unsaturated infiltration function is gained to research seep and stability of the slope of expansive soil. The rainfall infiltration and ability of slope considering multifarious factors are studied by analyzing fissure cause of Mengzi expansive soil. The mechanism of the slope disaster is brought forward by the double controlling effect between suction and fissure. Two new kinds of technique are put forward to resolve disaster of expansive soil and the technique of reinforcement on embankment is optimized, which gives a useful help to solving engineering trouble.

填土边坡水文作用机理与破坏过程研究—以深圳羊宝地填土滑坡为例

A full understanding of failure mechanism, critical hydrological condition, and process of mobilization and deposition of a landslide is essential for optimal design of stabilization measure and forecasting of landslide hazard. This requires a quantitative study of hydrological response of a slope to rainfall through field monitoring, laboratory test and numerical modelling. At 13:40 on September 18, 2002, a fill slope failed following a period of prolonged rain in Shenzhen, resulting in 5 fatalities and 31 injuries. The failed mass with a volume about 2.5×104m3 traveled about 140m on level ground. Field monitoring, laboratory test, theoretical analysis and numerical modelling were carried out to undestand the hydrological response and failure mechanism of this fill slope. This thesis mainly focuses on the following aspects: (1) The hydrological responses and failure processes of slopes under rainfall infiltration were reviewed. Firstly, the factors influencing on the hydrological responses of slopes were analysed. Secondly, the change of stress state of slope soil and modelling methods of slope failure under rainfall infiltration were reviewed. (2) The characteristics of the Yangbaodi landslide and associated rainfall triggering the failure were presented. The failure was characterized by shallow flowslide, due to an increase of ground water table caused by rainfall infiltration. (3) A fully automated instrumentation was carried out to monitor rainfall, and saturated – unsaturated hydrological response of the fill slope, using a raingauge, piezometers, tensiometers and moisture probes. A conceptual hydrogeological model was presented based on field monitoring and borehole data. Analysis of monitoring data showed that the high pore water pressure in fill slope was caused by upward flow of semiconfined groundwater in the moderately decomposed granite. (4) Laboratory and in-situ testing was performed to study the physical and mechanical properties of fills. Isotropically consolidated undrained compression tests and anisotropically consolidated constant shear stress tests were carried out to understand the failure mechanism of the fill slope. It is indicated that loosely compacted soil is of strain-softening behaviour under undrained conditions, accompanied with a rapid increase in excess pore water pressure. In anisotropically consolidated constant shear stress tests, a very small axial strain was required to induce the failure and the excess pore water pressure increased quickly at failure. This indicated that static liquefaction caused by rise in groundwater table due to rainfall infiltration occurred. (5) The hydraulic conductivity of the highly and moderately decomposed granite was estimated using monitering data of pore water pressure. A saturated – unsaturated flow was modeled to study the hydrological response of the fill slope using rainfall records. It was observed that the lagged failure was due to the geological conditions and the discrepancy of hydraulic conductivity of slope soils. The hydraulic conductivity of moderately decomposed granite is relatively higher than the other materials, resulting in a semiconfied groundwater flow in the moderately decomposed granite, and subsequent upward flow into the upper fill layer. When the ground water table in the fill layer was increased to the critical state, the fill slope failed. (6) Numerical exercises were conducted to replay the failure process of the fill slope, based on field monitoring, laboratory and in-situ testing. It was found that the fill slope was mobilized by a rapid transfer of the concentrated shear stress. The movement of failure mass was characterized by viscosity fluid with a gradual increase in velocity. The failure process, including mobilization and subsequent movement and deposition, was studied using numerical methods.

滑坡型泥石流转化机理研究—以四川省达县青宁乡滑坡为例

Slide-debris flow is debris flow which is transformed from landslide consecutively in a short time, it comprises of two phases: First, Landslide starts to slide; Second, landslide changes to debris flow. Slide-debris flow which brings great property and life loss happens frequently at home and abroad. In order to forecast the happening possibility and scope of slide-debris flow, transfromation mechanism of Slide-debris flow must be studied. Research on transformation mechanism of slide-debris flow is intersectant science of landslide kinetics and debris flow starting theory, It is a fringe problem as well as front problem of geological hazard. This paper takes Qingning slide-debris flow in Da County, Sichuan Province for example and has studied the mechanism of its instability and transfromation into debris flow through indoor test (including usual soil test and ring shear test) and digital modeling method.The research gets the following conclusions. Qingning Landslide took place mainly because of confined water head arising from rainfall infiltration. Before Landslide occurring, it rained continuously for 22 days, accumulated precipitation arrived at 521.6mm.Investigation shows that strata of Qingning Landslide contains quaternary loose accumulation, slip soil and highly weathered bedrock, which is a good condition for formation of confined water in the slope. Further more, groundwater seepage in the slope body and corresponding slope safety factor before landslide occurring have been computed through finite element method. The result shows that because of infiltration of rainfall, confined water head in the slope arose sharply, accordingly, the safety factor of the slope declined quickly. The result also shows that force put on the slide body by the rock mass detached from Dazhaiyan mountain was the direct factor for landslide occurring. Qingning slide-debris transformation mode has been summarized, the process the landslide changed into debris flow is divided into three phases in the prospective of macroscopic geological condition: landslide occurring, transformation and debris flow. Landslide occurring phase is from slope’ local creeping slide to Landslide occurring; transformation phase contains slide body sliding on the slide bed after slide occurring and sliding on the slope after shearing opening; debris flow phase is that slide body breaks up completely and flows downward into the ditches. The transformation mechanism of Qingning slide-debris flow has been studied through indoor ring shear test of slip soil. The result shows that transformation mechanism contains two points: first, during slide body sliding on the slide bed and slope after shearing opening, shearing shrinkage, grain crushing and grain layering brought about declining of its volume and produced excess pore water pressure, and because producing velocity of excess pore water pressure is much greater than its dissipating velocity, shear strength of slide body decreased sharply because of accumulated pore water pressure. Second, grains crushing and grains layering during slide body sliding brought about thick liquefied layer at the bottom of the slidebody, liquefied layer contained high water content and its shear strength was very low, its thickness increased as the sliding displacement increasing. Liquefied layer makes slide body sliding fast and easily break down to debris flow. Excess pore water pressure and liquefied layer made shear strength of slidebody became very low, furthermore, water in the pit of slope joining in the slidebody was also a facter that made slidebody accelerate the transformation. Influence of slide body thickness and fine grains content to transformation of slide-debris flow has been studied through ring shear test. The result reaches two conclusions. First, thickness of slide body affects transformation of slide-debris flow by two ways, porewater pressure and effect of “soft base” increases as thickness of slide body increasing.so the thicker slide body is ,the easier transformation is. Second, actual dissipating velocity of porewater pressure should be considered when studying the influence of fine grains content to tranformation of slide-debris flow. There should be a critical content of fine grains which makes the difference of producing and dissipating velocity of water pore pressre greatest, this value is the best for slide-debris transformation. The whole process of slide-debris flow transformation is reproduced through discrete element method. Transformation mechanism of slide-debris flow is studied through monitoring various parameters including pore water pressure, grain crushing and grain layering in the slide body during the transformation. The result confirms and supplements the transformation mechanism of slide-debris flow got from ring shear test well.

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

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.

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Isotope Geochemistry of Gold Ore Deposits in the Gezhen Shear Zone, Qiongxi, Hainan Island

Gold deposits hosted in the Gezhen shear zone at Qingxi, Hainan Island occur in the Precambrian metamorphic rock series and are regionally developed in the N-E direction along the tectonic zone. From northeast to southwest are distributed the Tuwaishan-Baoban gold mining district, the Erjia gold mining district and the Bumo gold mining district, making up the most industrially important gold metallogenesis zone on the Hainan Island. Isotope geochemical studies of the typical gold deposits in this metallogen

Experimental velocity measurements and effect of flow maldistribution on predicted permeator performances

The shell-side flow distribution of a parallel flow hollow-fiber gas permeator was characterized using a thermo-anemograph. The permeator has an internal diameter of 103 mm and contains 21000 fibers. The overall fiber packing fraction is 40.1%. Experimental results revealed that shell-side flow maldistributions exist in the operating conditions studied. The gas-flow velocity is the highest at the permeator center, but lowest near the shell wall. The effects of shell-side flow maldistribution on predicted permeator performances are discussed with a simple model. Model calculation results show that flow maldistributions can have considerable effect on permeation systems with relatively high separation factors and stage cuts. (C) 1998 Elsevier Science B.V.