鄂尔多斯盆地上古生界低渗透砂岩有效储层识别与预测

Sulige Gasfield, with a basically proven reserve as high as one trillion cubic meters, is one giant gas field discovered in China. The major gas -bearing layers are Upper Paleozoic strata with fluvial-lacustrine sedimentary facies. Generally, gas reservoirs in this field are characteristic by "five low" properties, namely low porosity, low permeability, low formation pressure, low productivity and low gas abundance. Reservoirs in this field also feature in a large distribution area, thin single sandbody thickness, poor reservoir physical properties, thin effective reservoir thickness, sharp horizontal and/or vertical changes in reservoir properties as well as poor connectivity between different reservoirs. Although outstanding achievements have been acquired in this field, there are still several problems in the evaluation and development of the reservoirs, such as: the relation between seismic attributes and reservoir property parameters is not exclusive, which yields more than one solution in using seismic attributes to predict reservoir parameters; the wave impedance distribution ranges of sandstone and mudstone are overlapped, means it is impossible to distinguish them through the application of post-stack impedance inversion; studies on seismic petrophysics, reservoir geophysical properties, wave reflection models and AVO features have a poor foundation, makes it difficult to recognize the specific differences between tight sandstone and gas-bearing sandstone and their distribution laws. These are the main reasons causing the low well drilling success rate and poor economic returns, which usually result in ineffective development and utilization of the field. Therefore, it is of great importance to perform studies on identification and prediction of effective reservoirs in low permeable sandstone strata. Taking the 2D and 3D multiwave-multicomponent seismic exploration block in Su6-Su5 area of Sulige field as a study area and He 8 member as target bed, analysis of the target bed sedimentary characteristics and logging data properties are performed, while criteria to identify effective reservoirs are determined. Then, techniques and technologies such as pre-stack seismic information (AVO, elastic impedance, wave-let absorption attenuation) and Gamma inversion, reservoir litological and geophysical properties prediction are used to increase the precision in identifying and predicting effective reservoirs; while P-wave and S-wave impedance, ratio of P/S wave velocities, rock elastic parameters and elastic impedance are used to perform sandstone gas-bearing property identification and gas reservoir thickness prediction. Innovative achievements are summarized as follows: 1. The study of this thesis is the first time that multiwave-multicomponent seismic data are used to identify and predict non-marine classic reservoirs in China. Through the application of multiwave-multicomponents seismic data and integration of both pre-stack and post-stack seismic data, a set of workflows and methods to perform high-precision prediction of effective reservoirs in low permeable sandstone is established systematically. 2. Four key techniques to perform effective reservoir prediction including AVO analysis, pre-stack elastic wave impedance inversion, elastic parameters inversion, and absorption attenuation analysis are developed, utilizing pre-stack seismic data to the utmost and increasing the correct rate for effective reservoir prediction to 83% from the former 67% with routine methods. 3. This thesis summarizes techniques and technologies used in the identification reservoir gas-bearing properties using multiwave-multicomponent seismic data. And for the first time, quantitative analysis on reservoir fluids such as oil, gas, and/or water are carried out, and characteristic lithology prediction techniques through the integration of pre-stack and post-stack seismic prediction techniques, common seismic inversion and rock elastic parameters inversion, as well as P-wave inversion and converted wave inversion is put forward, further increasing the correct rate of effective reservoir prediction in this area to 90%. 4. Ten seismic attribute parameters are selected in the 3D multi-wave area to perform a comprehensive evaluation on effective reservoirs using weighted-factor method. The results show that the first class effective reservoir covers an area of 10.08% of the study area, while the second and the third class reservoirs take 43.8% and 46% respectively, sharply increasing the success rate for appraisal and development wells.

Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

Influence of Electrostatic Interaction on Fibrinogen Adsorption on Gold Studied by Imaging Ellipsometry Combined with Electrochemical Methods

Imaging ellipsometry was combined with electrochemical methods for studying electrostatic interactions of protein and solid surfaces. The potential of zero charge for gold-coated silicon wafer/solution interfaces wad determined by AC impedance method. The potential of the gold-coated silicon wafer was controlled at the potential of zero charge, and the adsorption of fibrinogen on the potential-controlled and non-controlled surfaces was measured in real time at the same time by imaging ellipsometry The effect of electrostatic interaction was studied by comparing the difference between the potential of controlled adsorption and the Potential of noncontrolled adsorption. It was shown that the rate of fibrinogen adsorption on the potentiostatic surface was faster than that on the nonpotentiostatic surface. The electrostatic influence on fibrinogen adsorption on the gold-coated silicon wafer was weak, so the hydrophobic interaction should be the major affinity.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

Assessment of Microscale Test Methods of Peeling and Splitting along Surface of Thin-Film/Substrates

Peel test methods are assessed through being applied to a peeling analysis of the ductile film/ceramic substrate system. Through computing the fracture work of the system using the either beam bend model (BB model) or the general plane analysis model (GPA model), surprisingly, a big difference between both model results is found. Although the BB model can capture the plastic dissipation phenomenon for the ductile film case as the GPA model can, it is much sensitive to the choice of the peeling criterion parameters, and it overestimates the plastic bending effect unable to capture crack tip constraint plasticity. In view of the difficulty of measuring interfacial toughness using peel test method when film is the ductile material, a new test method, split test, is recommended and analyzed using the GPA model. The prediction is applied to a wedge-loaded experiment for Al-alloy double-cantilever beam in literature.

On the Gas-Liquid Two-Phase Flow Regimes and their Recognition Methods

气液两相流体系是一个复杂的多变量随机过程体系,流型的定义、流型过渡准则和判别方法等方面的研究是多相流学科目前研究的重点内容。本文就与气液两相流流型及其判别有关的研究状况进行了回顾和评述,力图反映近年来气液两相流流型及其判别问题研究的状态和趋势。

Practical Criterion for Universal Laws and Practical Methods for Construction of Co-variance Sets

给出相对论力学中普遍定律的实用判别法和协变集的实用构造法,还给出实现非普遍定律的“可导出性”的一种实用方法．

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic–plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a “hold” period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and “hold-at-the-maximum-indenter-displacement” for determining the instantaneous modulus using spherical indenters.

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.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic-plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a "hold" period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and "hold-at-the-maximum-indenter-displacement" for determining the instantaneous modulus using spherical indenters.

Pipeline Corrosion Analyzed by Fuzzy Reliability Methods

Fuzzy reliability methods are used to study the corrosion of pipelines. Three methods are used. They consist of using fracture failure modes, failure assessment diagram (FAD) and residual strength for establishing fuzzy reliability. Calculations are made by application of JC, improved GA-JC and Mente-carlo methods. Examples for oilfield injecting water pipeline show the residual strength well agree with field data. Mente-carlo methods appear to yield results that have better agreement with field data.

Analytical Methods of Residual Strength in Injecting Water Pipeline

In this paper, common criterions about residual strength evaluation at home and abroad are generalized and seven methods are acquired, namely ASME-B31G, DM, Wes-2805-97, CVDA-84, Burdekin, Irwin and J integral methods. BP neural network are Combined with Genetic Algorithm (GA) named by modified BP-GA methods to successfully predict residual strength and critical pressure of injecting water, corrosion pipelines. Examples are shown that calculation results of every kind of method have great difference and calculating values of Wes-2805-97 criterion, ASME-B31G criterion, CVDA-84 criterion and Irwin fracture mechanics model are conservative and higher than, those of J integral methods while calculating values of Burdiken model and DM fracture mechanics model are dangerous and less than those of J integral methods and calculating values of modified BP-GA methods are close and moderate to those of J integral methods. Therefore modified BP-GA methods and J integral methods are considered better methods to calculate residual strength and critical pressure of injecting water corrosion pipelines

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.