三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

辽河盆地东部凹陷热历史及其与油气关系

This dissertation focuses on the basin geothermal history, tectonothermal evolution and the relationship between geothermal field evolution and hydrocarbon generation. Based on the research of present-day geothermal field, geothermal history of Eastern Subdepression of Liaohe Basin was reconstructed with available data from drillings, loggings, seismic cross-sections, BHTs and thermal indicators. 12 heat flow density values were calculated. Ranging from 45.7 mW/m~2 to 70.0 mW/m~2, the mean value of these determinations exhibits 58.0(±5.83mW/m~22). The heat flow density in the uplift and ramp is greater than that in the sag. The main factors affecting the heat flow density are the morphostructure of basement and thickness of sedimentary cover. The Subdepression experienced a much higher heat flow period from 43 Ma to about 25 Ma. The heat flow increased gradually from Sha3 stage to Dongying stage, and reached the peak at the late of Dongying stage, then cooled down. Structural subsidence analysis shows that the subsidence of Eastern Subdepression can be divided into two phases: earlier (25-43Ma) initial (rift) and late (since 25Ma) thermal subsidence. The lower present-day heat flow and the higher palaeo-heat flow corresponding to structural subsidence stage as well as the typical rift subsidence style in Eastern Subdepression provide with some insights to the tectonic-thermal evolution of the basin. The source rocks of Sha3 member began to generate oil in the Shal stage, and entered oil-window at the late of Dongying stage. The source rocks of Shal began to generate oil at the late of Dongying stage, and being at the stage of lower maturation-maturation now. Most of Dongying source rocks are not mature now. The late of Dongying stage is the critical time for the oil system.

环渤海地区地壳和上地幔地震层析成像研究

In Circum-Bohai region (112°~124°E, 34°~42°N ), there exists rich gas-petroleum while inner-plate seismic activity is robust. Although the tectonic structure of this region is very complicated, plenty of geological, geophysical and geochemical researches have been carried out.In this paper, guided by the ideas of "One, Two, Three and Many" and "The depth controls the shallow, the regional constrains the local", I fully take advantage of previous results so as to establish a general image of this region. After collecting the arrival-time of P-wave phases of local events and tele-seismic events recorded by the stations within this region from 1966 to 2004, I process all these data and build an initial model. Then, a tomography image of crust and upper-mantle of this region is obtained. With reference to previous results, we compare the image of various depths and five cross-profiles traverse this region along different direction. And finally, a discussion and conclusion is made.The principle contents is listed as below: 1) in the first chapter, the purpose and meaning of this thesis, the advance in seismic tomography, and the research contents and blue-print is stated; 2) in the second chapter, I introduce the regional geological setting of Circum-Bohai region, describe the tectonic and evolutionary characteristics of principle tectonic units, including Bohai Bay Basin, Yanshan Fold Zone, Taihangshan Uplifted Zone, Jiao-Niao Uplifted Zone and Luxi Uplifted Zone, and primary deep faults; 3) In the third chapter, the previous geophysical researches, i.e., gravity and geomagnetic characters, geothermal flow, seismic activity, physical character of rocks, deep seismic sounding, and previous seismic tomography, are discussed; 4) in the fourth chapter, the fundamental theory and approach of seismic tomography is introduced; 5) in the fifth chapter, the technology and approaches used in this thesis, including collecting and pre-processing of data, the establishment of initial velocity model and relocation of all events; 6) in the sixth chapter, I discuss and analyze the tomography image of various depth and five cross-sections; 7)in the seventh chapter, I make a conclusion of the results, state the existing problems and possible solutions.

Broadband optical cavity absorption spectroscopy in the near-ultraviolet: applications in atmospheric chemistry

A novel spectroscopic method, incoherent broadband cavity enhanced absorption spectroscopy (IBBCEAS), has been modified and extended to measure absorption spectra in the near-ultraviolet with high sensitivity. The near-ultraviolet region extends from 300 to 400 nm and is particularly important in tropospheric photochemistry; absorption of near-UV light can also be exploited for sensitive trace gas measurements of several key atmospheric constituents. In this work, several IBBCEAS instruments were developed to record reference spectra and to measure trace gas concentrations in the laboratory and field. An IBBCEAS instrument was coupled to a flow cell for measuring very weak absorption spectra between 335 and 375 nm. The instrument was validated against the literature absorption spectrum of SO2. Using the instrument, we report new absorption cross-sections of O3, acetone, 2-butanone, and 2-pentanone in this spectral region, where literature data diverge considerably owing to the extremely weak absorption. The instrument was also applied to quantifying low concentrations of the short-lived radical, BrO, in the presence of strong absorption by Br2 and O3. A different IBBCEAS system was adapted to a 4 m3 atmosphere simulation chamber to record the absorption cross-sections of several low vapour pressure compounds, which are otherwise difficult to measure. Absorption cross-sections of benzaldehyde and the more volatile alkyl nitrites agree well with previous spectra; on this basis, the cross-sections of several nitrophenols are reported for the first time. In addition, the instrument was also used to study the optical properties of secondary organic aerosol formed following the photooxidation of isoprene. An extractive IBBCEAS instrument was developed for detecting HONO and NO2 and had a sensitivity of about 10-9 cm-1. This instrument participated in a major international intercomparison of HONO and NO2 measurements held in the EUPHORE simulation chamber in Valencia, Spain, and results from that campaign are also reported here.

Simulation of multigate SOI transistors with silicon, germanium and III-V channels

In this work by employing numerical three-dimensional simulations we study the electrical performance and short channel behavior of several multi-gate transistors based on advanced SOI technology. These include FinFETs, triple-gate and gate-all-around nanowire FETs with different channel material, namely Si, Ge, and III-V compound semiconductors, all most promising candidates for future nanoscale CMOS technologies. Also, a new type of transistor called “junctionless nanowire transistor” is presented and extensive simulations are carried out to study its electrical characteristics and compare with the conventional inversion- and accumulation-mode transistors. We study the influence of device properties such as different channel material and orientation, dimensions, and doping concentration as well as quantum effects on the performance of multi-gate SOI transistors. For the modeled n-channel nanowire devices we found that at very small cross sections the nanowires with silicon channel are more immune to short channel effects. Interestingly, the mobility of the channel material is not as significant in determining the device performance in ultrashort channels as other material properties such as the dielectric constant and the effective mass. Better electrostatic control is achieved in materials with smaller dielectric constant and smaller source-to-drain tunneling currents are observed in channels with higher transport effective mass. This explains our results on Si-based devices. In addition to using the commercial TCAD software (Silvaco and Synopsys TCAD), we have developed a three-dimensional Schrödinger-Poisson solver based on the non-equilibrium Green’s functions formalism and in the framework of effective mass approximation. This allows studying the influence of quantum effects on electrical performance of ultra-scaled devices. We have implemented different mode-space methodologies in our 3D quantum-mechanical simulator and moreover introduced a new method to deal with discontinuities in the device structures which is much faster than the coupled-mode-space approach.

Towards more accurate and reliable predictions for nuclear applications

The need for nuclear data far from the valley of stability, for applications such as nuclear as- trophysics or future nuclear facilities, challenges the robustness as well as the predictive power of present nuclear models. Most of the nuclear data evaluation and prediction are still performed on the basis of phenomenological nuclear models. For the last decades, important progress has been achieved in funda- mental nuclear physics, making it now feasible to use more reliable, but also more complex microscopic or semi-microscopic models in the evaluation and prediction of nuclear data for practical applications. In the present contribution, the reliability and accuracy of recent nuclear theories are discussed for most of the relevant quantities needed to estimate reaction cross sections and beta-decay rates, namely nuclear masses, nuclear level densities, gamma-ray strength, fission properties and beta-strength functions. It is shown that nowadays, mean-field models can be tuned at the same level of accuracy as the phenomenological mod- els, renormalized on experimental data if needed, and therefore can replace the phenomenogical inputs in the prediction of nuclear data. While fundamental nuclear physicists keep on improving state-of-the-art models, e.g. within the shell model or ab initio models, nuclear applications could make use of their most recent results as quantitative constraints or guides to improve the predictions in energy or mass domain that will remain inaccessible experimentally.

The HITRAN 2008 molecular spectroscopic database

This paper describes the status of the 2008 edition of the HITRAN molecular spectroscopic database. The new edition is the first official public release since the 2004 edition, although a number of crucial updates had been made available online since 2004. The HITRAN compilation consists of several components that serve as input for radiative-transfer calculation codes: individual line parameters for the microwave through visible spectra of molecules in the gas phase; absorption cross-sections for molecules having dense spectral features, i.e. spectra in which the individual lines are not resolved; individual line parameters and absorption cross-sections for bands in the ultraviolet; refractive indices of aerosols, tables and files of general properties associated with the database; and database management software. The line-by-line portion of the database contains spectroscopic parameters for 42 molecules including many of their isotopologues. © 2009 Elsevier Ltd.

The HITRAN molecular spectroscopic database: Edition of 2000 including updates through 2001

This paper describes the status circa 2001, of the HITRAN compilation that comprises the public edition available through 2001. The HITRAN compilation consists of several components useful for radiative transfer calculation codes: high-resolution spectroscopic parameters of molecules in the gas phase, absorption cross-sections for molecules with very dense spectral features, aerosol refractive indices, ultraviolet line-by-line parameters and absorption cross-sections, and associated database management software. The line-by-line portion of the database contains spectroscopic parameters for 38 molecules and their isotopologues and isotopomers suitable for calculating atmospheric transmission and radiance properties. Many more molecular species are presented in the infrared cross-section data than in the previous edition, especially the chlorofluorocarbons and their replacement gases. There is now sufficient representation so that quasi-quantitative simulations can be obtained with the standard radiance codes. In addition to the description and justification of new or modified data that have been incorporated since the last edition of HITRAN (1996), future modifications are indicated for cases considered to have a significant impact on remote-sensing experiments. © 2003 Elsevier Ltd. All rights reserved.

Phase-accuracy comparisons and improved far-field estimates for 3-D edge elements on tetrahedral meshes

Edge-element methods have proved very effective for 3-D electromagnetic computations and are widely used on unstructured meshes. However, the accuracy of standard edge elements can be criticised because of their low order. This paper analyses discrete dispersion relations together with numerical propagation accuracy to determine the effect of tetrahedral shape on the phase accuracy of standard 3-D edgeelement approximations in comparison to other methods. Scattering computations for the sphere obtained with edge elements are compared with results obtained with vertex elements, and a new formulation of the far-field integral approximations for use with edge elements is shown to give improved cross sections over conventional formulations.

Single and double detachment from H-

Absolute cross sections for single and double detachment from H– following electron impact have been measured over a range of collision energies from the thresholds to 170 eV. The measurements were made using a magnetic storage ring. The ions in the ring were merged with a monoenergetic electron beam and neutral and positively charged fragments were detected. We cover larger energy ranges than in many of the previous experiments, and this is the first time both single and double detachment have been measured simultaneously. This allows us to present accurate ratios between the single and double detachment cross sections. On the basis of these ratio measurements we discuss possible mechanisms leading to double detachment.

Effect of excitation wavelength on the Raman spectroscopy of the porcine photoreceptor layer from the area centralis.

Raman microscopy, based upon the inelastic scattering (Raman) of light by molecular species, has been applied as a specific structural probe in a wide range of biomedical samples. The purpose of the present investigation was to assess the potential of the technique for spectral characterization of the porcine outer retina derived from the area centralis, which contains the highest proportion of cone:rod cell ratio in the pig retina. METHODS: Retinal cross-sections, immersion-fixed in 4% (w/v) PFA and cryoprotected, were placed on salinized slides and air-dried prior to direct Raman microscopic analysis at three excitation wavelengths, 785 nm, 633 nm, and 514 nm. RESULTS: Raman spectra of each of the photoreceptor inner and outer segments (PIS, POS) and of the outer nuclear layer (ONL) of the retina acquired at 785 nm were dominated by vibrational features characteristic of proteins and lipids. There was a clear difference between the inner and outer domains in the spectroscopic regions, amide I and III, known to be sensitive to protein conformation. The spectra recorded with 633 nm excitation mirrored those observed at 785 nm excitation for the amide I region, but with an additional pattern of bands in the spectra of the PIS region, attributed to cytochrome c. The same features were even more enhanced in spectra recorded with 514 nm excitation. A significant nucleotide contribution was observed in the spectra recorded for the ONL at all three excitation wavelengths. A Raman map was constructed of the major spectral components found in the retinal outer segments, as predicted by principal component analysis of the data acquired using 633 nm excitation. Comparison of the Raman map with its histological counterpart revealed a strong correlation between the two images. CONCLUSIONS: It has been demonstrated that Raman spectroscopy offers a unique insight into the biochemical composition of the light-sensing cells of the retina following the application of standard histological protocols. The present study points to the considerable promise of Raman microscopy as a component-specific probe of retinal tissue.

Fe XIII emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for FeXIII are used to generate emission-line ratios involving 3s23p2-3s3p3 and 3s23p2-3s23p3d transitions in the 170-225 and 235-450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new FeXIII emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several FeXIII transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s23p3d1D2 as their upper level. The most useful FeXIII electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in Ne over the range 108-1011cm-3. It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170-225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235-450 Å.

Positronium scattering by atomic hydrogen with inclusion of target excitation channels

Calculations are reported for positronium (Ps) scattering by atomic hydrogen (H) in the energy range 0-6.5 eV in a coupled- pseudostate approximation in which excitation and ionization channels of both the Ps and the H are taken into account. The approximation contains an accurate representation of the van der Waals coefficient. Results are presented for phase shifts, scattering lengths, effective ranges, and various cross sections including partial wave, total, and ortho-para conversion cross sections. An analysis of the possible spin transitions is provided and the energy of the positronium hydride (PsH) bound state is determined. Substantial differences are found from earlier work within the frozen target approximation, now clearly confirming the importance of target excitation channels. Good agreement is obtained with recent calculations of S-wave phase shifts and scattering lengths using the stabilization method. Convergence to the exact binding energy for PsH appears to be slow. Resonances corresponding to unstable states of the positron orbiting H- are seen in the electronic spin singlet partial waves. The importance of the H- formation channel is discussed.

Second-order distorted wave calculation for electron impact ionization of hydrogen

The triple differential cross sections for ionization of atomic hydrogen by electron impact are analysed in the case of coplanar, asymmetric geometry within the framework of second- order distorted wave theory. Detailed calculations are performed without making any approximations (other than numerical) in the evaluation of the second-order amplitude. The present results are compared with experimental measurements and other theoretical calculations for incident energies of 250, 150 and 54.4 eV. It is found that the second-order calculations represent a marked improvement over the results obtained from first-order theories for impact energies of 150 eV and higher. The close agreement between the present second-order plane wave calculation and those of Byron et al calculated using the closure approximation at an incident energy of 250 eV implies that the closure approximation is valid for this energy. The large difference between the present second-order distorted wave calculations and experiment at an incident energy of 54.4 eV suggests that higher order effects are important for incident energies less than 100 eV.

Electron/positron collisions with positronium

We report results for e(+/-)-Ps(Is) scattering in the energy range up to 80 eV calculated in 9-state and 30-state coupled pseudostate approximations. Cross-sections are presented for elastic scattering, ortho-para conversion, discrete excitation, ionization and total scattering. Resonances associated with the Ps(n = 2) threshold are also examined and their positions and widths determined. Very good agreement is obtained with the variational calculations of Ward et al. [J. Phys. B 20 (1987) 127] below 5.1 eV. (C) 2004 Elsevier B.V. All rights reserved.