岩性油藏储层预测的地震理论与方法研究

At present the main object of the exploration and development (E&D) of oil and gas is not the structural oil-gas pools but the subtle lithological oil-gas reservoir. Since the last 90's, the ratio of this kind of pools in newly-added oil reserves is becoming larger and larger, so is the ratio in the eastern oilfields. The third oil-gas resource evaluation indicates the main exploration object of Jiyang depression is the lithological oil-gas pools in future. However, lack of effective methods that are applied to search for this kind of pool makes E&D difficult and the cost high. In view of the urgent demand of E&D, in this paper we deeply study and analyze the theory and application in which the seismic attributes are used to predict and describe lithological oil-gas reservoirs. The great results are obtained by making full use of abundant physics and reservoir information as well as the remarkable lateral continuity involved in seismic data in combination with well logging, drilling-well and geology. ①Based on a great deal of research and different geological features of Shengli oilfield, the great progresses are made some theories and methods of seismic reservoir prediction and description. Three kinds of extrapolation near well seismic wavelet methods-inverse distance interpolation, phase interpolation and pseudo well reflectivity-are improved; particularly, in sparse well area the method of getting pseudo well reflectivity is given by the application of the wavelet theory. The formulae for seismic attributes and coherent volumes are derived theoretically, and the optimal method of seismic attributes and improved algorithms of picking up coherent data volumes are put forward. The method of making sequence analysis on seismic data is put forward and derived in which the wavelet transform is used to analyze not only qualitatively but also quantitatively seismic characteristics of reservoirs.② According to geologic model and seismic forward simulation, from macro to micro, the method of pre- and post-stack data synthetic analysis and application is put forward using seismic in close combination with geology; particularly, based on making full use of post-stack seismic data, "green food"-pre-stack seismic data is as possible as utilized. ③ In this paper, the formative law and distributing characteristic of lithologic oil-gas pools of the Tertiary in Jiyang depression, the knowledge of geological geophysics and the feasibility of all sorts of seismic methods, and the applied knowledge of seismic data and the geophysical mechanism of oil-gas reservoirs are studied. Therefore a series of perfect seismic technique and software are completed that fit to E&D of different categories of lithologic oil-gas reservoirs. ④ This achievement is different from other new seismic methods that are put forward in the recent years, that is multi-wave multi-component seismic, cross hole seismic, vertical seismic, and time-lapse seismic etc. that need the reacquisition of seismic data to predict and describe the oil-gas reservoir. The method in this paper is based on the conventional 2D/3D seismic data, so the cost falls sharply. ⑤ In recent years this technique that predict and describe lithologic oil-gas reservoirs by seismic information has been applied in E&D of lithologic oil-gas reservoirs on glutenite fans in abrupt slop and turbidite fans in front of abrup slop, slump turbidite fans in front of delta, turbidite fans with channel in low slope and channel sanbody, and a encouraging geologic result has been gained. This achievement indicates that the application of seismic information is one of the most effective ways in solving the present problem of E&D. This technique is significant in the application and popularization, and positive on increasing reserves and raising production as well as stable development in Shengli oilfield. And it will be directive to E&D of some similar reservoirs

复杂断块油藏形成机制及剩余油预测

The central uplift in the Huimin depression is famous for its large amounts of faults and small-scale fault-block area, and it is the famed typical complicated fault-block group oil & gas field in the whole world. After many years of rolling exploration and exploitation, many complex oil &gas field have been discovered in the central uplift, and won the splendent fruit. With the gradual deepening and development of the rolling exploitation, the exploration faces more and more difficulties. Therefore, it is important to reveal the forming mechanism and distributing rule of the complex fault-block reservoir, and to realize the forecast of the complex fault-block reservoir, sequentially, expedite the exploration step. This article applies the new multi-subject theory, method and technique such as structure geometry, kinematics, dynamics, structural stress field, fluid potential field, well logging record and constrained inversion of seismic records, coherence analysis, the seal mold and seal history of oil-bounded fault etc, and try to reveal the forming mechanism and distributing law of the complex fault-block reservoir, in result, implements the forecast of the fault-block reservoir and the remaining oil distributing. In order to do so, this article synthetically carries out structural estimate, reservoir estimate, fault sealing history estimate, oil-bearing properties estimate and residual. This article also synthetically researches, describes and forecast the complex fault-block in Huimin depression by use of the techniques, e.g. seismetic data post-stack processing technique, multi-component demarcating technique, elaborate description technique for the fault-block structure, technique of layer forecasting, fault sealing analysis technique, comprehensive estimate technique of fault-block, comprehensive analysis and estimate technique of remaining oil etc. The activities of the faults varies dramatically in the Huimin depression, and most of the second-class and the third-class faults are contemporaneous faults, which control the macroscopical distribution of the reservoir in the Huimin depression. The fourth-class faults cause the complication between the oil & gas among the fault-blocks. The multi-period strong activities of the Linyi fracture resulted in the vertical migration of large amount of oil & gas along with the faults. This is the main reason for the long vertical distribution properties near the Linyi fracture in the Huimin depression. The sealing ability of the fault is controlled by the property,size and direction of the main stress, the contact relationship of the both sides of the fault, the shale polluting factor, and the configuration relationship between the fault move period and the migration period of oil & gas. The article suggest four fault-sealing modes in the research zone for the first time, which establishes the foundation for the further forecast of the complex fault-block reservoir. Numerical simulation of the structural stress field reveals the distribution law and the evolvement progress of the three-period stress field from the end of the Dongying period to the Guantao period to nowadays. This article puts forward that the Linyi and Shanghe regions are the low value of the maximum main stress data. This is combined with the fault sealing history estimate, then multi-forming-reservoir in the central uplift is put forward. In the Shanghe oilfield, the article establishes six reservoir geological modes and three remaining oil distributing modes(the plane, the inside layer and the interlayer), then puts forward six increase production measure to enhance the remaining oil recovery ratio. Inducting the exploitation of oilfield, it wins notable economic effects and social effects.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

Automated Placement of Cameras in a Floorplan to Satisfy Task-Specific Constraints

In many multi-camera vision systems the effect of camera locations on the task-specific quality of service is ignored. Researchers in Computational Geometry have proposed elegant solutions for some sensor location problem classes. Unfortunately, these solutions utilize unrealistic assumptions about the cameras' capabilities that make these algorithms unsuitable for many real-world computer vision applications: unlimited field of view, infinite depth of field, and/or infinite servo precision and speed. In this paper, the general camera placement problem is first defined with assumptions that are more consistent with the capabilities of real-world cameras. The region to be observed by cameras may be volumetric, static or dynamic, and may include holes that are caused, for instance, by columns or furniture in a room that can occlude potential camera views. A subclass of this general problem can be formulated in terms of planar regions that are typical of building floorplans. Given a floorplan to be observed, the problem is then to efficiently compute a camera layout such that certain task-specific constraints are met. A solution to this problem is obtained via binary optimization over a discrete problem space. In preliminary experiments the performance of the resulting system is demonstrated with different real floorplans.

Automated Camera Layout to Satisfy Task-Specific and Floorplan-Specific Coverage Requirements

In many multi-camera vision systems the effect of camera locations on the task-specific quality of service is ignored. Researchers in Computational Geometry have proposed elegant solutions for some sensor location problem classes. Unfortunately, these solutions utilize unrealistic assumptions about the cameras' capabilities that make these algorithms unsuitable for many real-world computer vision applications: unlimited field of view, infinite depth of field, and/or infinite servo precision and speed. In this paper, the general camera placement problem is first defined with assumptions that are more consistent with the capabilities of real-world cameras. The region to be observed by cameras may be volumetric, static or dynamic, and may include holes that are caused, for instance, by columns or furniture in a room that can occlude potential camera views. A subclass of this general problem can be formulated in terms of planar regions that are typical of building floorplans. Given a floorplan to be observed, the problem is then to efficiently compute a camera layout such that certain task-specific constraints are met. A solution to this problem is obtained via binary optimization over a discrete problem space. In experiments the performance of the resulting system is demonstrated with different real floorplans.

Optimisation of Smart Grid performance using centralised and distributed control techniques

A massive change is currently taking place in the manner in which power networks are operated. Traditionally, power networks consisted of large power stations which were controlled from centralised locations. The trend in modern power networks is for generated power to be produced by a diverse array of energy sources which are spread over a large geographical area. As a result, controlling these systems from a centralised controller is impractical. Thus, future power networks will be controlled by a large number of intelligent distributed controllers which must work together to coordinate their actions. The term Smart Grid is the umbrella term used to denote this combination of power systems, artificial intelligence, and communications engineering. This thesis focuses on the application of optimal control techniques to Smart Grids with a focus in particular on iterative distributed MPC. A novel convergence and stability proof for iterative distributed MPC based on the Alternating Direction Method of Multipliers is derived. Distributed and centralised MPC, and an optimised PID controllers' performance are then compared when applied to a highly interconnected, nonlinear, MIMO testbed based on a part of the Nordic power grid. Finally, a novel tuning algorithm is proposed for iterative distributed MPC which simultaneously optimises both the closed loop performance and the communication overhead associated with the desired control.

The effect of hydrogen bonding on the diffusion of water in n-alkanes and n-alcohols measured with a novel single microdroplet method.

While the Stokes-Einstein (SE) equation predicts that the diffusion coefficient of a solute will be inversely proportional to the viscosity of the solvent, this relation is commonly known to fail for solutes, which are the same size or smaller than the solvent. Multiple researchers have reported that for small solutes, the diffusion coefficient is inversely proportional to the viscosity to a fractional power, and that solutes actually diffuse faster than SE predicts. For other solvent systems, attractive solute-solvent interactions, such as hydrogen bonding, are known to retard the diffusion of a solute. Some researchers have interpreted the slower diffusion due to hydrogen bonding as resulting from the effective diffusion of a larger complex of a solute and solvent molecules. We have developed and used a novel micropipette technique, which can form and hold a single microdroplet of water while it dissolves in a diffusion controlled environment into the solvent. This method has been used to examine the diffusion of water in both n-alkanes and n-alcohols. It was found that the polar solute water, diffusing in a solvent with which it cannot hydrogen bond, closely resembles small nonpolar solutes such as xenon and krypton diffusing in n-alkanes, with diffusion coefficients ranging from 12.5x10(-5) cm(2)/s for water in n-pentane to 1.15x10(-5) cm(2)/s for water in hexadecane. Diffusion coefficients were found to be inversely proportional to viscosity to a fractional power, and diffusion coefficients were faster than SE predicts. For water diffusing in a solvent (n-alcohols) with which it can hydrogen bond, diffusion coefficient values ranged from 1.75x10(-5) cm(2)/s in n-methanol to 0.364x10(-5) cm(2)/s in n-octanol, and diffusion was slower than an alkane of corresponding viscosity. We find no evidence for solute-solvent complex diffusion. Rather, it is possible that the small solute water may be retarded by relatively longer residence times (compared to non-H-bonding solvents) as it moves through the liquid.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

Three Hundred Years in Annapolis: Phase III Archaeological Investigations of the Anne Arundel County Courthouse Site (18AP63), Annapolis, Maryland

During the summer of 1994, Archaeology in Annapolis conducted archaeological investigations of the city block bounded by Franklin, South and Cathedral Streets in the city of Annapolis. This Phase III excavation was conducted as a means to identify subsurface cultural resources in the impact area associated with the proposed construction of the Anne Arundel County Courthouse addition. This impact area included both the upper and lower parking lots used by Courthouse employees. Investigations were conducted in the form of mechanical trenching and hand excavated units. Excavations in the upper lot area yielded significant information concerning the interior area of the block. Known as Bellis Court, this series of rowhouses was constructed in the late nineteenth century and was used as rental properties by African-Americans. The dwellings remained until the middle of the twentieth century when they were demolished in preparation for the construction of a Courthouse addition. Portions of the foundation of a house owned by William H. Bellis in the 1870s were also exposed in this area. Construction of this house was begun by William Nicholson around 1730 and completed by Daniel Dulany in 1732/33. It was demolished in 1896 by James Munroe, a Trustee for Bellis. Excavations in the upper lot also revealed the remains of a late seventeenth/early eighteenth century wood-lined cellar, believed to be part of the earliest known structure on Lot 58. After an initially rapid deposition of fill around 1828, this cellar was gradually covered with soil throughout the remainder of the nineteenth century. The fill deposit in the cellar feature yielded a mixed assemblage of artifacts that included sherds of early materials such as North Devon gravel-tempered earthenware, North Devon sgraffito and Northem Italian slipware, along with creamware, pearlware and whiteware. In the lower parking lot, numerous artifacts were recovered from yard scatter associated with the houses that at one time fronted along Cathedral Street and were occupied by African- Americans. An assemblage of late seventeenth century/early eighteenth century materials and several slag deposits from an early forge were recovered from this second area of study. The materials associated with the forge, including portions of a crucible, provided evidence of some of the earliest industry in Annapolis. Investigations in both the upper and lower parking lots added to the knowledge of the changing landscape within the project area, including a prevalence of open space in early periods, a surprising survival of impermanent structures, and a gradual regrading and filling of the block with houses and interior courts. Excavations at the Anne Arundel County Courthouse proved this to be a multi-component site, rich in cultural resources from Annapolis' Early Settlement Period through its Modern Period (as specified by Maryland's Comprehensive Historic Preservation Plan (Weissman 1986)). This report provides detailed interpretations of the archaeological findings of these Phase III investigations.

A comprehensive copper stockpile leach model: background and model formulation

The first phase in the sign, development and implementation of a comprehensive computational model of a copper stockpile leach process is presented. The model accounts for transport phenomena through the stockpile, reaction kinetics for the important mineral species, oxgen and bacterial effects on the leach reactions, plus heat, energy and acid balances for the overall leach process. The paper describes the formulation of the leach process model and its implementation in PHYSICA+, a computational fluid dynamic (CFD) software environment. The model draws on a number of phenomena to represent the competing physical and chemical features active in the process model. The phenomena are essentially represented by a three-phased (solid liquid gas) multi-component transport system; novel algorithms and procedures are required to solve the model equations, including a methodology for dealing with multiple chemical species with different reaction rates in ore represented by multiple particle size fractions. Some initial validation results and application simulations are shown to illustrate the potential of the model.

On the continuum modelling of segregation of granular mixtures during hopper emptying in core flow mode

In this paper, the application of a continuum model is presented, which deals with the discharge of multi-component granular mixtures in core flow mode. The full model description is given (including the constitutive models for the segregation mechanism) and the interactions between particles at the microscopic level are parametrised in order to predict the development of stagnant zone boundaries during core flow discharges. Finally, the model is applied to a real industrial problem and predictions are made for the segregation patterns developed during mixture discharge in core flow mode.

Non-invasive in vivo imaging in small animal research

Non-invasive real time in vivo molecular imaging in small animal models has become the essential bridge between in vitro data and their translation into clinical applications. The tremendous development and technological progress, such as tumour modelling, monitoring of tumour growth and detection of metastasis, has facilitated translational drug development. This has added to our knowledge on carcinogenesis. The modalities that are commonly used include Magnetic Resonance Imaging (MRI), Computed Tomography (CT), Positron Emission Tomography (PET), bioluminescence imaging, fluorescence imaging and multi-modality imaging systems. The ability to obtain multiple images longitudinally provides reliable information whilst reducing animal numbers. As yet there is no one modality that is ideal for all experimental studies. This review outlines the instrumentation available together with corresponding applications reported in the literature with particular emphasis on cancer research. Advantages and limitations to current imaging technology are discussed and the issues concerning small animal care during imaging are highlighted.

Identical extraction behavior and coordination of trivalent or hexavalent f-element cations using ionic liquid and molecular solvents

The extraction of both UO22+ and trivalent lanthanide and actinide ions (Am3+, Nd3+, Eu3+) by dialkylphosphoric or dialkylphosphinic acids from aqueous solutions into the ionic liquid, 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide has been studied and compared to extractions into dodecane. Radiotracer partitioning measurements show comparable patterns of distribution ratios for both the ionic liquid/aqueous and dodecane/aqueous systems, and the limiting slopes at low acidity indicate the partitioning of neutral complexes in both solvent systems. The metal ion coordination environment, elucidated from EXAFS and UV-visible spectroscopy measurements, is equivalent in the ionic liquid and dodecane solutions with coordination of the uranyl cation by two hydrogen-bonded extractant dimers, and of the trivalent cations by three extractant dimers. This is the first definitive report of a system where both the biphasic extraction equilibria and metal coordination environment are the same in an ionic liquid and a molecular organic solvent.

Notes on the Asymmetric Hydrogenation of Methyl Acetoacetate in Neoteric Solvents

Asymmetric hydrogenation of methyl acetoacetate to methyl (R)-3-hydroxybutyrate by [(R)-RuCl(binap)( p-cymen)] Cl has been studied in methanol-ionic liquid and methanol-dense CO(2) solvent systems. The ionic pairs triethylhexylammonium and 1-methylimidazolium with bis(trifluoromethane sulfonyl) imide and hexafluorophosphates were used. The role of ionic pairs on the kinetic parameters and (enantio) selectivity has been demonstrated. Although the CO(2) expanded methanol system suffered from a reduction in both reaction rate and product selectivity, this changed in the presence of water. The high selectivity of the optimized methanol-CO(2)-water-halide system was designed as a consequence of observed additive effects.

Low frequency localized wavepackets in dusty plasmas with opposite charge polarity dust components

The reductive perturbation technique is employed to investigate the modulational instability of dust-acoustic (DA) waves propagating in a four-component dusty plasma. The dusty plasma consists of both positive- and negative-charge dust grains, characterized by a different mass, temperature and density, in addition to a background of Maxwellian electrons and ions. Relying on a multi-fluid plasma model and employing a multiple scales technique, a nonlinear Schrodinger type equation (NLSE) is obtained for the electric potential amplitude perturbation. The occurrence of localized electrostatic wavepackets is shown, in the form of oscillating structures whose modulated envelope is modelled as a soliton (or multi-soliton) solution of the NLSE. The DA wave characteristics, as well as the associated stability thresholds, are studied analytically and numerically. The relevance of these theoretical results with dusty plasmas observed in cosmic and laboratory environments is analysed in detail, by considering realistic multi-component plasma configurations observed in the polar mesosphere, as well as in laboratory experiments.