辽东湾隐蔽油气藏地球物理识别技术研究

In China and world, more than half the recent basin discovered reserves involve lithologic hydrocarbon reservoir reserves. The major target for further hydrocarbon basin exploration is the subtle reservoir. The Liaodong Bay prospect is much important in Bohai Sea, which includes Liaoxi low uplift, Liaodong uplift, Liaoxi sag and Liaozhong sag. After dozens years’ exploration in Liaodong Bay, few unexplored big-and-middle-sized favorable structural traps are remained and most of the stock structure targets are bad for fragmentary. Thus seeking for new prospect area and making a breakthrough, have become the unique way to relieve the severe exploration condition in Liaodong Bay. Technique Route Based on the petrophysical property of target area, the seismic forward inference of typical subtle trap model is expanded with analysis of logging, seismic and geologic data. According to petrophysical characteristics and forward inference and research on seismic response of actual seismic data in target area, the optimization of geophysical technique is used in subtle trap identification and the geophysical identification technique system of subtle reservoir is formed. The Key Research ① Petrophysical Model The petrophysical parameter is the basic parameter for seismic wave simulation. The seismic response difference of rocks bearing different fluids is required. With the crossplot of log data, the influence of petrophysical parameters on rock elastic properties of target area is analyzed, such as porosity, shale index, fluid property and saturation. Based on the current research on Biot-Gassmann and Kuster-Toksoz model, the petrophysical parameter calculator program which can be used for fluid substitution is established. ② S-wave evaluation based on conventional log data The shear velocity is needed during forward inference of AVO or other elastic wave field. But most of the recent conventional log data is lack of shear wave. Thus according to the research on petrophysical model, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. ③ AVO forward modeling based on well data For 6 wells in JZ31-6 block and 9 wells in LD22-1 block, the AVO forward modeling recording is made by log curve. The classification of AVO characteristics in objective interval is made by the lithologic information. ④ The 2D parameter model building and forward modeling of subtle hydrocarbon trap in target area. According to the formation interpretation of ESS03D seismic area, the 2D parameter model building and seismic wave field forward modeling are carried on the given and predicted subtle hydrocarbon trap with log curve. ⑤ The lithology and fluid identification of subtle trap in target area After study the seismic response characteristics of lithology and fluid in given target area, the optimization of geophysical technique is used for lithology identification and fluid forecast. ⑥The geophysical identification technique system of subtle reservoir The Innovative Points of this Paper ① Based on laboratory measurement and petrophysical model theory, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. Then the fluid substitution method based on B-G and K-T theory is provided. ② The method and workflow for simulating seismic wave field property of subtle hydrocarbon trap are established based on the petrophysical model building and forward modeling of wave equation. ③ The description of subtle trap structural feature is launched. According to the different reflection of frequency wave field structural attribute, the fluid property of subtle trap can be identified by wave field attenuation attribute and absorption analysis. ④ It’s the first time to identify subtle trap by geophysical technique and provide exploration drilling well location. ⑤ The technique system of subtle reservoir geophysical identification is formed to provide available workflow and research ideas for other region of interest.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

Z-scan theory based on a diffraction model

Based on Fresnel-Mrchhoff diffraction theory, a diffraction model of nonlinear optical media interacting with a Gaussian beam has been set up that can interpret the Z-scan phenomenon in a new way. This theory not only is consistent with the conventional Z-scan theory for a small nonlinear phase shift but also can be used for larger nonlinear phase shifts. Numerical computations indicate that the shape of the Z-scan curve is greatly affected by the value of the nonlinear phase shift. The symmetric dispersionlike Z-scan curve is valid only for small nonlinear p base shifts (\Deltaphi(0)\ < pi), but, with increasingly larger nonlinear phase shifts, the valley of the transmittance is severely suppressed and the peak is greatly enhanced. The power output through the aperture will oscillate with increasing nonlinear phase shift caused by the input laser power. The aperture transmittance will attenuate and saturate with increasing Kerr constant. (C) 2003 Optical Society of America.

Biomimetic (topological) pattern recognition - A new model of pattern recognition theory and its application

A new theoretical model of Pattern Recognition principles was proposed, which is based on "matter cognition" instead of "matter classification" in traditional statistical Pattern Recognition. This new model is closer to the function of human being, rather than traditional statistical Pattern Recognition using "optimal separating" as its main principle. So the new model of Pattern Recognition is called the Biomimetic Pattern Recognition (BPR)(1). Its mathematical basis is placed on topological analysis of the sample set in the high dimensional feature space. Therefore, it is also called the Topological Pattern Recognition (TPR). The fundamental idea of this model is based on the fact of the continuity in the feature space of any one of the certain kinds of samples. We experimented with the Biomimetic Pattern Recognition (BPR) by using artificial neural networks, which act through covering the high dimensional geometrical distribution of the sample set in the feature space. Onmidirectionally cognitive tests were done on various kinds of animal and vehicle models of rather similar shapes. For the total 8800 tests, the correct recognition rate is 99.87%. The rejection rate is 0.13% and on the condition of zero error rates, the correct rate of BPR was much better than that of RBF-SVM.

Strain gradient theory based on a new framework of non-local model

A new framework of non-local model for the strain energy density is proposed in this paper. The global strain energy density of the representative volume element is treated as a non-local variable and can be obtained through a special integral of the local strain energy density. The local strain energy density is assumed to be dependent on both the strain and the rotation-gradient. As a result of the non-local model, a new strain gradient theory is derived directly, in which the first and second strain gradients, as well as the triadic and tetradic stress, are introduced in the context of work conjugate. For power law hardening materials, size effects in thin metallic wire torsion and ultra-thin cantilever beam bend are investigated. It is found that the result predicted by the theoretical model is well consistent with the experimental data for the thin wire torsion. On the other hand, the calculation result for the micro-cantilever beam bend clearly shows the size effect.

A ionization theory based on united and separated atom model

The direct Coulomb ionization process can be generally well described by the ECPSSR theory, which bases on the perturbed-stationary- state(PSS) and accounts for the energy-loss, Coulomb-deflection, and relativistic effects. But the ECPSSR calculation has significant deviations for heavy projectile at low impinging energies. In this paper we propose a new modified ECPSSR theory, i.e. MECUSAR, in which PSS is replaced by an united and separated atom model, and molecule-orbit effect is considered. The MECUSAR calculations give better agreement with the experimental data at lower impinging energies, and agree with the ECPSSR calculations at high energies. By using OBKN (Oppenheimer-Brinkman-Kramers formulas of Nikolaev) theory to describe the contribution of the electron capture, we further modified the proposed MECUSAR theory, and calculated the target ionization cross sections for different charge states of the projectile.

BRST invariant theory of a generalized 1+1 Dimensional nonlinear sigma model with topological term

We give a generalized Lagrangian density of 1 + 1 Dimensional O( 3) nonlinear sigma model with subsidiary constraints, different Lagrange multiplier fields and topological term, find a lost intrinsic constraint condition, convert the subsidiary constraints into inner constraints in the nonlinear sigma model, give the example of not introducing the lost constraint. N = 0, by comparing the example with the case of introducing the lost constraint, we obtain that when not introducing the lost constraint, one has to obtain a lot of various non-intrinsic constraints. We further deduce the gauge generator, give general BRST transformation of the model under the general conditions. It is discovered that there exists a gauge parameter beta originating from the freedom degree of BRST transformation in a general O( 3) nonlinear sigma model, and we gain the general commutation relations of ghost field.

多分量对偶超导理论和口袋模型；Multi-component Dual Superconductor Theory and Bag Model

基于联络新参数化方案研究了多分量对偶超导模型。给出了多分量Ginzburg-Landau模型中的自对偶解,并研究了磁通量子数趋于无穷大时的墙涡旋解,以及与口袋模型之间的联系。

Deformation Monitoring and Forecasting Model Study Based on Grey System Theory

The grey system theory studies the uncertainty of small sample size problems. This paper using grey system theory in the deformation monitoring field, based on analysis of present grey forecast models, developed the spatial multi-point model. By using residual modification, the spatial multi-point residual model eras developed in further study. Then, combined with the sedimentation data of Xiaolangdi Multipurpose Dam, the results are compared and analyzed, the conclusion has been made and the advantages of the residual spatial multi-point model has been proved.

Application of the Sanchez-Lacombe lattice fluid theory to the system pvme/ps and model calculations

Cloud point curves reported in the literature for five representatives of the system poly(vinyl methyl ether)/polystyrene were evaluated theoretically by means of the Sanchez-Lacombe lattice fluid theory. The measured phase separation behavior can be described within experimental error using only one adjustable parameter (quantifying the interaction between the unlike mers). The Flory-Huggins interaction parameters calculated from this theoretical description depend in good approximation linearly on composition (volume fractions) and on the inverse temperature. An evaluation of these data yields a maximum heat effect which is almost one order of magnitude less (ca. -0.25 J/cm(3)) than obtained via Hess's cycle (dissolution of the components and of the blend) from calorimetric measurements. Model calculations on the basis of the present theory demonstrate that the critical points shift to a different extent upon a certain relative change in the molar mass of the blend components. The sensitivity of the calculated phase diagrams against changes in the scaling parameter decreases in the following order: interaction energies between unlike mers, differences in the scaling temperatures, pressures and densities.

A meso elastoplastic constitutive model for polycrystalline metals based on equivalent slip systems with latent hardening

A meso material model for polycrystalline metals is proposed, in which the tiny slip systems distributing randomly between crystal slices in micro-grains or on grain boundaries are replaced by macro equivalent slip systems determined by the work-conjugate principle. The elastoplastic constitutive equation of this model is formulated for the active hardening, latent hardening and Bauschinger effect to predict macro elastoplastic stress-strain responses of polycrystalline metals under complex loading conditions. The influence of the material property parameters on size and shape of the subsequent yield surfaces is numerically investigated to demonstrate the fundamental features of the proposed material model. The derived constitutive equation is proved accurate and efficient in numerical analysis. Compared with the self-consistent theories with crystal grains as their basic components, the present theory is much simpler in mathematical treatment.

Peierls-Nabarro model of interfacial misfit dislocation: An analytic solution

We propose a method to treat the interfacial misfit dislocation array following the original Peierls-Nabarro's ideas. A simple and exact analytic solution is derived in the extended Peierls-Nabarro's model, and this solution reflects the core structure and the energy of misfit dislocation, which depend on misfit and bond strength. We also find that only with beta < 0.2 the structure of interface can be represented by an array of singular Volterra dislocations, which conforms to those of atomic simulation. Interfacial energy and adhesive work can be estimated by inputting ab initio calculation data into the model, and this shows the method can provide a correlation between the ab initio calculations and elastic continuum theory.

Source Model for Blasting Vibration

By analyzing and comparing the experimental data, the point source moment theory and the cavity theory, it is concluded that the vibrating signals away from the blasting explosive come mainly from the natural vibrations of the geological structures near the broken blasting area. The source impulses are not spread mainly by the inelastic properties (such as through media damping, as believed to be the case by many researchers) of the medium in the propagation pass, but by this structure. Then an equivalent source model for the blasting vibrations of a fragmenting blasting is proposed, which shows the important role of the impulse of the source's time function under certain conditions. For the purpose of numerical simulation, the model is realized in FEM, The finite element results are in good agreement with the experimental data.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.