渤海湾盆地济阳坳陷东部走滑构造特征及其对油气成藏的影响研究

研究区位于郯庐断裂中段与济阳坳陷的构造结合部，区内走滑构造广泛发育，主要的走滑断裂有7条，分别是郯庐断裂带的东西两支、垦东断层、孤东断层、长堤断层、埕东断层和发育于垦东凸起中部的浅层走滑构造带。走滑构造带与油气富集带有着明显的对应关系。 通过对研究区内二维、三维地震测线和平面构造图的精细解释和分析，分别揭示了各走滑断裂在平面、剖面和三维空间上的构造形态。根据走滑断裂及其伴生构造的平面和剖面上的几何学特征，将研究区内的走滑断裂划分为三种类型：成熟型走滑断裂、隐伏型走滑断裂、不连续型的走滑断裂。 从理论模式研究入手，推导了拉分盆地中盆地的走滑速率与沉降速率之间的关系，证实了走滑速率同盆地的几何形状参数、最大沉降深度和盆地的沉降速率存在着稳定的数值关系。通过对莱州湾地区潍北凹陷基底沉降历史的分析，建立了潍北凹陷沉降速率与郯庐断裂中段走滑速率之间的经验关系式，进而求出郯庐断裂中段新生代右行走滑位移量的大小为40km。 运用2DMove软件，对研究区内四条典型剖面进行构造复原，计算出了各条剖面每个时期的伸展参数，对研究区构造活动强度进行了定量分析，揭示了研究区的构造演化规律。通过运用Ansys软件进行有限元模拟，恢复了晚白垩世晚期-古近纪早期研究区内的构造应力场和应变场，揭示了扭张作用是研究区内走滑断层开始走滑的主要原因。 通过上述分析，结合对究区内近几年勘探开发成功和失败的实例分析，全面探讨了走滑活动对于油气成藏“生”、“储”、“盖”、“圈”、“运”、“保”各因素的影响。

新疆东昆仑蛇绿岩岩石学、地球化学特征及其大地构造意义

The East Kunlun area of Xinjiang (briefly EKAX) is the western part of broadly speaking East Kunlun orogenic zone. The absence of geological data (especially ophiolites) on this area has constrained our recognition to its geology since many years. Fund by National 305 Item (96-915-06-03), this paper, by choosing the two ophiolite zones (Muztag and Southwestern Margin of Aqikekule Lake ophiolite zones) exposed at EKAX as the studied objects and by the analysis of thin section, electron probe, XRF, ICP-MS, SEM and Sm-Nd isotope, totally and sys ematically dealt with the field geological, petrological, minerological, petrochemical and geochemical characteristics (including trace, rare earth element and Sm-Nd isotope) and the tectonic setting indicated by them for each ophilite zone. Especially, this paper discussed the trace and rare earth element patterns for metamorphic peridotites, their implications and related them to the other components of ophiolite in order to totally disclose ophiolite origins. Besides, this paper also studied the petrological, geochemical and paleobiological characteristics for the cherts coexsisted with the Muztag ophiolite and the tectonic setting indicated by them. Based on these, the author discussed the tectonic evolution from Proterozoic to Permian for this area. For Muztag ophiolite, their field geological, petrological, minerological, petrochemical and geochemical characteristics show that: ① outcropped along the Muztag-Jingyuhu fault with west-to-east strike, the ophiolite is composed of such three components as metamorphic peridotites, cumulates and volcanic rocks; ② metamophic peridotites consist of such types as lherzolites, serpentinized lherzolites and serpentinites, only pyroxenites is seen of cumulates and volcanic rocks include basalts, basaltic andesites and andesites; ③ mineralogical data on this ophiolite suggest it formed in supra-subduction zone (SSZ)environment, and its mantle wedge is heterogeneous; ④ whole-rock TiO_2 and Al_2O_3 of metamorphic peridotites indicate their original environment with the MORB and SSZ characteristics; ⑤ metamorphic peridotites have depleted LREE and flat REE patterns and volcanic rocks have enriched LREE patterns; ⑥ trace element characteristics of metamorphic peridotites imply that they had undergone Nb and Ta enrichment event after partial melting; ⑦ trace element characteristics of volcanic rocks and their tectonic diagrams show they are formed in the spreading and developed island arc environment with back-arc basin, such as rifted island arc, which is supported by the ε_(Nd)(t) -2.11～+3.44. In summary, the above evidence implies that Muztag ophiolite is formed in SSZ environment, where heterogeneous mantle wedge was metasomatised by the silica-enriched melt from subducted sediments and/or oceanic crust, which makes the mantle wedge enriched again, and this enriched mantle wedge later partially melted to form the volcanic rocks. For Southwestern Margin of Aqikekule Lake ophiolite, their field geological, petrological, minerological, petrochemical and geochemical characteristics show that: ① it outcropped as tectonic slices along the near west-to-east strike Kunzhong fault and is composed of metamorphic perodotties, cumulates and volcanic rocks, in which, chromites are distributed in the upper part of metamorphic peridotites as pods, or in the lower part of cumulates as near-strata; ② metamorphic peridotites include serpentinites, chromite-bearing serpentinites, thlorite-epidote schists and chromitites, of which, chromitites have nodular and orbicular structure, and cumulates include pyroxenits, serpentinites, chromite-bearing serpentinites, chromites and metamorphically mafic rocks and only basalts are seen in volcanic rocks; ③ Cr# of chromites suggest that they formed in the SSZ and Al_2O_3 and TiO_2 of metamorphic peridotites also suggest SSZ environment; ④metamorphic peridotites have V type and enriched LREE patterns, cumulates have from strongly depleted LREE, flat REE to enriched LREE patterns with universally striking positive Eu anomalies and basalts show flat REE or slight enriched LREE patterns with no Eu anomalies; ⑤ trace element and Sm-Nd isotope characteristics of metamorphic peridotites imply their strikingly heterogeneous mantle character(ε_(Nd)(t)+4.39～+26.20) and later Nb, Ta fertilization; ⑥ trace element characteristics of basalts and their tectonic diagrams show they probably formed in the rifted island arc or back-arc basin enviromnent. In summary, the above evidence shows that this ophiolite formed in the SSZ environment and melts from subudcted plate are joined during its formation. Rare earth element, whole-rock and sedimentary characteristics of cherts with the Muztag ophiolite show that they formed in the continental margin environment with developed back-arc basin, and radiolarias in the cherts indicate that the upper age of Muztag ophiolite is early carboniferous. Based on the accreted wedge models of Professor Li Jiliang for Kunlunshan Mountain and combined with study on the two typical ophiolite profiles of EKAX, the author discussed the tectonic evolution of EKAX from Proterzoic to Permian.

Effects of 'sinking in' and 'piling up' on estimating the contact area under load in indentation

The phenomena of the 'piling up' and 'sinking-in' of surface profiles in conical indentation in elastic-plastic solids with work hardening are studied using dimensional and finite-element analysis. The degree of sinking in and piling up is shown to depend on the ratio of the initial yield strength Y to Young's modulus E and on the work-hardening exponent n. The widely used procedure proposed by Oliver and Pharr for estimating contact depth is then evaluated systematically. By comparing the contact depth obtained directly from finite-element calculations with that obtained from the initial unloading slope using the Oliver-Pharr procedure, the applicability of the procedure is discussed.

Relationships between precipitation, soil water and groundwater at Chongling catchment with the typical vegetation cover in the Taihang mountainous region, China

Vegetation cover plays an important role in the process of evaporation and infiltration. To explore the relationships between precipitation, soil water and groundwater in Taihang mountainous region, China, precipitation, soil water and water table were observed from 2004 to 2006, and precipitation, soil water and groundwater were sampled in 2004 and 2005 for oxygen-18 and deuterium analysis at Chongling catchment. The soil water was sampled at three sites covered by grass (Carex humilis and Carex lanceolata), acacia and arborvitae respectively. Precipitation is mainly concentrated in rainy seasons and has no significant spatial variance in study area. The stable isotopic compositions are enriched in precipitation and soil water due to the evaporation. The analysis of soil water potential and isotopic profiles shows that evaporation of soil water under arborvitae cover is weaker than under grass and acacia, while soil water evaporation under grass and acacia showed no significant difference. Both delta O-18 profiles and soil water potential dynamics reveal that the soil under acacia allows the most rapid infiltration rate, which may be related to preferential flow. In the process of infiltration after a rainstorm, antecedent water still takes up over 30% of water in the topsoil. The soil water between depths of 0-115 cm under grass has a residence time of about 20 days in the rainy season. Groundwater recharge from precipitation mainly occurs in the rainy season, especially when rainstorms or successive heavy rain events happen.

General Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters of Arbitrary Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expression is well known for indentation in elastic and in elastic-plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. Furthermore, the same expression holds true for both fast loading and unloading. These results should provide a sound basis for using the relationship for determining properties of viscoelastic solids using indentation techniques.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

Acoustic Emission Experiments of Rock Failure under Load Simulating the Hypocenter Condition

A series of acoustic emission (AE) experiments of rock failure have been conducted under cyclic load in tri-axial stress tests. To simulate the hypocenter condition the specimens are loaded by the combined action of a constant stress, intended to simulate

Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Comparison of Various Adhesion Contact Theories and the Influence of Dimensionless Load Parameter

Three models, JKR (Johnson, Kendall and Roberts), DMT (Derjaguin, Muller, and Toporov) andMD (Maugis-Dugdale),are compared with the Hertz model in dealing with nano-contact problems. It has been shown that both the dimensionless load parameter, P D P=.1/4

Fractals and Scaling in Fracture Induced by Microcrack Coalescence

A numerical model is proposed to simulate fracture induced by the coalescence of numerous microcracks, in which the condition for coalescence between two randomly nucleated microcracks is determined in terms of a load-sharing principle. The results of the simulation show that, as the number density of nucleated microcracks increases, stochastic coalescence first occurs followed by a small fluctuation, and finally a newly nucleated microcrack triggers a cascade coalescence of microcracks resulting in catastrophic failure. The fracture profiles exhibit self-affine fractal characteristics with a universal roughness exponent, but the critical damage threshold is sensitive to details of the model. The spatiotemporal distribution of nucleated microcracks in the vicinity of critical failure follows a power-law behaviour, which implies that the microcrack system may evolve to a critical state.

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

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

The instability analysis of saturated soil under shear load

A theoretical description of shear instability is presented in a system of equations. It is shown that two types of instability may exist. One of them is dominated by pore pressure softening while the other by strain softening. A criterion combining pore pressure softening, strain hardening, and volume strain coefficient is obtained and practical implications are discussed.

Scaling, Dimensional Analysis, and Indentation Measurements

We provide an overview of the basic concepts of scaling and dimensional analysis, followed by a review of some of the recent work on applying these concepts to modeling instrumented indentation measurements. Specifically, we examine conical and pyramidal indentation in elastic-plastic solids with power-law work-hardening, in power-law creep solids, and in linear viscoelastic materials. We show that the scaling approach to indentation modeling provides new insights into several basic questions in instrumented indentation, including, what information is contained in the indentation load-displacement curves? How does hardness depend on the mechanical properties and indenter geometry? What are the factors determining piling-up and sinking-in of surface profiles around indents? Can stress-strain relationships be obtained from indentation load-displacement curves? How to measure time dependent mechanical properties from indentation? How to detect or confirm indentation size effects? The scaling approach also helps organize knowledge and provides a framework for bridging micro- and macroscales. We hope that this review will accomplish two purposes: (1) introducing the basic concepts of scaling and dimensional analysis to materials scientists and engineers, and (2) providing a better understanding of instrumented indentation measurements.