塔里木盆地寒武－奥陶系白云岩类型、成因机理及其孔隙性研究

Cambrian-Ordovician dolostones in Tarim Basin are hydrocarbon reservoir rocks of vital importance. Under the guidance of the theories of sedimentology and the sedimentology of carbonate reservoir, based on the first-hand qualitative and quantitative data especially, combined with micro-study, geochemical and reservoir capacity analysis, and precursor research, the origin and reservoir characteristics of the dolostones were discussed. Based on detailed petrographic investigations, four types of dolostone have been recognized, which are, respecitively, mud-silt-sized dolostones, algal laminated dolostones (ALD), prophyritic dolostone, and neomorphic dolostone. Mud-silt-sized dolostones always presents as laminas together with evaporated signatures, its REE patterns and ΣREE are all close to that of the finely crystalline limestone. This kind of dolomite probably experienced relatively low fluid-rock ratio during diagenesis was formed in hypersaline and oxidizing environment and involved fast dolomitization process. It was dolomitized by evaporated seawater in sabkha environment.The main primary fabrics of algal lamination in algal laminated dolomite (ALD) can still be identified and its ΣREE (21.37) is very close to that of algae. This reveals that ALD was dolomitized during early diagenesis and algae possibly played an important role. The ALD was formed under mediation of organic matter and Mg2+ were supplied by magnesium concentrated algal laminites and sea water. Prophyritic dolostones presents mainly as patchy occurrence and yield the lowest δ13C and Z value. Its ΣREE is much less than that of the finely crystalline limestone. These characteristics reveal that the cloudy cores were dolomitized in shallow early diagenetic environments by pore fluids riched in Mg2+. Whereas the clear rims were likely formed in subsequent burial into deeper subsurface environments, and the Mg2+ needed for further dolomitization possibly was supplied by the transformations of clay minerals. Neomorphic dolostones consist of coarse, turbid crystals and exhibits sucrosic and mosaic textures. It has highest Fe2+ contents and average homogeneous temperature (110.2℃). Collectively, these characteristics demonstrate that the neomorphic dolostones was likely formed by recrystallization of pre-existing dolomites during deep burial. The ΣREE of the four types of dolostone distinctly differentiates from each other. However, their REE patterns are all enriched in LREE, depleted in HREE and have Eu negative anomalies. Its ΣREE 13.64 ppm, less than 1/4 of finely crystalline limestone, and ranks the lowest in the 4 types.These characteristics are comparable to those of finely crystalline limestone, and are mainly infuenced by the sea water. These four types of dolostone show similar REE mobility behaviour and no significant fractionation, althouth they have been subjected to evidently different diageneses. Seven main pore types are identified in the dolostones , which are fenestral, moldic, intercrystal, dissolved，breccia, dissolved breccia and stylolite pores. Fenestral pores are primary and the others are secondary. The dissolved pores and intercrystal pores are the most important reservoir spaces and followed by breccias and dissolved breccia pores, and the moldic and fenestral pores are less important. Stylolites can enhance permeability of reservoir rocks in one hand, for the other hand, the capacity of reservoir and permeability are enhanced and then better reservoir rocks can be formed when they are combined with patchy dolostones. The relationship between porosity and the type of dolostones is that the dissolved neomorphic dolostones have the highest porosity of 3.65%, than followed by dissolved Mud-silt-sized dolostones of 3.35%. The mud-silt-sized dolostones without dissolution have the lowest porosity of 0.90%. Moreover, the porosity of prophyritic dolostones and the neomorphic dolostones without dissolution are lower, respectively 1.675% and 1.41%. Although algal laminated dolostones consist of euhedral crystals and riched in intercrystal pores, its porosity just yields 1.20%. The relationship between permeability and the type of dolostones is that that algal laminated dolostones have the highest permeability of 0.462mD and followed by 0.065mD of prophyritic dolostones. Dissolution have no significant influence on the permeability of neomorphic dolostones and this presented by the permeability of dissolved and non-dissolved are very close, respectively 0.043mD and 0.062mD. No matter dissolved or not, mud-silt-sized dolostones are much less permeable. The permeability of non-dissolved and dissolved are 0.051mD and 0.016mD. Collectively, in the study area, neomorphic dolostones can form high quality reservoir.

各向同性与各向异性介质中地震波旅行时正演模拟与反演成像

By seismic tomography, interesting results have been achieved not only in the research of the geosphere with a large scale but also in the exploration of resources and projects with a small scale since 80'. Compared with traditional inversion methods, seismic tomography can offer more and detailed information about subsurface and has been being paid attention by more and more geophysicists. Since inversion based on forward modeling, we have studied and improved the methods to calculate seismic traveltimes and raypaths in isotropic and anisotropic media, and applied the improved forward methods to traveltime tomography. There are three main kinds of methods to calculate seismic traveltime field and its ray path distribution, which are ray-tracing theory, eikonal equation by the finite-difference and minimum traveltime tree algorithm. In ray tracing, five methods are introduced in the paper, including analytic ray tracing, ray shooting, ray bending, grid ray tracing and rectangle grid ray perturbation with three points. Finite-difference solution of eikonal equation is very efficient in calculation of seismic first-break, but is awkward in calculation of reflection traveltimes. We have put forward a idea to calculate traveltimes of reflected waves using a combining way of eikonal equation method and other one in order to improve its capability of dealing with reflection waves. The minimum traveltime tree algorithm has been studied with emphases. Three improved algorithms are put forward on the basis of basic algorithm of the minimum traveltime tree. The first improved algorithm is called raypath tracing backward minimum traveltime algorithm, in which not only wavelets from the current source but also wavelets from upper source points are all calculated. The algorithm can obviously improve the speed of calculating traveltimes and raypaths in layered or blocked homogeneous media and keep good accuracy. The second improved algorithm is raypath key point minimum traveltime algorithm in which traveltimes and raypaths are calculated with a view of key points of raypaths (key points of raypths mean the pivotal points which determine raypaths). The raypath key point method is developed on the basis of the first improved algorithm, and has better applicability. For example, it is very efficient even for inhomogeneous media. Another improved algorithm, double grid minimum traveltime tree algorithm, bases upon raypath key point scheme, in which a model is divided with two kinds of grids so that the unnecessary calculation can be left out. Violent undulation of curved interface often results in the phenomenon that there are no reflection points on some parts of interfaces where there should be. One efficacious scheme that curved interfaces are divided into segments, and these segments are treated respectively is presented to solve the problem. In addition, the approximation to interfaces with discrete grids leads to large errors in calculation of traveltimes and raypaths. Noting the point, we have thought a new method to remove the negative effect of mesh and to improve calculation accuracy by correcting the traveltimes with a little of additional calculation, and obtained better results.

利用重力与地形数据探讨我国大陆若干地区岩石圈有效弹性厚度

The lithosphere behaves as an elastic thin plate overlying a weak fluid. When loaded by topography, it deflects to compensate the topography and isostasy is achieved. The response of the plate is often characterized by the flexural rigidity, or equivalently, by the effective elastic thickness(EET). The relationships between gravity and topography have commonly been employed to investigate quantitively the isostatic compensation of the lithosphere. Cross-spectral techniques have been developed to estimate the admittance and coherence between gravity and topography. The observed admittance and coherence functions, when compared with theoretical admittance and coherence computed for various models and parameters, can provide an estimate of the effective elastic thickness. As for coherence, the wavelength at which the coherence drops from approaching unity at long wavelengths to value approaching zero at short wavelengths is a measure of the rigidity of the lithosphere. This research takes advantage of the high-resolution gravity and topography data in Erdos. Using the coherence technique, we have estimated the effective elastic thickness of this region. Subsurface loads are also considered in our calculation. When the admittance between topography and gravity is obtained, we can filter the topography to give an estimate of gravity anomaly. This would be very important when the gravity data is scarce. Several other regions of China have been selected to investigate the effective elastic thickness of lithosphere. We compare those result with lithosphere thickness obtained through seismological technique and heat flow of the region. We find the effective elastic thickness is always smaller than the seismogenic thickness. This is due to what the effective elastic thickness reflects is only the up elastic part of the lithosphere. And we also find there is some degree of correlation between effective elastic thickness and heat flow. This suggests EET is probably controlled by the thermal state of the lithosphere and correlates with their tectonic age. Thus the estimations of the effective elastic thickness of lithosphere can help to investigate geophysical features of active tectonics of continental lithosphere.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

松辽盆地河流-三角洲体系高分辩率层序地层学、储层构形及非均质模型研究

This dissertation tries to combine the new theories of high-resolution sequence stratigraphy and reservoir architecture with fine sedimentology to form a integral theory system -"high-resolution sequence sedimentology", which can be applied widely ranging from the early petroleum exploration to the tertiary recovery stage in marine and terrestrial basin. So the west slope area in south of Songliao basin, in which, early-fine exploration have been developed, and Xingnan area of Daqing placantictine in high water-bearing and tertiary recovery stage, are selected as target areas to research and analyze. By applying high-resolution sequence stratigraphy theory as well as analysis of source area-facies, the west slope area has been divided into two source areas and three drainage systems and the following conclusions have been drawn: three high values sandstone areas, two sandstone pinchout zones and one stratigraphic pinchout overlap; the facies between Baicheng-Tongyu drainage system is frist ascertained as large-scale argillaceous filled plain facies; fine-grained braided channel-delta depositional system has been found; plane sedimentary facies and microfacies maps of different-scale sequence have been completed, and then twenty-eight lithologic traps have been detected in the east of Taobao-Zhenlai reverse fault zone; In no exploration area of the west, large-scale stratigraiphic overlap heavy oil reservoirs has first been found, which has become an important prograss. In Xingnan area, in the view of high-resolution sequence stratigraphy, the surface of unconformity (the bottom of SSC13) in P I group has been identified, and the following method and technique have been advanced: the division and correlation methods of short-scale base-level cycle sequence (SSC); the comprehensive research methods of SSC plane microfacies; the division technique of hierarchy and type of flow unit, the origin of large-scale composite sandbody and flow unit; And ,on the basis of these, 103 monosandstone bodies and 87 flow units of the third levels have been identified, and four levels of flow units model of five sandstone-bodies types have been established. Because it is a very difficult task all over the world to research architecture in subsurface monosandstone body, brings forward a series of techniques as follows: technique of researching architecture of thin interbed in subsurface monosandstone body; classification, type and liquid-resisting mechanism of thin interbed; multiple-remember vertical subsequence model of remaining oil in monosadstone body. Models of heterogeneity and architecture of thin interbed in five types of monosandstone body have been established. Applying these techniques, type and distribution of remaining oil in different types of monosandstone bodies have been predicated.

油气藏自然电位异常与油气赋存的相关性研究

The Stack Spontaneous Potential (SSP) is a direct hydrocarbon location technology and a new hydrocarbon detection method with independent intellectual property. A subsurface hydrocarbon accumulation associated with the upward hydrocarbon micro seepage induces a relatively strong negative potential abnormal zone, of which the anomaly can be measured on the surface with specially designed instruments through careful field measuring procedures. With special software programmed according to a unique geochemical and geophysical model, the original data are analyzed, processed and interpreted on the computer, and then on a series of resulting anomaly distribution maps and/or profiles, the favorable surface locations of the hydrocarbon accumulations can be easily identified. The study of the SSP has been started since 1989, and especially from 1996 to 1997, both profile and area tests were conducted in the Daqing Oilfield. On the testing line of 15kms, there are 6 wells in total, among which some are oil-producing wells, and some are water-producing wells. The final matching ratio of the favorable oil well locations and the possible water well locations predicted by the SSP to those of known wells was up to 83 percent. In the area test, of which the acreage is 800 km2, the matching ratio compared with the existing wells was 87 percent; furthermore, regarding to wells subsequently drilled after the test, the matching ratio was 85 percent. The matching ratio in the development area is more about 10 percent than those of in exploration area. The reason is that, comparing the exploration area, the development area acreage is less and the container rocks are more simplex. In development area there is not so much interference of SSP also. Since 1997 the SSP has been tested and applied all over China to a number of hydrocarbon bearing basins and known oil fields, including the Daqing, Jiangsu, Changqing, Shengli, Nanyang, Jianghan and Zhongyuan Oilfields, only to name a few. The SSP surveys in total areas of over 10,000km2 in more than 30 regions in China so far have been completed in various exploration and development stages, the satisfactory outcomes of which have further evidenced that the dependence between the SP anomaly and abundance of hydrocarbon. Up to date, a substantial amount of successful tests and actual surveys finished in exploration and development practices have evidenced that the SSP is significantly more reliable in comparison with any other similar direct hydrocarbon indication technique generally known to the oil industry, such as the Redox. The SSP can be applied to search for almost all kinds of hydrocarbon accumulations, regardless of the type of traps, such as structural, stratigraphic, buried hill traps, and so on; however, it is interesting to be noted that the SSP seems to be particularly effective in detecting the stratigraphic oil traps according to our practices. On the other hand, there is virtually no surface geographical constrains in terms of field data acquisition, except for those water covered areas, because of the inherent characteristics of the technology itself. Furthermore, utilizing the SSP requires no special considerations to subsurface geological conditions in regard to formation resistivity, since the SSP measurements will not be influenced by either overly high or overly low resistivity of formations lying above the hydrocarbon accumulations. There are two kind of theory, of which, as we know one is called hypbyssal theory such as "Redox"[61 the other is call plutonic theory such as cracking of hydrocarbon [8][9] and natural polarization [3], to describe the mechanism of SP anomaly of oil reservoir and to indicate that the dependence between the SP abnormality and abundance of hydrocarbon has be existed theoretically/The quantitative dependence, which has not been founded due to the complicity of container rocks, be discovered during the exploration and development practices is the crux to the quantitative analysis of SP Anomaly processing. Based on the thorough study of the complex of collector rocks, every kind of thickness of collector rock can be conversed to be a standard effective thickness; the thickness is called apparent effective thickness (AET). The conversation coefficient (ai, 1=1,2,3) could be determined by the variety of every collector rock storability (CRS). The discoveration of quantitative: dependence between AET and the amplitude of SSP, in the practices of exploration and development, is a promotion for the SSP supplied in the oil exploration, and make the data analysis forward to the quantitative stage.

VSP克希霍夫法多波联合成像

Geophones being inside the well, VSP can record upgoing and downgoing P waves, upgoing and downgoing S waves simultaneously.Aiming at overcoming the shortages of the known VSP velocity tomography , attenuation tomography , inverse Q filtering and VSP image method , this article mainly do the following jobs:CD; I do the common-source-point raytracing by soving the raytracing equations with Runge-Kutta method, which can provide traveltime , raypath and amplitude for VSP velocity tomography , attenuation tomography and VSP multiwave migration.(D. The velocity distribution can be inversed from the difference between the computed traveltime and the observed traveltime of the VSP downgoing waves. I put forward two methods: A. VSP building-velocity tomography method that doesn't lie on the layered model from which we can derive the slowness of the grids' crunodes . B. deformable layer tomography method from which we can get the location of the interface if the layer's velocity is known..(3). On the basis of the velocity tomography , using the attenuation information shown by the VSP seismic wave , we can derive the attenuation distribution of the subsurface. I also present an algorithm to solve the inverse Q filtering problem directly and accurately from the Q modeling equation . Numerical results presented have shown that our algorithm gives reliable results . ?. According to the theory that the transformed point is the point where the four kinds of wave come into being , and where the stacked energy will be the largest than at other points . This article presents a VSP multiwave Kirchhoff migration method . Application on synthetic examples and field seismic records have shown that the algorithm gives reliable results . (5). When the location of the interface is determined and the velocity of the P wave and S wave is known , we can obtain the transmittivity and reflection coefficient 5 thereby we can gain the elastic parameters . This method is also put into use derive good result.Above all, application on models and field seismic records show that the method mentioned above is efficient and accurate .

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

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.

地震反射界面形态与介质各向异性属性联立成像：地震反射波旅行时与偏振信息综合反演

It has been a difficult problem faced by seismologists for long time that how exactly to reconstruct the earth's geometric structure and distribution of physical attributes according to seismic wave's kinematical and dynamic characteristics, obtained in seismological observation. The jointing imaging of seismic reflector and anisotropy attributes in the earth interior is becoming the research hot spot. The limitation of shoot and observation system makes that the obtained seismic data are too scarce to exactly reconstruct the geological objects. It is popular that utilizing only seismic reflection traveltimes or polarizations information make inversion of the earth's velocity distribution by fixing seismic reflector configuration (vice versa), these will lead to the serious non-uniqueness reconstruction due to short of effective data, the non-uniqueness problem of reconstructing anisotropy attributes will be more serious than in isotropy media. Obviously it is not enough to restrict the media structure only by information of seismic reflection traveltimes or polarizations, which even sometimes will lead to distorted images and misinterpretation of subsurface structure. So we try to rebuild seismic reflection structure (geometry) and media anisotropic structure (physics) in the earth interior by jointing data of seismic wave kinematics and dynamics characteristics, we carry out the new experiment step by step, and the research mainly comprises of two parts: one is the reconstruction of P-wave vertical velocity and anisotropic structure(Thomsen parameter s and 8) in the transversely isotropic media with vertical symmetrical axis(VTI) by fixing geometrical structure, and the other is the simultaneous inversion of the reflector surface conformation and seismic anisotropic structure by jointing seismic reflection traveltimes and polarizations data. Simulated annealing method is used to the first research part, linear inversion based on BG theory and Simulated annealing are applied to the second one. All the research methods are checked by model experiments, then applied to the real data of the wide-angle seismic profile from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. The results are as following The inversion results based on jointing seismic PP-wave or PSV-wavereflection traveltimes and polarizations data are more close to real model than themodels based simply on one of the two data respectively. It is shown that the methodwe present here can effectively reconstruct the anisotropy attributes in the earth'sinterior when seismic reflector structure is fixed.The layer thickness, P-wave vertical velocity and Thomsen anisotropicparameters {s and 8) could be resolved simultaneously by jointing inversion ofseismic reflection traveltimes and polarizations with the linear inversion methodbased on BG theory.The image of the reflector structure, P-wave vertical velocity and theanisotropy parameters in the crust could be obtained from the wide-angle seismicprofile from Tunxi (in Anhui Province), to Wenzhou (in Zhejiang Province). Theresults reveal the difference of the reflector geometrical structure and physicalattributes in the crust between Yangtze block and Cathaysia block, and attempt tounderstand the characteristics of the crustal stress field in the areas.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

济阳地区深部结构与含油气性讨论

Oil and gas exploration of marine strata in China's Pre-Cenozoic residual basins is regarded as a worldwide puzzle because of existent problems and cruxes. Objectively speaking, the subsurface geologic structure is complicated, and the surface conditions of some areas are tough. On the other hand, there are still many problems to be solved in oil and gas exploration technologies of Pre-Cenozoic marine fades, and theoretic cognition about petroleum geology is not profound yet. Therefore, it is principal to explore integrated geophysical research ways of Pre-Cenozoic residual basins. Seismic prospecting and geophysical integrated interpretation technologies aimed at middle Paleozoic marine facies with deeper burial and complicated geologic conditions have not formed due to bad quality of deep strata data. Pre-Cenozoic strata, and especially extension, thickness and internal structure of Paleozoic strata can not be recognized from seismic profiles, so it is hard to systematically cognize structural features and oil-gas resources prospect of Pre-Cenozoic basins. To further investigation of fabric and structural features, basin prototype, formation and evolution pattern of Pre-Cenozoic basins, and also their control over formation, migration and aggregation of oil and gas, will play a guiding and promotive role in developing new surveying areas, selecting advantageous zones and predicting oil-gas resources.This paper follows the modem macrocontrol theory of "Region constrains local, deep strata controls shallow ones", and uses the integrated geophysical method of "One guide, two hinges, three combinations and multi feedbacks'*. Based on several years of geological and geophysical results of the Shengli Oilfield, and 14 newly-joint regional seismic profiles, deep structure and oil-gas bearing capacity of the Jiyang area are discussed and new cognitions are drawn as below.Seismic identification marks Tr, Tg, Tgl and Tg2 are established for importantPre-Cenozoic geological interfaces, and promoted to the whole Jiyang area.Through area-wide tracking and clogging of important seismic reflection marker,the isochronic framework of pre-Tertiary basin is set up in the Jiyang area for the firsttime, which is vital for basin research.Integrated with geological and geophysical research results, the Jiyang area isdivided into four first-order tectonic sequences- basement, lower tectonic layer,upper tectonic layer, and top tectonic layer. The basement and lower tectonic sequence which are related to Pre-Cenozoic are studied with emphasis.Through the research of regional seismic profiles, the point of view is given thatthe Kongdian Formation of Jiyang is structural transition period. The top-bottomunconformable interface of the Kongdian Formation is found out for the first time,and the basin model is determined primarily, which lay a basis for prototype basinresearch of the Jiyang Kongdian Formation.The distribution status of Middle-Paleozoic is delineated in the Jiyang area.The maximum thickness of Paleozoic lies in the top of the south declivity of half-graben. The thickness gets thinner towards the center of Mesozoic and Cenozoic half-graben basin, and even disappears. Structural action in the west-north affects the distribution of Paleozoic residual strata.6. The features of second-order tectonic sequence of the Jiyang depression isstudied and its evolution history of is rebuilt.Combined with the 5-stage evolution history of the China continent and structure evolution features of the Jiyang area, the structure sedimentary process since Paleozoic is divided into 5 periods - basement forming , Indosinian orogenic, Yanshan negative reversal, Himalayan extension and Neogene subsidence period.Combined with the research results of gravity, magnetic surveying and regionalprofiles, this paper brings forward the idea for the first time that the western boundaryof the Jiyang depression is the Ningjin-Yangpan fracture zone, and forms aside-column assemblage with the Wudi fracture zone.The opinion that under Middle-Cenozoic basins in the middle Jiyang area theremight superimpose an old residual basin is given for the first time. And if it is provedto be true, a new exploration space will be pioneered for Jiyang and even north China.There exists many types of tectonic-stratigraphic traps formed under piezotropy,extension and compound action in Pre-Cenozoic Jiyang. Therein all kinds of burialhills are the most important oil-gas trap type of Pre-Cenozoic, which should besurveyed layeredly according to the layout of oil sources.As such a new challenging project and field, the paper systematically analyses different geophysical responses of the Jiyang area, frames the deep structure of the area, and preliminarily recognizes the Pre-Cenozoic residual basins. It breaks through to a certain extent in both theory and practice, and is expected to provide new geophysical and geotectonic clues for deep exploration in Shengli.

Restructuring and redispersion of silver on SiO2 under oxidizing/reducing atmospheres and its activity toward CO oxidation

The effects of oxygen-hydrogen pretreatments of nanosilver catalysts in cycle mode on the structure and particle size of silver particles, and subsequently the activity of the catalyst toward CO oxidation (or CO selective oxidation in the presence of H-2) are reported in this paper. Ag/SiO2 catalyst with silver particle sizes of ca. 6 similar to 8 nm shows relatively high activity in the present reaction system. The adopting of a cycle of oxidation/reduction pretreatment has a marked influence on the activity of the catalyst. Oxygen pretreatment at 500 degrees C results in the formation of subsurface oxygen and activates the catalyst. As evidenced by in-situ XRD and TEM, the following H-2 treatment at low temperatures (100 similar to 300 degrees C) causes surface faceting and redispersing of the silver particles without destroying the subsurface oxygen species. The subsequent in-situ FTIR and catalytic reaction results show that CO oxidation occurs at -75 degrees C and complete CO conversion can be obtained at 40 degrees C over such a nanosilver catalyst pretreated with oxygen at 500 degrees C followed by H-2 at 100 degrees C. However, prolonged hydrogen treatment at high temperatures (> 300 degrees C) after oxygen pretreatment at 500 degrees C induces the aggregation of silver particles and also depletes so much subsurface oxygen species that the pathway of CO oxidation by the subsurface oxygen species is inhibited. Meanwhile, the ability of the catalyst to adsorb reactants is greatly depressed, resulting in a 20 similar to 30% decrease in the activity toward CO oxidation. However, the activity of the catalyst pretreated with oxygen at 500 degrees C followed by hydrogen treatment at high temperatures (> 300 degrees C) is still higher than that directly pretreated with H,. This kind of catalytic behavior of silver catalyst is associated with physical changes in the silver crystallites because of surface restructuring and crystallite redispersion during the course of oxygen-hydrogen pretreatment steps.

PARAMETERS AFFECTING THE IMAGING AND DETECTION OF MICROBUBBLES IN BLOOD

Stabilized micron-sized bubbles, known as contrast agents, are often injected into the body to enhance ultrasound imaging of blood flow. The ability to detect such bubbles in blood depends on the relative magnitude of the acoustic power backscattered from the microbubbles (‘signal’) to the power backscattered from the red blood cells (‘noise’). Erythrocytes are acoustically small (Rayleigh regime), weak scatterers, and therefore the backscatter coefficient (BSC) of blood increases as the fourth power of frequency throughout the diagnostic frequency range. Microbubbles, on the other hand, are either resonant or super-resonant in the range 5-30 MHz. Above resonance, their total scattering cross-section remains constant with increasing frequency. In the present thesis, a theoretical model of the BSC of a suspension of red blood cells is presented and compared to the BSC of Optison® contrast agent microbubbles. It is predicted that, as the frequency increases, the BSC of red blood cell suspensions eventually exceeds the BSC of the strong scattering microbubbles, leading to a dramatic reduction in signal-to-noise ratio (SNR). This decrease in SNR with increasing frequency was also confirmed experimentally by use of an active cavitation detector for different concentrations of Optison® microbubbles in erythrocyte suspensions of different hematocrits. The magnitude of the observed decrease in SNR correlated well with theoretical predictions in most cases, except for very dense suspensions of red blood cells, where it is hypothesized that the close proximity of erythrocytes inhibits the acoustic response of the microbubbles.

ACOUSTO-OPTIC IMAGING IN DIFFUSE MEDIA USING PULSED ULTRASOUND AND THE PHOTOREFRACTIVE EFFECT

Acousto-optic imaging (AOI) in optically diffuse media is a hybrid imaging modality in which a focused ultrasound beam is used to locally phase modulate light inside of turbid media. The modulated optical field carries with it information about the optical properties in the region where the light and sound interact. The motivation for the development of AOI systems is to measure optical properties at large depths within biological tissue with high spatial resolution. A photorefractive crystal (PRC) based interferometry system is developed for the detection of phase modulated light in AOI applications. Two-wave mixing in the PRC creates a reference beam that is wavefront matched to the modulated optical field collected from the specimen. The phase modulation is converted to an intensity modulation at the optical detector when these two fields interfere. The interferometer has a high optical etendue, making it well suited for AOI where the scattered light levels are typically low. A theoretical model for the detection of acoustically induced phase modulation in turbid media using PRC based interferometry is detailed. An AOI system, using a single element focused ultrasound transducer to pump the AO interaction and the PRC based detection system, is fabricated and tested on tissue mimicking phantoms. It is found that the system has sufficient sensitivity to detect broadband AO signals generated using pulsed ultrasound, allowing for AOI at low time averaged ultrasound output levels. The spatial resolution of the AO imaging system is studied as a function of the ultrasound pulse parameters. A theoretical model of light propagation in turbid media is used to explore the dependence of the AO response on the experimental geometry, light collection aperture, and target optical properties. Finally, a multimodal imaging system combining pulsed AOI and conventional B- mode ultrasound imaging is developed. B-mode ultrasound and AO images of targets embedded in both highly diffuse phantoms and biological tissue ex vivo are obtained, and millimeter resolution is demonstrated in three dimensions. The AO images are intrinsically co-registered with the B-mode ultrasound images. The results suggest that AOI can be used to supplement conventional B-mode ultrasound imaging with optical information.

HEATING IN VASCULAR TISSUE AND FLOW-THROUGH TISSUE PHANTOMS INDUCED BY FOCUSED ULTRASOUND

High intensity focused ultrasound (HIFU) can be used to control bleeding, both from individual blood vessels as well as from gross damage to the capillary bed. This process, called acoustic hemostasis, is being studied in the hope that such a method would ultimately provide a lifesaving treatment during the so-called "golden hour", a brief grace period after a severe trauma in which prompt therapy can save the life of an injured person. Thermal effects play a major role in occlusion of small vessels and also appear to contribute to the sealing of punctures in major blood vessels. However, aggressive ultrasound-induced tissue heating can also impact healthy tissue and can lead to deleterious mechanical bioeffects. Moreover, the presence of vascularity can limit one’s ability to elevate the temperature of blood vessel walls owing to convective heat transport. In an effort to better understand the heating process in tissues with vascular structure we have developed a numerical simulation that couples models for ultrasound propagation, acoustic streaming, ultrasound heating and blood cooling in Newtonian viscous media. The 3-D simulation allows for the study of complicated biological structures and insonation geometries. We have also undertaken a series of in vitro experiments, in non-uniform flow-through tissue phantoms, designed to provide a ground truth verification of the model predictions. The calculated and measured results were compared over a range of values for insonation pressure, insonation time, and flow rate; we show good agreement between predictions and measurements. We then conducted a series of simulations that address two limiting problems of interest: hemostasis in small and large vessels. We employed realistic human tissue properties and considered more complex geometries. Results show that the heating pattern in and around a blood vessel is different for different vessel sizes, flow rates and for varying beam orientations relative to the flow axis. Complete occlusion and wall- puncture sealing are both possible depending on the exposure conditions. These results concur with prior clinical observations and may prove useful for planning of a more effective procedure in HIFU treatments.