层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

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

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

板北断块油藏储层均质特征与剩余油分布研究

With the growing development and perfection of reservoir describing technology, its research achievements have played an increasingly important role in old oilfields in recent years. Reservoir description quantitatively describes, characterizes and predicts every kind of reservoir characters in 3D space. The paper takes Banbei block reservoir as an object, studies the reservoir characters and residual oil distributing characteristics of gravity flow genetic reservoir, and definitudes potential adjustment direction of reservoir development. Main achievements are gained as follows. Through fine correlation of strati graphic sequence, the classification of layers and single sands of main payzones in Banbei block is ascertained, the classifying methods of sedimentary unit in gravity flow reservoir characterized with picked cyclical marker bed are formed. On the basis of comprehensive logging evaluation, depositional characters of Banbei block are studied, and classifying methods of sedimentary microfacies in gravity flow reservoir are described. The sedimentary background of main oil layers in Banbei block is open lake with shallow water, and belongs to lacustrine underwater gravity flow- lacustrine phase depositional system. Main microfacies types are underwater water course^ water course side-wing, underwater floodplain, between two water courses, and lacustrine mud, etc. Reservoir sands mainly are underwater water course sands. Influenced by distributing characters of gravity flow underwater water course, sand shapes in plane mainly are stripe, finger-shape, tongue-shape. Sand distribution shows obvious split property. Sands overlap each other. According to comprehensive analysis of lithologic data, logging parameters, and dynamic production data, the researching threads and methods of reservoir heterogeneous characters are perfected. The depositional characters of gravity flow underwater water course in Banbei block determine its high reservoir heterogeneity. Macroscopic heterogeneity is studied in many aspects such as the scale of layers, the scale of single sands, in-situ scale, the distribution of interlayer types, the interlayer scale, and heterogeneity in plane. Thus, heterogeneous characters of reservoir are thoroughly analyzed. Through microscopic research of reservoir, the types of porous structure and related parameters are determined. According to the analysis of dynamic production data, the reaction and inner influential factors of reservoir heterogeneity in waterflood development are further revealed. Started with the concept and classifying methods of flow unit, clustering classification which can better meet the requirements of production is formed. The flow unit of Banbei block can be classified into four types. According to comprehensive evaluation, the first and second type of flow unit have better percolating capability and reserving capability. Research thread of 3D model-building and reservoir numerical simulation combined as an integral is adopted. The types and characters of residual oil distribution are determined. Residual oil of Banbei block mainly distributes in the boundary of sands, near the faults, areas with non-perfect injection-production well pattern , undeveloped sands, vertically poor developed layers. On the basis of comprehensive reservoir study, the threads and methods of improving development effect towards reservoir with high water cut, high recovery percent, serious heterogeneity are ascertained. The whole waterflood development effect of Banbei block reservoir is good. Although its water cut and recovery percent is relatively high, there is still some potential to develop. According to depositional characters of gravity flow and actual production situation? effective means of further improving development level are as follows. We should drill new wells in every kind of areas abounding with residual oil, implement comprehensive measures such as increasing liquid discharge, cyclic waterflood, changing fluid direction when injection-production well pattern is perfected, improve water quality, enhance displacement efficiency in flooding.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Experimental Study on Vortex-Induced Vibrations of a Pipeline Near a Plane Boundary in Steady Flow

对单向水流作用下近壁管道横向涡激振动进行了实验模拟,重点探讨了管道与壁面间隙比(e/D)对管道涡激振动幅值和涡激振动频率响应特性的影响规律.实验结果表明,管道与壁面间隙宽度对管道涡激振动特性有较明显影响.在较大间隙比(e/D＞0.66)下,管道振幅随着Vr数的增大先快速增长到最大值,然后平缓下降;在振动初期(即Vr数较小时),管道振动频率变化基本符合Strouhal规律;在振动中后期(即Vr数较大时),管道振动频率变化不符合Strouhal规律,而在管道固有频率附近缓慢增长.在较小间隙比(e/D＜0.30)下,管道振幅随Vr数的增大先平缓上升到最大值,随后较快速下降;在振动初期,管道振动频率变化不遵循Strouhal规律;在整个振动范围内,与较大间隙比情况相比,随着Vr数增加,管道振动频率增长幅度明显较大.

Mathematical Modeling of Near-Wall Flows of Two-Phase Mixture with Evaporating Droplets

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream. 在双连续介质理论框架下,采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程,给出了控制汽雾流的相似判据。对于沿曲面的流动,边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体验点附近蒸汽。液滴边界层的数值计算结果。它们表明：由于蒸发,在边界层内近壁处形成了一个无液滴区域；在该区上边界处,液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加,甚至在自由来流中液滴质量浓度很低时此效应依然存在。

Stability of ZrTiCuNiBe Bulk Metallic Glass Upon Isothermal Annealing Near the Glass Transition Temperature

The stability of Zr41Ti14Cu12.5Ni10Be22.5 bulk metallic glass (BMG) upon isothermal annealing near the glass transition temperature has been investigated by using x-ray diffraction, differential scanning calorimetry, and the pulse echo overlap method. The density, elastic constants, and thermodynamic parameters as well as their annealing time dependence have been determined. The microstructural and properties changes of the annealed BMG were checked by acoustic measurement. Obvious structural and property changes were observed with prolonged annealing of the BMG near the glass transition temperature.

Experimental Study on Vortex-Induced Vibrations of Submarine Pipeline Near Seabed Boundary in Ocean Currents

Unlike most previous studies on vortex- induced vibrations of a cylinder far from a boundary, this paper focuses On the influences of close proximity of a submarine pipeline to a rigid seabed boundary upon the dynamic responses of the pipeline in ocean currents. The effects of gap-to-diameter ratio and those of the stability parameter on the amplitude and frequency responses of a pipeline are investigated experimentally with a novel hydro-elastic facility. A comparison is made between the present experimental results Of the amplitude and frequency responses for the pipes with seabed boundary effects and those for wall-free cylinders given by Govardhan and Williamson (2000) and Anand ( 1985). The comparison shows that the close proximity of a pipeline to seabed has much influence on the vortex- induced vibrations of the pipeline. Both the width of the lock-in ranges in terms of V, and the dimensionless amplitude ratio A(max)/D become larger with the decrease of the gap-to-diameter ratio e/D. Moreover, the vibration of the pipeline becomes easier to occur and its amplitude response becomes more intensive with the decrease of the stability parameter, while the pipeline frequency responses are affected slightly by the stability parameter.

An In-Situ Observation On Initial Aggregation Process Of Colloidal Particles Near Three-Phase Contact Line Of Air, Water And Vertical Substrate

The self-assembling process near the three-phase contact line of air, water and vertical substrate is widely used to produce various kinds of nanostructured materials and devices. We perform an in-situ observation on the self-assembling process in the vicinity of the three phase contact line. Three kinds of aggregations, i.e. particle-particle aggregation, particle-chain aggregation and chain-chain aggregation, in the initial stage of vertical deposition process are revealed by our experiments. It is found that the particle particle aggregation and the particle-chain aggregation can be qualitatively explained by the theory of the capillary immersion force and mirror image force, while the chain-chain aggregation leaves an opening question for the further studies. The present study may provide more deep insight into the self-assembling process of colloidal particles.

Experimental Study Of Vortex-Induced Vibrations Of A Pipeline Near An Erodible Sandy Seabed

Based on similarity analyses, a series of experiments have been conducted with a newly established hydro-elastic facility to investigate the transverse vortex-induced vibrations (VIVs) of a submarine pipeline near an erodible sandy seabed under the influence of ocean currents. Typical characteristics of coupling processes between pipe vibration and soil scour in the currents have been summarized for Case 1: pipe is laid above seabed and Case 11: pipe is partially embedded in seabed on the basis of the experimental observations. Pipe vibration and the corresponding local scour are usually two coupled physical processes leading to an equilibrium state. The influence of initial gap-to-diameter ratio (e(0)/D) on the interaction between pipe vibration and local scour has been studied. Experimental results show that the critical values of V-r for the initiation of VIVs of the pipe near an erodible sand bed get bigger with decreasing initial gap-to-diameter ratio within the examined range of e(0)/D (-0.25 < e(0)/D < 0.75). The comparison of the pipe vibrations near an erodible soil with those near a rigid boundary and under wall-free conditions indicates that the vibration amplitudes of the pipe near an erodible sand bed are close to the curve fit under wall-free conditions; nevertheless, for the same stability parameter, the maximum amplitudes for the VIV coupled with local scour increase with the increase of initial gap-to-diameter ratio. (c) 2007 Elsevier Ltd. All rights reserved.

Transition waves and nonlinear interactions in the near wake of a circular cylinder

Transition waves and interactions between two kinds of instability-vortex shedding and transition wave in the near wake of a circular cylinder in the Reynolds number range 3 000-10 000 are studied by a domain decomposition hybrid numerical method. Based on high resolution power spectral analyses for velocity new results on the Reynolds-number dependence of the transition wave frequency, i.e. f(t)/f(s) similar to Re-0.87 are obtained. The new predictions are in good agreement with the experimental results of Wei and Smith but different from Braza's prediction and some early experimental results f(t)/f(s) similar to Re-0.5 given by Bloor et nl. The multi-interactions between two kinds of vortex are clearly visualized numerically. The strong nonlinear interactions between the two independent frequencies (f(t), f(s)) leading to spectra broadening to form the coupling mf(s) +/- nf(t) are predicted and analyzed numerically, and the characteristics of the transition are described. Longitudinal variations of the transition wave and its coupling are reported. Detailed mechanism of the flow transition in the near wake before occurrence of the three-dimensional evolution is provided.

Character of simplified Navier-Stokes (SNS) equations near separation point for two-dimensional laminar flow over a flat plate

It is proved that the simplified Navier-Stokes (SNS) equations presented by Gao Zhi[1], Davis and Golowachof-Kuzbmin-Popof (GKP)[3] are respectively regular and singular near a separation point for a two-dimensional laminar flow over a flat plate. The order of the algebraic singularity of Davis and GKP equation[2,3] near the separation point is indicated. A comparison among the classical boundary layer (CBL) equations, Davis and GKP equations, Gao Zhi equations and the complete Navier-Stokes (NS) equations near the separation point is given.

Effect of grain size on slow fatigue crack propagation and plastic deformation near crack tip

Fatigue crack growth and its threshold are investigated at a stress ratio of 0.5 for the three-point bend specimen made of Austenitic stainless steel. The effect of grain size on the crack tip plastic deformation is investigated. The results show that the threshold value Δkth increases linearly with the square root of grain size d and the growth rate is slower for materials with larger grain size. The plastic zone size and ratio for different grain sizes are different at the threshold. The maximum stress intensity factor is kmax and σys is the yield strength. At the same time, the characteristics of the plastic deformation development is discontinuous and anti-symmetric as the growth rate is increased from 2·10—8 to 10−7 mm/cycle.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.