牡丹江碰撞带的组成及演化

Jiamusi Massif is an important tectonic unit in Northeast China. It’s significant for understanding the evolution of Paleo-Asian Ocean and reconstruction of the tectonic framework of Northeast China. Mudanjiang area is located in the southern margin of Jiamusi Massif and is the key to understand the evolution of Jiamusi Massif. However, the detailed geological research for Mudanjiang area has long been deficient in many important problems, such as the tectonic components of the Mudanjiang collision zone (MCZ), the age of collisional complexes and the scenario of tectonic evolution. Based on the lithology, geochemistry and the SHRIMP zircon U-Pb geochronology in Mudanjiang area, our new data and results come to some constraints for the tectonic reconstruction of MCZ as follows: 1） It is identified that the former suggestion, which the so-called “Heilongjiang Group” in Mudanjiang area is the vestige of oceanic crust, is correct. The oceanic relics represent the Neo-Proterozoic-Early Paleozoic oceanic basins based on the SHRIMP zircon U-Pb geochronology. 2） One sheet of gabbroic complex with oceanic island-type geochemical signature was discovered by this work in Mudanjiang area. 3） It is verified that the Proterozoic concordant U-Pb ages of the migmatites developed along the southern margin of Jiamusi massif, which represent the events of magmatic intrusion, as the direct evidence for the existence of the Proterozoic crystalline basements of the Jiamusi Massif. Based on geochronology, we suggest that the migmatization and coeval S-type granite magmatism of the southern margin of Jiamusi Massif took place about 490Ma. 4） The island arc complex has been found in the Heilongjiang Group, and the oceanic relics was found distributing on both sides, as provided important constraint for the tectonic reconstruction of the MCZ. 5） ~440Ma metamorphic event and coeval post-collisional granite magmatism have been firmly identified in the MCZ and its southern neighboring area. Together with previous data obtained by other researchers, our conclusions on the reconstruction of the tectonic architecture and evolution of the MCZ as follows: 1） The orogenic assemblages developed in the Mudanjiang collisional zone are featured by a sequence of ancient active continental margins and ensuing orogenic processing. The Mashan Group is the reworking basement of Jiamusi Massif, whereas the Heilongjiang Group represents arc and oceanic complexes characterized by imbricate deep-seated sliced and slivering sheets due to multi-phases of thrusting and nappe stacking. 2） The northern sub-belt of MCZ is probably the arc-continent collisional boundary related to the closure of main oceanic basin. The collisional age can be constrained by the events of syn-orogenic migmatization of migmatite, coeval S-type granite magmatism and the related granulite-facies metamorphism. Therefore, we suggested the collisional age of northern sub-belt is probably Cambrian-Early Ordovician. The extensive granulite-facies metamorphism of the Mashan Group in Jiamusi Massif, as affirmed by former works, was probably related with the collisional event. 3） The southern sub-belt of the MCZ was possibly related with the closure of back-arc basin. We presumed that the collisional age of southern sub-belt is about Ordovician-Early Silurian according to the ~440Ma extensive metamorphism and the occurrence of coeval post-collisional granite magmatism. 4） The extant structural architecture of the MCZ is related to the multi-phases of intra-continental superimposition, which is characterized by the Mesozoic nappe structure.

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

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.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

预应力结构抗滑桩设计计算方法及其在滑坡治理中的应用研究

Landslide is a kind of serious geological hazards and its damage is very great. In recent years, landslides become more and more frequent along with increase of scale of engineering constructions and cause greater loss. Consequently, how to protect landslides has become important research subject in the engineering field. This paper improves the method how to compute landslide thrust and solves the irrational problem in the design of piles because of the irrational landslide thrust according to the theory and technology of existed anti-slide piles and pre-stressed cable anti-slide piles. Modern pre-stressing technology has been introduced and load balancing method has been used to improve the stressing behavior of anti-slide piles. Anchor cables, anti-slide piles and modern pre-stressing technology have been used to prevention complicated landslide. It is an important base to select values for the landslide thrust. An improved method to calculate design thrust of anti-slide piles has been presented in this paper on the base of residual thrust method by comparing existing methods to select values of landslide thrust in the design of anti-slide piles. In the method, residual landslide thrust behind the anti-slide piles and residual skid resistance before the piles has been analyzed, equitable distribution of residual landslide thrust behind the piles has been realized, and the method to select value of design thrust becomes more reasonable. The pre-stressed cable anti-slide piles are developed from the common anti-slide piles and are common method to prevent landslide. Their principle is that internal force of anti-slide piles is adjusted and size of section is diminished by changing constraint conditions of anti-slide piles. For landslides with deep slip surface and large scale of slopes, limitation of the method appears. Such landslides are in need of long piles and anchor cables which are not only non-economic but also can generate larger deformation and leave potential danger after prevention. For solving the problem, a new kind of anti-slide piles, inner pre-stressing force anti-slide piles, is presented in this paper, and its principle is that an additional force, which is generated in the inner anti-slide piles by arranging pre-stressed reinforcement or tight wire in a certain form in interior of anti-slide piles and stretching the steel reinforcement or tight wire, may balance out the internal force induced by landslide thrust whole or partly (load balancing method). The method will change bending moment which anti-slide piles are not good at bearing into compressive stress which piles are good at bearing, improve stressing performance of anti-slide piles greatly, diminish size of section, and make anti-slide piles not fissured in the natural service or postpone appearance of the fissures, and improve viability of anti-slide piles. Pre-stressed cable anti-slide piles and inner pre-stressing force anti-slide piles go by the general name of pre-stressed structure anti-slide piles in the paper, and their design and calculation method is also analyzed. A new calculation method is provided in the paper for design of anti-slide piles. For pre-stressed structure anti-slide piles, a new computation mode is firstly presented in the paper on the foundation of cantilever piles. In the mode, constraint form of load-bearing section of the anti-slide piles should be confirmed according to reservoir conditions in order to figure out amount of pre-stress of the anchor cables, and internal force should be analyzed for the load-bearing section of pre-stressed structure anti-slide piles so as to confirm anchorage section of anti-slide piles. Pre-stressed cables of the pre-stressed cable anti-slide piles can be arranged as required. This paper analyzes the load-bearing section of single-row and double-row pre-stressed cable anti-slide piles and provides a calculation method for design of the pre-stressed cable anti-slide piles. Inner pre-stressing force anti-slide piles are a new kind of structural style. Their load-bearing section is divided into four computation modes according to whether pre-stressed cables are applied for exterior of the anti-slide piles, and whether single-row or double-row exterior pre-stressed cables are applied. The load balancing method is used to analyze the computation modes for providing a method to design the inner pre-stressing force anti-slide piles rationally. Pre-stressed cable anti-slide piles and inner pre-stressing force anti-slide piles are applied to research on Mahe landfall in Yalong Lenggu hydropower station by the improved method to select value of design thrust of anti-slide piles. A good effect is obtained in the analysis.

笼形水簇和溶解甲烷之间平均力势能的约束型分子动力学计算及其在天然气水合物成核研究中的意义

That the dodecahedral water cluster (DWC) can adsorb dissolved methane molecules, an important phenomenon related to the hydrate nucleation study, has been observed through molecular dynamics simulations, but it has not been explained satisfactorily [Guang-Jun Guo; Yi-Gang Zhang; Hua Liu. J. Phys. Chem. C, 2007, 111, 2595]. In order to explain this phenomenon by using the potential of mean force (PMF) between the DWC and the dissolved methane, we perform several series of constrained molecular dynamics simulations in the methane-water system. The distance between the center of DWC and the methane molecule is constrained from 5 Å to 18 Å by adding 0.2 Å every time. For each fixed distance, we perform 20 independent simulations to improve the statistical precision. We first get the constraint force between the DWC and the dissolved methane in each simulation and then calculate the PMF by integrating these forces. Subsequently, the radial distribution function (RDF) is obtained from the PMF through an equation of statistical mechanics. The results show that the RDF has a sharp peak at about 6.2 Å, successfully explaining why the DWC adsorbs dissolved methane molecules. The preferential binding coefficient is a positive value (=2.05±0.5), indicates that the DWC tends to adsorb dissolved methane rather than water molecules in methane aqueous solutions. The curve of PMF for the DWC encaging a methane almost coincides that for the empty DWC, meaning that it is the DWC rather than the encaged methane who could adsorb dissolved methane molecules. By comparing the curves of PMF for different directions of the DWC relative to the dissolved methane, we find that it is the cage face rather than the cage edge or vertex that plays an essential role when the DWC adsorbing dissolved methane. This research sheds light on the driving force for the methane adsorption, and it is helpful in understanding the nucleation process of methane hydrate.

地震资料叠前储层预测综合研究--以辽中凹陷东二下段浊积扇气藏为例

Reservoir prediction techniques from prestack seismic are among the most important ones for exploration of lithologic hydrocarbon reservoir. In this paper, we set the turbidite fan sandstone reservoir in Liao-Zhong depress as our researching target, and aims to solve the apllication difficulties on pre-stack inversion in the area, where the drilling data is scarce and the reservoir is lateral varied. Meanwhile, AVO analysis and pre-stack inversion for gas-bearing detection is systematically researched. The seismic reflection characters of gas-bearing sandstone in turbidite fan with different fluid content are defined, after analyzing results from AVO seismic simulation and porous fluid replacement of real log data, and under the guides of the seismic characters from classical gas-bearing sandstone reservoir and numerical simulation for complicate gas-bearing sandstone. It is confirmed that detecting gas-bearing sandstone in turbidite fan via AVO technologies is feasible. In terms of AVO analysis, two AVO characters, fluid detection factor and product of intercept and gradient, can effectively identify top and bottom boundaries and lateral range of tuibidite gas sand by comparing real drilling data. Cross-plotting of near and far angle stack data could avoid the correlation existing in P-G analysis. After comparing the acoustic impedance inversions with routine stacked data and AVO intercept, impedance derived from AVO intercept attribute could reduce the acoustic impedance estimating error which is caused by AVO. On the aspect of elastic impedance inversion, the AVO information in the pre-stack gathers is properly reserved by creating partial angle stack data. By the far angle elastic impedance alone, the gas sand, with abnormally low range of values, can be identified from the background rocks. The boundary of gas sand can also be clearly determined by cross-plotting of near and far angle elastic impedances. The accuracy of far angle elastic impedance is very sensitive to the parameter K, and by taking the statistical average of Vp/Vs on the targeted section in key wells, the accuracy of low frequency trends is gurranteed; the intensive absorsion within the area of the gas sand, which tends to push the spectral of seismic data to the lower end, will cause errors on the inversion result of elastic impedance. The solution is to confine the inversion on the interested area by improving the wavelet. On the aspect of prestack AVA simultaneous inversion, the constraint of local rock-physical trends between velocities of P-wave、S-wave and density successfully removes the instability of inversion, thus improves the precision of the resulting elastic parameters. Plenty of data on rock properties are derived via AVO analysis and prestack seismic data inversion. Based on them, the fluid anomaly is analysized and lithological interpretation are conducted. The distribution of gas sand can be consistently determined via various of ways, such as cross-plotting of P and G attributes, near and far partial angle stack data, near and far angle elastic impedances, λρ and Vp/Vs, etc. The shear modulo and density are also reliable enough to be used for lithological interpretation. We successfully applied the AVO analysis and pre-stack inversion techniques to gas detecting for turbidite fan sand reservoir in Liao-Zhong depression.

叠前地震多角度三参数同步反演方法

Pre-stack seismic inversion has become the emphasis and hotspot owing to the exploration & exploitation of oil field and the development of seismic technology. Pre-stack seismic inversion has the strongpoint of making the most of amplitude versus offset compared with the post-stack method. In this dissertation, the three parameters were discussed from multi-angle reflectance of P-wave data based on Zoeppritz’s and Aki & Richard’s equation, include P-wave velocity, S-wave velocity, and density. The three parameters are inversed synchronously from the pre-stack multi-angle P-wave data, based on rockphysics model and aimed at the least remnant difference between model simulation and practical data. In order to improve the stability of inversion and resolution to thin bed, several techniques were employed, such as the wavelet transform with multi-scale function, adding the Bayesian soft constraint and hard constraints (the horizon, structure and so on) to the inversion process. Being the result, the uncertainty of the resolution is reduced, the reliability and precision are improved, the significance of parameters becomes clearer. Meeting to the fundamental requirement of pre-stack inversion, some research in rockphysics are carried out which covered the simulation and inversion of S-wave velocity, the influence of pore fluids to geophysical parameters, and the slecting and analyzing of sensitive parameters. The difference between elastic wave equation modeling and Zoeppritz equation method is also compared. A series of key techniques of pre-stack seismic inversion and description were developed, such as attributes optimization, fluid factors, etc. All the techniques mentioned above are assembled to form a technique sets and process of synchronous pre-stack seismic inversion method of the three parameters based on rock physics and model simulation. The new method and technology were applied in many areas with various reservoirs, obtained both geological and economic significance, which proved to be valid and rational. This study will promote the pre-stack inversion technology and it’s application in hidden reservoirs exploration, face good prospects for development and application.

胶东地区金矿床成矿时代及其成矿地球动力学背景-兼论壳幔相互作用与成岩成矿

Jiaodong Peninsula is the largest repository of gold in China. Varieties of studies have been involved in the mechanism of metallogenesis. This thesis is a part of the project "Study of basic geology related to the prespecting of the supra-large deposits" which supported by National Climbing Program of China to Prof. Zhou. One of the key scientific problems is to study the age and metallogenic dynamics of ore deposit and to understand how interaction between mantle and crust constrains on metallogenesis and lithogenesis. As Jiaodong Peninsula to be study area, the Rb-Sr, Sm-Nd and Pb isotopic systematics of pyrite and altered rocks are measured to define the age and origin of gold. The elemental and Sr-Nd-Pb isotopic compositions of dikes and granites was studied to implicate the source and lithogenesis of the dike and granite and removal of lithosphere and the interaction between mantle and crust in the Jiaodong Peninsula. Considering the tectonic of Jiaodong Peninsula, basic on the time and space, this thesis gives a metallogenic dynamics of gold mineralization and discusses the constraints of the interaction between mantle and crust on the metallogenesis and lithogenesis. This thesis reports the first direct Rb-Sr dating of pyrites and ores using sub-sampling from lode gold deposit in Linglong, Jiaodong Peninsula and the results demonstrate this as a useful geochronological technique for gold mineralization with poor age constraint. The Rb-Sr data of pyrites yields an isochron age of (121.6-122.7) Ma, whereas, those of ore and ore-pyrite spread in two ranges from 120.0 to 121.8 Ma and 110.0-111.7 Ma. Studies of characteristic of gold deposit, microscopy of pyrite and quartz indicate that the apparent ages of ore and ore-pyrite are not isochron ages, it was only mixed by two end members, i.e., the primitive hydrothermal fluids and wall rocks. However, the isochron age of pyrite samples constrains the age of gold mineralization, i.e., early Cretaceous, which is in good consistence with the published U-Pb ages of zircon by using the SHRIMP technique. The whole rock Rb-Sr isochron age of altered rocks indicates that the age of gold mineralizing in the Xincheng gold deposit is 116.6 ± 5.3 Ma. The Sr, Nd and Pb isotopic compositions of pyrite and altered rocks indicate that the gold and relevant elements were derived from multi-sources, i.e. dikes derived from enriched lithospheric mantle and granites, granodiorites and metamorphic rocks outcropped on the crust. It also shows that the hydrothermal fluids derived from mantle magma degassing had play an important role in the gold mineralizing. The major and trace elements, Sr-Nd-Pb isotopic data of granites and granodiorites suggest that the Linglong Granite and Kunyushan Granite were derived from partial melting of basement rocks in the Jiaodong Peninsula at post-collision of North China Craton with South China Craton. Guojialing Granodiorite was considered to be derived from a mixture source, that is, mixed by magmas derived from an enriched lithospheric mantle and crust during the delamination of lithosphere induced by the subduction of Izanagi Plate and the movement of Tancheng-Lujiang Fault. There are kinds of dikes occurred in the Jiaodong Peninsula, which are accompanying with gold mineralization in time and space. The dikes include gabrro, diabase, pyroxene diorite, gabrrophyre, granite-porphyry, and aplite. The whole rock K-Ar ages give two age intervals: 120-124 Ma for the dikes that erupted at the gold mineralizing stage, and <120 Ma of the dikes that intruded after gold mineralizing. According to the age and the relationship between the dikes and gold mineralizing, the dikes could be divided into two groups: Group I (t = 120-124 Ma) and Group II (t < 120Ma). Group I dikes show the high Mg and K, low Ti contents, negative Nb anomalies and positive Eu anomalies, high ~(87)Sr/~(86)Sr and negative εNd(t) values and an enrichment in light rare earth elements, large ion lithosphile elements and a depletion in high field strength elements. Thus the elemental and isotopic characteristics of the Group I dikes indicate that they were derived from an enriched lithospheric mantle perhaps formed by metasomatism of the melt derived from the recycled crustal materials during the deep subduction of continent. In contrast, the Group II dikes have high Ti, Mg and K contents, no negative Nb anomalies, high ~(87)Sr/~(86)Sr and positive or little negative εNd(t) values, which indicate the derivation from a source like OIB-source. The geochemical features also give the tectonic constraints of dikes, which show that Group I dikes were formed at continental arc setting, whereas Group II dikes were formed within plate background. Considering the tectonic setting of Jiaodong Peninsula during the period of gold mineralizing, the metallogenic dynamics was related to the subduction of Izanagi Plate, movement of Tancheng-Lujiang Fault and removal of lithopheric mantle during Late Mesozoic Era.

断陷湖盆缓坡带成藏模式研究

Halfgraben-like depressions have multiple layers of subtle traps, multiple coverings of oil-bearing series and multiple types of reservoirs. But these reservoirs have features of strong concealment and are difficult to explore. For this reason, many scholars contribute efforts to study the pool-forming mechanism for this kind of basins, and establish the basis for reservoir exploration and development. However, further study is needed. This paper takes HuiMin depression as an example to study the pool-forming model for the gentle slope belts of fault-depression lake basins. Applying multi-discipline theory, methods and technologies including sedimentary geology, structural geology, log geology, seismic geology, rock mechanics and fluid mechanics, and furthermore applying the dynamo-static data of oil reservoir and computer means in maximum limitation, this paper, qualitatively and quantitatively studies the depositional system, structural framework, structural evolution, structural lithofacies and tectonic stress field, as well as fluid potential field, sealing and opening properties of controlling-oil faults and reservoir prediction, finally presents a pool-forming model, and develops a series of methods and technologies suited to the reservoir prediction of the gentle slope belt. The results obtained in this paper richen the pool-forming theory of a complex oil-gas accumulative area in the gentle slope belt of a continental fault-depression basin. The research work begins with the study of geometric shape of fracture system, then the structural form, activity stages and time-space juxtaposition of faults with different level and different quality are investigated. On the basis of study of the burial history, subsidence history and structural evolution history, this paper synthesizes the studied results of deposition system, analyses the structural lithofacies of the gentle slope belt in the HuiMing Depression and its controlling roles to oil reservoir in the different structural lithofacies belts in time-space, and presents their evolution patterns. The study of structural stress field and fluid potential field indicates that the stress field has a great change from the Dong Ying stages to nowadays. One marked point among them is that the Dong Ying double peak- shaped nose structures usually were the favorable directional area for oil and gas migration, while the QuDi horst became favorable directional area since the GuanTao stage. Based on the active regular of fractures and the information of crude oil saturation pressure, this paper firstly demonstrates that the pool-forming stages of the LingNan field were prior to the stages of the QuDi field, whici provides new eyereach and thinking for hydrocarbon exploration in the gentle slope belt. The BeiQiao-RenFeng buried hill belt is a high value area with the maximum stress values from beginning to end, thus it is a favorable directional area for oil and gas migration. The opening and sealing properties of fractures are studied. The results obtained demonstrate their difference in the hydrocarbon pool formation. The seal abilities relate not only with the quality, direction and scale of normal stress, with the interface between the rocks of two sides of a fault and with the shale smear factor (SSF), but they relate also with the juxtaposition of fault motion stage and hydrocarbon migration. In the HuiMin gentle slope belt, the fault seal has difference both in different stages, and in different location and depth in the same stage. The seal extent also displays much difference. Therefore, the fault seal has time-space difference. On the basis of study of fault seal history, together with the obtained achievement of structural stress field and fluid potential field, it is discovered that for the pool-forming process of oil and gas in the studied area the fault seal of nowadays is better than that of the Ed and Ng stages, it plays an important role to determine the oil column height and hydrocarbon preservation. However, the fault seal of the Ed and Ng stages has an important influence for the distribution state of oil and gas. Because the influential parameters are complicated and undefined, we adopt SSF in the research work. It well reflects synthetic effect of each parameter which influences fault seal. On the basis of the above studies, three systems of hydrocarbon migration and accumulation, as well as a pool-forming model are established for the gentle slope belt of the HuiMin depression, which can be applied for the prediction of regular patterns of oil-gas migration. Under guidance of the pool-forming geological model for the HuiMin slope belt, and taking seismic facies technology, log constraint evolution technology, pattern recognition of multiple parameter reservoir and discrimination technology of oil-bearing ability, this paper develops a set of methods and technologies suited to oil reservoir prediction of the gentle slope belt. Good economic benefit has been obtained.

地震波散射非线性反演问题不动点理论及方法的研究

How to create a new method to solve the problem or reduce the influence of that the result of the seismic waves scattering nonlinear inversion is not uniqueness is a main purpose of this research work in the paper. On the background of research into the seismic inversion, new progress of the nonlinear inversion is introduced at the first chapter in this paper. Especially, the development, basic theories and assumptions on some major theories of seismic inversion are analyzed, discussed and summarized in mathematics and physics. Also, the problems faced by the mathematical basis of investigations of the seismic inversion are discussed, and inverse questions of strongly seismic scattering due to strong heterogeneous media in the Earth interior are analyzed and viewed. What the kernel of paper is that gathers all our attention making a new nonlinear inversion method of seismic scattering. The paper provides a theory and method of how to introduce the fixed-point theory into the nonlinear seismic scattering inversion and how to obtain the solution, and gives the actually method to create a serials of contractive mappings of velocity parameter's in the mapping space of wave. Therefore, the results testify the existence of fixed point of velocity parameter and give the method the find it. Further, the paper proves the conclusion that the value obtained by taking the fixed point of velocity parameter into wave equation is the fixed point of the wave of the contractive mapping. Thence, the fixed point is the global minima since the stabilities quality of the fixed point. Based on the new theory, in the chapter three, many inverse results are obtained in the numerical value test. By analysis the results one could find a basic facts that all the results, which are inversed by the different initial model, are tended to the true value in theoretical true model. In other words, the new method mostly eliminates the non-uniqueness that which is existed in seismic waves scattering nonlinear inversion in degree. But, since the test results are quite finite now, more test is need here to positive our theory. As a new theoretical method, it must be existed many weaken in it. The chapter four points out all the questions which is bother us. We hope more people to join us to solve the problem together.

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

利用分频与反演提高地震资料信噪比与分辨率方法研究

This thesis mainly talks about the wavelet transfrom and the frequency division method. It describes the frequency division processing on prestack or post-stack seismic data and application of inversion noise attenuation, frequency division residual static correction and high resolution data in reservoir inversion. This thesis not only describes the frequency division and inversion in theory, but also proves it by model calculation. All the methods are integrated together. The actual data processing demonstrates the applying results. This thesis analyzes the differences and limitation between t-x prediction filter and f-x prediction filter noise attenuation from wavelet transform theory. It considers that we can do the frequency division attenuation process of noise and signal by wavelet frequency division theory according to the differences of noise and signal in phase, amplitude and frequency. By comparison with the f-x coherence noise, removal method, it approves the effects and practicability of frequency division in coherence and random noise isolation. In order to solve the side effects in non-noise area, we: take the area constraint method and only apply the frequency division processing in the noise area. So it can solve the problem of low frequency loss in non-noise area. The residual moveout differences in seismic data processing have a great effect on stack image and resolutions. Different frequency components have different residual moveout differences. The frequency division residual static correction realizes the frequency division and the calculation of residual correction magnitude. It also solves the problems of different residual correction magnitude in different frequency and protects the high frequency information in data. By actual data processing, we can get good results in phase residual moveout differences elimination of pre-stack data, stack image quality and improvement of data resolution. This thesis analyses the characters of the random noises and its descriptions in time domain and frequency domain. Furthermore it gives the inversion prediction solution methods and realizes the frequency division inversion attenuation of the random noise. By the analysis of results of the actual data processing, we show that the noise removed by inversion has its own advantages. By analyzing parameter's about resolution and technology of high resolution data processing, this thesis describes the relations between frequency domain and resolution, parameters about resolution and methods to increase resolution. It also gives the processing flows of the high resolution data; the effect and influence of reservoir inversion caused by high resolution data. Finally it proves the accuracy and precision of the reservoir inversion results. The research results of this thesis reveal that frequency division noise attenuation, frequency residual correction and inversion noise attenuation are effective methods to increase the SNR and resolution of seismic data.

Analysis of Differential and Matching Methods for Optical Flow

Several algorithms for optical flow are studied theoretically and experimentally. Differential and matching methods are examined; these two methods have differing domains of application- differential methods are best when displacements in the image are small (<2 pixels) while matching methods work well for moderate displacements but do not handle sub-pixel motions. Both types of optical flow algorithm can use either local or global constraints, such as spatial smoothness. Local matching and differential techniques and global differential techniques will be examined. Most algorithms for optical flow utilize weak assumptions on the local variation of the flow and on the variation of image brightness. Strengthening these assumptions improves the flow computation. The computational consequence of this is a need for larger spatial and temporal support. Global differential approaches can be extended to local (patchwise) differential methods and local differential methods using higher derivatives. Using larger support is valid when constraint on the local shape of the flow are satisfied. We show that a simple constraint on the local shape of the optical flow, that there is slow spatial variation in the image plane, is often satisfied. We show how local differential methods imply the constraints for related methods using higher derivatives. Experiments show the behavior of these optical flow methods on velocity fields which so not obey the assumptions. Implementation of these methods highlights the importance of numerical differentiation. Numerical approximation of derivatives require care, in two respects: first, it is important that the temporal and spatial derivatives be matched, because of the significant scale differences in space and time, and, second, the derivative estimates improve with larger support.

Olympic Robot Building Manual

The 1989 AI Lab Winter Olympics will take a slightly different twist from previous Olympiads. Although there will still be a dozen or so athletic competitions, the annual talent show finale will now be a display not of human talent, but of robot talent. Spurred on by the question, "Why aren't there more robots running around the AI Lab?", Olympic Robot Building is an attempt to teach everyone how to build a robot and get them started. Robot kits will be given out the last week of classes before the Christmas break and teams have until the Robot Talent Show, January 27th, to build a machine that intelligently connects perception to action. There is no constraint on what can be built; participants are free to pick their own problems and solution implementations. As Olympic Robot Building is purposefully a talent show, there is no particular obstacle course to be traversed or specific feat to be demonstrated. The hope is that this format will promote creativity, freedom and imagination. This manual provides a guide to overcoming all the practical problems in building things. What follows are tutorials on the components supplied in the kits: a microprocessor circuit "brain", a variety of sensors and motors, a mechanical building block system, a complete software development environment, some example robots and a few tips on debugging and prototyping. Parts given out in the kits can be used, ignored or supplemented, as the kits are designed primarily to overcome the intertia of getting started. If all goes well, then come February, there should be all kinds of new members running around the AI Lab!