Biopotentials of marine sponges from China oceans: past and future

An extensive literature survey of over 17 Journals was carried out on Chinese sponges and their natural products in the period from 1980 to 2001. This review is thus intended to provide the first thorough overview of research on marine sponges from China Ocean territories. Information is provided about the rather-limited taxonomic study of Chinese marine sponges, with an analysis on their distribution and diversity. Research findings on the natural products and their bioactivity screening from Chinese sponges are summarized. The weaknesses, gaps and problems in the past R&D program of Chinese sponges are identified, which point to the future opportunities in exploiting these huge untapped sponge resources. The report is expected to serve as an entry point for understanding Chinese sponges and for furthering R&D on their bioactive compounds for new drug development. (C) 2003 Elsevier Science B.V. All rights reserved.

ELECTROCHEMISTRY AND IN-SITU IR MICROSPECTROELECTROCHEMISTRY IN A VERSATILE THIN-LAYER REFLECTION FOURIER-TRANSFORM IR MICROSPECTROELECTROCHEMICAL CELL

A simple, convenient and versatile thin layer reflection Fourier transform IR microspectroelectrochemical (FTIRMSEC) cell has been described and characterized. Electrochemistry and in situ FTIR microspectroscopy were studied by using the hexacyanoferrate redox couple in aqueous sulphate solution, indicating that this type of cell is characteristic of both micro- or ultramicroelectrode and thin layer spectroelectrochemistry. Furthermore, the application of this FTIRMSEC cell to IR for characterization of the products of electrochemical reactions was carried out for the oxidation of (mesotetraphenylporphinato)manganese(III) perchlorate in dichloromethane + tetrabutylammonium perchlorate solution. Finally, the advantages and problems of this type of cell compared with a conventional optically transparent thin layer FTIR spectroelectrochemical cell were discussed.

Biopotentials of marine sponges from China oceans: past and future

An extensive literature survey of over 17 Journals was carried out on Chinese sponges and their natural products in the period from 1980 to 2001. This review is thus intended to provide the first thorough overview of research on marine sponges from China Ocean territories. Information is provided about the rather-limited taxonomic study of Chinese marine sponges, with an analysis on their distribution and diversity. Research findings on the natural products and their bioactivity screening from Chinese sponges are summarized. The weaknesses, gaps and problems in the past R&D program of Chinese sponges are identified, which point to the future opportunities in exploiting these huge untapped sponge resources. The report is expected to serve as an entry point for understanding Chinese sponges and for furthering R&D on their bioactive compounds for new drug development. (C) 2003 Elsevier Science B.V. All rights reserved.

藏东南巨厚松散堆积体发育条件、成灾模式及防治对策——以川藏公路然乌至鲁朗段为例

The unique geologic, geomorphic and climatic conditions of southeast Tibet have made the region to develop the multi-style and frequently occurring geologic hazards, especially the collapses and landslides and debris flows along the section of Ranwu-Lulang in Sichuan-Tibet highway. However, most of those geologic hazards have close relationship with the loose accumulations. That is, the loose accumulations are the main carrier of most geologic hazards. Thereof, the huge-thick accumulations along the highway is regarded as the objective in the thesis to study the geologic background, hazarding model and mitigation methods comprehensively, based on the multi-disciplinary theories and former materials. First of all, in the paper, based on field engineering geologic investigations, the genetic type and the characteristics of spatiotemporal distribution of the huge-thick loose accumulations along the highway, have been analysized from the factors of regional geology and geomorphy and climate, as well as the coupling acting of those factors with inoculation and eruption of the loose accumulations geologic hazards. The huge-thick loose accumulations has complex genetic types and specific regulations of spatiotemporal distribution, closely controlled by the outer environment of the region. The accumulations are composed of earth and boulder, with disorder structure and poor sorting, specific forming environments and depositing conditions. And its physical and mechanic properties are greatly distinguished from rock and common earth inland. When Sichuan-Tibet highway was firstly constructed along the north bank of Purlung Tsangpo River, the huge-thick loose accumulations was cut into many high and steep slopes. Through the survey to the cut-slopes and systematic investigation to their failures, the combination of height and angle of the accumulations slope has been obtained. At the same time, the types of genetic structure of those cut-slopes are also analysized and concluded, as well as their failure models. It is studied in the paper that there are piaster, duality, multielement and complexity types in genetic structure, and rip-dump-repose, rip-shear-slip and weathering-flake types in failure models. Moreover, it is briefly introduced present engineering performance methods and techniques dealing with the deformation and failure of the accumulations cut-slope. It is also suggested that several new techniques of slope enforcement and the method of landslide and rockfall avoiding should be applied. The research of high and steep cut-slope along the highway has broadened the acknowledgement of the combination of cut-slope height and angle. Especially, the dissertation also has made the monographic studies about the geologic background and hazarding models and prevention methods of some classic but difficult accumulations geologic hazards. They are: (1) Research of the engineering geologic background of the 102 landslide group and key problems about the project of tunnel. The 102 landslide group is a famous accumulational one composed of glacial tills and glaciofuvial deposit. The tunnel project is a feasible and optional one which can solve the present plight of “sliding after just harnessing” in the 102 section. Based on the glacial geomorphy and its depositing character, distribution of seepage line, a few drillhole materials and some surveying data, the position of contact surface between gneiss and accumulations has been recognized, and the retreating velocities of three different time scales (short, medium and long term) have been approximately calculated, and the weathering thickness of gneiss has also been estimated in the paper. On the basis of above acknowledgement, new engineering geomechnic mode is established. Numerical analysis about the stability of the No.2 landslide is done by way of FLAC program, which supplies the conclusion that the landslide there develops periodically. Thereof, 4 projects of tunnel going through the landslide have been put forwards. Safety distance of the tunnel from clinohefron has been numerically analysized. (2) Research of the geologic setting and disaster model and hazard mitigation of sliding-sand-slope. From the geologic setting of talus cone, it is indicated that the sliding-sand-slope is the process of the re-transportation and re-deposit of sand under the gravity action and from the talus cone. It is the failure of the talus cone essentially. The layering structure of the sliding-sand-slope is discovered. The models of movement and failure of the sliding-sand-slope has been put forwards. The technique, “abamurus+grass-bush fence+degradable culture pan”, is suggested to enforcement and green the sliding-sand-slope. (3) Characteristics and hazarding model and disaster mitigation of debris flow. The sources of solid material of three oversize debris flows have been analysized. It is found that a large amount of moraine existing in the glacial valley and large landslide dam-break are the two important features for oversize debris flow to be taken place. The disaster models of oversize and common debris flows have been generalized respectively. The former model better interpret the event of the Yigong super-large landslide-dam breaking. The features of common debris flow along the highway section, scouring and silting and burying and impacting, are formulated carefully. It is suggested that check dam is a better engineering structure to prevent valley from steeply scouring by debris flow. Moreover, the function of check dam in enforcing the slope is numerically calculated by FLAC program. (4) Songzong ancient ice-dammed lake and its slope stability. The lacustrine profile in Songzong landslide, more than 88 meters thick, is carefully described and measured. The Optical Simulated Luminescence (OSL) ages in the bottom and top of the silty clay layer are 22.5±3.3 kaB.P., 16.1±1.7 kaB.P., respectively. It is indicated by the ages that the lacustrine deposits formed during the Last Glacial Maximum ranging from 25ka B.P. to 15ka B.P. The special characteristics of the lacustrine sediment and the ancient lake line in Songzong basin indicated that the lacustrine sediment is related to the blocking of the Purlung Tsangpo River by the glacier in Last Glacial Maximum from Dongqu valley. The characteristics of the lacustrine profile also indicate that the Songzong ice-dammed lake might run through the Last Glacial Maximum. Two dimensional numerical modeling and analysis are done to simulate the slope stability under the conditions of nature and earthquake by FLAC program. The factor of safety of the lacusrtine slope is 1.04, but it will take place horizontal flow under earthquake activity due to the liquefaction of the 18.33 m silt layer. The realign to prevent the road from landslide is suggested.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

冀东油田柳南地区新近系油藏剩余油形成及分布规律研究

This thesis is based on the research project of Study on the Geological Characteristics and Remaining Oil Distribution Law of Neogene Reservoirs in Liunan Area, which is one of the key research projects set by PetroChina Jidong Oilfield Company in 2006. The determination of remaining oil distribution and its saturation changes are the most important research contents for the development and production modification of oilfields in high water-cut phases. Liunan oilfield, located in Tangshan of Hebei Province geographically and in Gaoliu structural belt of Nanpu sag in Bohai Bay Basin structurally, is one of the earliest fields put into production of Jidong oilfield. Focusing on the development problems encountered during the production of the field, this thesis establishes the fine geological reservoir model through the study of reservoir properties such as fine beds correlation, sedimentary facies, micro structures, micro reservoir architecture, flow units and fluid properties. Using routine method of reservoir engineering and technology of reservoir numerical modeling, remaining oil distribution in the target beds of Liunan area is predicted successfully, while the controling factors of remaining oil distribution are illustrated, and the model of remaining oil distribution for fault-block structure reservoirs is established. Using staged-subdivision reservoir correlation and FZI study, the Strata in Liunan Area is subdivided step by step; oil sand body data-list is recompiled; diagram databases are established; plane and section configuration of monolayer sandstone body, and combination pattern of sandstone bodys are summarized. The study of multi-level staged subdivision for sedimentary micro-facies shows that the Lower member of Minghuazhen formation and the whole Guantao formation in Liunan Area belong to meandering river and braided river sedimentary facies respectively, including 8 micro facies such as after point bar, channel bar, channel, natural levee, crevasse splay, abandoned channel, flood plain and flood basin. Fine 3D geological modeling is performed through the application of advanced software and integration of geological, seismic logging and reservoir engineering data. High resolution numerical simulation is performed with a reserve fitting error less than 3%, an average pressure fitting fluctuation range lower than 2Mpa and an accumulate water cut fitting error less than 5%. In this way, the distribution law of the target reservoir in the study area is basically recognized. Eight major remaining oil distribution models are established after analysis of production status and production features in different blocks and different layers. In addition, fuzzy mathematics method is used to the integreted evaluation and prediction of abundant remaining oil accumulation area in major production beds and key sedimentary time units of the shallow strata in Liunan Area and corresponding modification comments are put forward. In summary, the establishment of fine reservoir geological model, reservoir numerical simulation and distribution prediction of remaining oil make a sound foundation for further stimulation of oilfield development performance.

徐家围子断陷连片叠前时间偏移处理方法研究

In this paper, the detailed analysis of fundamental seismic data and theoretical method are given, and the tests of some new technologies are performed. For seismic data processing assembly, some key technologies are developed and applied, such as global static correction, amplitude consistency processing, wavelet consistency shaping, fine velocity model establishing and prestack time migration. These technologies can efficiently settle the problems during the course of multiple- block– jointed prestack time migration processing, and it is highly significant for holding the oil output of 40,000,000 tons for Daqing oilfield. Through the research of this dissertation, the following important contributions are shown: (1) The combination of near-surface model method and refraction static correction method is developed, and is applied to solve global static correction for the whole merging area. (2) Prestack amplitude normalization processing method based on fold is developed. The method eliminates the effects of fold on amplitude uniformity, and solves the problem of energy uniformity for tie-area prestack migration processing. (3) Wavelet consistency is investigated. For multiple survey blocks existing in the area, the optimum method of wavelet shaping is developed, which removes the waveform variance between two adjacent blocks. (4) Controlled velocity inversion (CVI) technique is used to establish migration velocity field. It can largely shorten the period of velocity modeling, and improve velocity analysis precision. (5) Float datum level technique is employed, and is able to guarantee prestack migration results of subsurface shallow layers. (6) The static partition of seismic data volume relating to migration aperture is firstly developed. And the precious imaging for huge data volume by prestack time migration is realized. (7) The numerical forward simulation and prestack migration processing is primarily combined to discuss the migration technique for a complex geology structure from practical field information. The combination of numerical simulation and prestack migration is a feasible way to solve the fine imaging of complex volcanic structure. And the combination approach can help to select appropriate migration parameters.

基于精细结构描述及数值试验的断续节理岩体力学参数智能选取

Evaluating the mechanical properties of rock masses is the base of rock engineering design and construction. It has great influence on the safety and cost of rock project. The recognition is inevitable consequence of new engineering activities in rock, including high-rise building, super bridge, complex underground installations, hydraulic project and etc. During the constructions, lots of engineering accidents happened, which bring great damage to people. According to the investigation, many failures are due to choosing improper mechanical properties. ‘Can’t give the proper properties’ becomes one of big problems for theoretic analysis and numerical simulation. Selecting the properties reasonably and effectively is very significant for the planning, design and construction of rock engineering works. A multiple method based on site investigation, theoretic analysis, model test, numerical test and back analysis by artificial neural network is conducted to determine and optimize the mechanical properties for engineering design. The following outcomes are obtained: (1) Mapping of the rock mass structure Detailed geological investigation is the soul of the fine structure description. Based on statistical window，geological sketch and digital photography，a new method for rock mass fine structure in-situ mapping is developed. It has already been taken into practice and received good comments in Baihetan Hydropower Station. (2) Theoretic analysis of rock mass containing intermittent joints The shear strength mechanisms of joint and rock bridge are analyzed respectively. And the multiple modes of failure on different stress condition are summarized and supplied. Then, through introducing deformation compatibility equation in normal direction, the direct shear strength formulation and compression shear strength formulation for coplanar intermittent joints, as well as compression shear strength formulation for ladderlike intermittent joints are deducted respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. (3) Model test of rock mass containing intermittent joints Model tests are adopted to study the mechanical mechanism of joints to rock masses. The failure modes of rock mass containing intermittent joints are summarized from the model test. Six typical failure modes are found in the test, and brittle failures are the main failure mode. The evolvement processes of shear stress, shear displacement, normal stress and normal displacement are monitored by using rigid servo test machine. And the deformation and failure character during the loading process is analyzed. According to the model test, the failure modes quite depend on the joint distribution, connectivity and stress states. According to the contrastive analysis of complete stress strain curve, different failure developing stages are found in the intact rock, across jointed rock mass and intermittent jointed rock mass. There are four typical stages in the stress strain curve of intact rock, namely shear contraction stage, linear elastic stage, failure stage and residual strength stage. There are three typical stages in the across jointed rock mass, namely linear elastic stage, transition zone and sliding failure stage. Correspondingly, five typical stages are found in the intermittent jointed rock mass, namely linear elastic stage, sliding of joint, steady growth of post-crack, joint coalescence failure, and residual strength. According to strength analysis, the failure envelopes of intact rock and across jointed rock mass are the upper bound and lower bound separately. The strength of intermittent jointed rock mass can be evaluated by reducing the bandwidth of the failure envelope with geo-mechanics analysis. (4) Numerical test of rock mass Two sets of methods, i.e. the distinct element method (DEC) based on in-situ geology mapping and the realistic failure process analysis (RFPA) based on high-definition digital imaging, are developed and introduced. The operation process and analysis results are demonstrated detailedly from the research on parameters of rock mass based on numerical test in the Jinping First Stage Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Then the applicable fields are figured out respectively. (5) Intelligent evaluation based on artificial neural network (ANN) The characters of both ANN and parameter evaluation of rock mass are discussed and summarized. According to the investigations, ANN has a bright application future in the field of parameter evaluation of rock mass. Intelligent evaluation of mechanical parameters in the Jinping First Stage Hydropower Station is taken as an example to demonstrate the analysis process. The problems in five aspects, i. e. sample selection, network design, initial value selection, learning rate and expected error, are discussed detailedly.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

非均匀性对岩石破裂过程影响的研究

Rock heterogeneity plays an important role in rock fracturing processes. However, because fracturing is a dynamic process and it is very difficult to quantify materials' heterogeneity, most of the theories dealing with local failure were based on the homogeneity assumption, very few involving stress distribution heterogeneity and successive local failure due to rock heterogeneity. Therefore, based on various references, the author studied the laws and mechanism of influences of heterogeneity on rock fracturing processes, under the frame of the project "Study on Associate Mechanism between Rock Mass Fracture and Strength Failure", funded by Nation Natural Science Fund. the research consists of such aspects as size effect correction to rock fracture parameters, SEM (Scanning Electron Microscope) real-time observation on rock samples under different loads, micro-hardness testing, and numerical simulating based on microstructure. There are some important research results as followed: 1. Unifying formula for nonlinear and non-singularity correction, simplifying the complex process of correcting size effect on rock fracture toughness. 2. Using the methods of micro-hardness testing mineral grain and random jointing micrograph digitizing mineral slice, preliminarily solving the problems of numerical simulating and quantitatively describing the heterogeneous strength and its distribution rules, which has certain innovation and better practicability. 3. Based on SEM real-time observation, studying the micro-process of fracturing in marble, sandstone, granite, and mushroom stone samples with premanufactured cracks under tension, pure-shear and compression-shear conditions. Strength Failure was observed: there was some kind failure occurred before Fracture Failure in marble and sandstone samples with double cracks under pure-shearing. It is believed that the reason of strength failure developing is that stress concentrations is some locations are larger than that near the end of pre-manufactured cracks. 4. Based on the idea that rock macro-constitute is composed of complex microstructure, the promising method used to handle heterogeneity considers not only the heterogeneity of the rock medium, but also the heterogeneity of the rock structure. 5. Putting forward two types of rock strength failure: medium strength failure induced by heterogeneity of rock medium and structure strength failure induced by heterogeneity rock structure. 6. By evaluating potential fracture cell with proper failure priority, the numerical simulating method solved the problem of simulating the coextensive strength failure and fracture failure with convention strength failure rules. The result of numerical analysis shows that the influence of heterogeneity on rock fracturing processes is evident. The sinuosity of the rock fracture-propagation path, and the irregular fluctuation of loading displacement curve, is mainly controlled by the heterogeneity of rock medium.

改善三维波动方程叠前深度偏移处理效果的方法研究

3D wave equation prestack depth migration is the effective tool for obtaining the exact imaging result of complex geology structures. It's a part of the 3D seismic data processing. 3D seismic data processing belongs to high dimension signal processing, and there are some difficult problems to do with. They are: How to process high dimension operators? How to improve the focusing? and how to construct the deconvolution operator? The realization of 3D wave equation prestack depth migration, not only realized the leap from poststack to prestack, but also provided the important means to solve the difficult problems in high dimension signal processing. In this thesis, I do a series research especially for the solve of the difficult problems around the 3D wave equation prestack depth migration and using it as a mean. So this thesis service for the realization of 3D wave equation prestack depth migration for one side and improve the migration effect for another side. This thesis expatiates in five departs. Summarizes the main contents as the follows: In the first part, I have completed the projection from 3D data point area to low dimension are using de big matrix transfer and trace rearrangement, and realized the liner processing of high dimension signal. Firstly, I present the mathematics expression of 3D seismic data and the mean according to physics, present the basic ideal of big matrix transfer and describe the realization of five transfer models for example. Secondly, I present the basic ideal and rules for the rearrange and parallel calculate of 3D traces, and give a example. In the conventional DMO focusing method, I recall the history of DM0 process firstly, give the fundamental of DMO process and derive the equation of DMO process and it's impulse response. I also prove the equivalence between DMO and prestack time migration, from the kinematic character of DMO. And derive the relationship between DMO base on wave equation and prestack time migration. Finally, I give the example of DMO process flow and synthetic data of theoretical models. In the wave equation prestak depth migration, I firstly recall the history of migration from time to depth, from poststack to prestack and from 2D to 3D. And conclude the main migration methods, point out their merit and shortcoming. Finally, I obtain the common image point sets using the decomposed migration program code.In the residual moveout, I firstly describe the Viterbi algorithm based on Markov process and compound decision theory and how to solve the shortest path problem using Viterbi algorithm. And based on this ideal, I realized the residual moveout of post 3D wave equation prestack depth migration. Finally, I give the example of residual moveout of real 3D seismic data. In the migration Green function, I firstly give the concept of migration Green function and the 2D Green function migration equation for the approximate of far field. Secondly, I prove the equivalence of wave equation depth extrapolation algorithms. And then I derive the equation of Green function migration. Finally, I present the response and migration result of Green function for point resource, analyze the effect of migration aperture to prestack migration result. This research is benefit for people to realize clearly the effect of migration aperture to migration result, and study on the Green function deconvolution to improve the focusing effect of migration.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.

井间地震技术研究和应用

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.

Supporting Reuse and Evolution in Software Design

Program design is an area of programming that can benefit significantly from machine-mediated assistance. A proposed tool, called the Design Apprentice (DA), can assist a programmer in the detailed design of programs. The DA supports software reuse through a library of commonly-used algorithmic fragments, or cliches, that codifies standard programming. The cliche library enables the programmer to describe the design of a program concisely. The DA can detect some kinds of inconsistencies and incompleteness in program descriptions. It automates detailed design by automatically selecting appropriate algorithms and data structures. It supports the evolution of program designs by keeping explicit dependencies between the design decisions made. These capabilities of the DA are underlaid bya model of programming, called programming by successive elaboration, which mimics the way programmers interact. Programming by successive elaboration is characterized by the use of breadth-first exposition of layered program descriptions and the successive modifications of descriptions. A scenario is presented to illustrate the concept of the DA. Technques for automating the detailed design process are described. A framework is given in which designs are incrementally augmented and modified by a succession of design steps. A library of cliches and a suite of design steps needed to support the scenario are presented.