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

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

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

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

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

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

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

A novel method to synthesize amorphous silica-alumina materials with mesoporous distribution without using templates and pore-regulating agents

In this study, a novel sol-gel method is used to synthesize amorphous silica-alumina materials with a narrow mesoporous distribution and various Si/Al molar ratios without using any templates and pore-regulating agents. During the preparation procedure, only inexpensive inorganic salts were used as raw materials, instead of expensive and harmful alkoxides. The precursor sol was dried at room temperature in a vacuum box kept at 60 mmHg until it began to form the gel. The results of a nitrogen sorption experiment indicate that the synthesized materials with different Si/Al molar ratios have similar mesoporous distributions (within 2-12 nm). Moreover, it was found that the material's pore size distribution remains at a similar value during the heat treatment from room temperature to 550 degreesC. On the basis of the nitrogen sorption, TEM, and AFM characterization results, a formation mechanism of mesopores which accounts for the experimental data is also suggested. This suggested mechanism involves rearrangement of the primary particles during the drying process to form the precursors of the similarly sized mesopores. The synthesized materials were characterized by XRD, thermal analysis (TG/DTA), Al-27 and Si-29 MAS NMR spectroscopy, SEM, TEM, and AFM. The results of Al-27 and 29Si MAS NMR indicate that the distribution of silicon and aluminum in the synthesized materials is more uniform and homogeneous than that in the mixed oxides prepared via the traditional sol-gel method even at high alumina contents. The type and density of the acid sites were studied using pyridine adsorption-desorption FTIR spectroscopy. It was shown that the acidity of the synthesized materials is higher than that of the silica-alumina materials prepared by conventional methods.

Catalytic oxidation of cyclohexane to cyclohexanol and cyclohexanone over Co3O4 nanocrystals with molecular oxygen

Co3O4 nanocrystals with average particle sizes of 30 and 50 run were synthesized using cobalt nitrate as precursor, and were characterized by X-ray diffraction (XRD), nitrogen adsorption, transmission electron microscopy (TEM), and Fourier transform infrared (FT-IR) spectroscopy. Catalytic oxidation of cyclohexane with molecular oxygen was studied over Co3O4 nanocrystals. These catalysts showed obviously higher activities as compared to Co3O4 prepared by the conventional methods, Co3O4/Al2O3, or homogeneous cobalt catalyst under comparable reaction conditions. The 89.1% selectivity to cyclohexanol and cyclohexanone at 7.6% conversion of cyclohexane was realized over 50 nm sized Co3O4 nanocrystals at 393 K for 6 h. (c) 2005 Elsevier B.V. All rights reserved.

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

Assessment of bursting detection methods with LES database of turbulent channel flows

Validated by comparison with DNS, numerical database of turbulent channel flows is yielded by Large Eddy Simulation (LES). Three conventional techniques: uv quadrant 2, VITA and mu-level techniques for detecting turbulent bursts are applied to the identification of turbulent bursts. With a grouping parameter introduced by Bogard & Tiedemann (1986) or Luchik & Tiederman (1987), multiple ejections detected by these techniques which originate from a single burst can be grouped into a single-burst event. The results are compared with experimental results, showing that all techniques yield reasonable average burst period. However, uv quadrant 2 and mu-level are found to be superior to VITA in having large threshold-independent range.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Improved optical heterodyne methods for measuring. frequency responses of photodetectors

An improved optical self-heterodyne method utilizing a distributed Bragg reflector (DBR) tunable laser and an optical fiber ring interferometer is presented in this paper. The interference efficiency can be increased by 7 dB compared with the scheme using the conventional Mach-Zehnder interferometer. The unsteady process that the beating frequency experiences in each tuning period is investigated. According to the measurement results, the wavelength and optical power of the tunable laser will be steady when the square-wave frequency is lower than 300 kHz. It has been shown that when a square-wave voltage is applied to the phase section of the tunable laser, the laser linewidths vary in a wide range, and are much larger than that under dc voltage tuning. The errors caused by the variations in the linewidth of the beat signal and optical power can be eliminated using the proposed calibration procedures, and the measurement accuracy can, therefore, be significantly improved. Experiments show that the frequency responses obtained using our method agree well with the data provided by the manufacturer, and the improved optical self-heterodyne method is as accurate as the intensity noise technique.

New algorithms of the TSM and TOM methods for calibrating microwave test fixtures

Based on the conventional through-short-match (TSM) method, an improved TSM method has been proposed in this Letter. This method gives an analytical solution and has almost all the advantages of conventional TSM methods. For example, it has no phase uncertainty and no bandwidth limitation. The experimental results show that the accuracy can be significantly improved with this method. The proposed theory can be applied to the through-open-match (TOM) method. (C) 2002 Wiley Periodicals. Inc.

INVESTIGATION OF C-60 SINGLE-CRYSTAL BY X-RAY-METHODS

C-60 Single crystals grown by a single-temperature-gradient technique were characterized by synchrotron radiation white beam x-ray topography and x-ray double crystal diffraction with Cu K-alpha 1 radiation on conventional x-ray source. The results show that the crystal is rather well crystallized, The x-ray topographies give an evidence of dendritic growth mechanism of C-60 Single crystal, and x-ray double crystal diffraction rocking curve shows that there are mosaic structural defects in the sample. A phase transition st 249+/-1.5% K from a simple cubic to a face centered cubic structure is confirmed by in situ observation of synchrotron radiation white beam x-ray topography with the temperature varing from 230 to 295 K.

Effect of drying methods on the structure and catalytic combustion activity of Mn-substituted hexaaluminate catalysts

Conventional oven drying (COD) and supercritical drying (SCD) methods were applied to the preparation of Mn-substituted hexaaluminate (BaMnA(11)O(19-alpha)) catalysts. The effect of drying methods on phase composition, specific surface area, pore structure and combustion activity of the samples was investigated. The samples obtained by SCD have higher surface area, narrower pore size distribution, and higher combustion activity than those obtained by COD.

COMPARATIVE-STUDY OF CA4Y6(SIO4)(6)O-A PHOSPHORS PREPARED BY SOL-GEL AND DRY METHODS (A=PB2+,EU3+, TB3+,DY3+)

Ca4Y6(SiO4)(6)O:A (A = Pb2+, Eu3+, Tb3+, Dy3+) phosphors have been prepared by two methods: the sol-gel method and the conventional dry method. The crystallization processes and the luminescence characteristics of the phosphors were studied, The sol-gel method features low-temperature formation of the phosphor, leading to successful preparation of Pb2+-activated phosphors which could not be prepared by the dry method at high temperature. The (4f)(8-)(4f)(7)(5d)(1) absorption band of Tb3+ and the charge-transfer (CT) band of Eu3+ have higher energies and narrower half-widths in the sol-gel-derived phosphors than in the phosphors prepared by the dry method, respectively. The Tb3+ and Dy3+ ions show stronger emission in the former than in the latter. Both the yellow-to-blue intensity ratio (Y:B) of Dy3+ and the red-to-orange intensity ratio (R:O) of Eu3+ in the sol-gel-derived phosphors are smaller than those for the phosphors derived by the dry method.