VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.

基于波动方程的层间多次波预测方法研究

In exploration geophysics，velocity analysis and migration methods except reverse time migration are based on ray theory or one-way wave-equation. So multiples are regarded as noise and required to be attenuated. It is very important to attenuate multiples for structure imaging, amplitude preserving migration. So it is an interesting research in theory and application about how to predict and attenuate internal multiples effectively. There are two methods based on wave-equation to predict internal multiples for pre-stack data. One is common focus point method. Another is inverse scattering series method. After comparison of the two methods, we found that there are four problems in common focus point method: 1. dependence of velocity model; 2. only internal multiples related to a layer can be predicted every time; 3. computing procedure is complex; 4. it is difficult to apply it in complex media. In order to overcome these problems, we adopt inverse scattering series method. However, inverse scattering series method also has some problems: 1. computing cost is high; 2. it is difficult to predict internal multiples in the far offset; 3. it is not able to predict internal multiples in complex media. Among those problems, high computing cost is the biggest barrier in field seismic processing. So I present 1D and 1.5D improved algorithms for reducing computing time. In addition, I proposed a new algorithm to solve the problem which exists in subtraction, especially for surface related to multiples. The creative results of my research are following: 1. derived an improved inverse scattering series prediction algorithm for 1D. The algorithm has very high computing efficiency. It is faster than old algorithm about twelve times in theory and faster about eighty times for lower spatial complexity in practice; 2. derived an improved inverse scattering series prediction algorithm for 1.5D. The new algorithm changes the computing domain from pseudo-depth wavenumber domain to TX domain for predicting multiples. The improved algorithm demonstrated that the approach has some merits such as higher computing efficiency, feasibility to many kinds of geometries, lower predictive noise and independence to wavelet； 3. proposed a new subtraction algorithm. The new subtraction algorithm is not used to overcome nonorthogonality, but utilize the nonorthogonality's distribution in TX domain to estimate the true wavelet with filtering method. The method has excellent effectiveness in model testing. Improved 1D and 1.5D inverse scattering series algorithms can predict internal multiples. After filtering and subtracting among seismic traces in a window time, internal multiples can be attenuated in some degree. The proposed 1D and 1.5D algorithms have demonstrated that they are effective to the numerical and field data. In addition, the new subtraction algorithm is effective to the complex theoretic models.

基于精细结构描述及数值试验的断续节理岩体力学参数智能选取

Evaluating the mechanical properties of rock masses is the base of rock engineering design and construction. It has great influence on the safety and cost of rock project. The recognition is inevitable consequence of new engineering activities in rock, including high-rise building, super bridge, complex underground installations, hydraulic project and etc. During the constructions, lots of engineering accidents happened, which bring great damage to people. According to the investigation, many failures are due to choosing improper mechanical properties. ‘Can’t give the proper properties’ becomes one of big problems for theoretic analysis and numerical simulation. Selecting the properties reasonably and effectively is very significant for the planning, design and construction of rock engineering works. A multiple method based on site investigation, theoretic analysis, model test, numerical test and back analysis by artificial neural network is conducted to determine and optimize the mechanical properties for engineering design. The following outcomes are obtained: (1) Mapping of the rock mass structure Detailed geological investigation is the soul of the fine structure description. Based on statistical window，geological sketch and digital photography，a new method for rock mass fine structure in-situ mapping is developed. It has already been taken into practice and received good comments in Baihetan Hydropower Station. (2) Theoretic analysis of rock mass containing intermittent joints The shear strength mechanisms of joint and rock bridge are analyzed respectively. And the multiple modes of failure on different stress condition are summarized and supplied. Then, through introducing deformation compatibility equation in normal direction, the direct shear strength formulation and compression shear strength formulation for coplanar intermittent joints, as well as compression shear strength formulation for ladderlike intermittent joints are deducted respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. (3) Model test of rock mass containing intermittent joints Model tests are adopted to study the mechanical mechanism of joints to rock masses. The failure modes of rock mass containing intermittent joints are summarized from the model test. Six typical failure modes are found in the test, and brittle failures are the main failure mode. The evolvement processes of shear stress, shear displacement, normal stress and normal displacement are monitored by using rigid servo test machine. And the deformation and failure character during the loading process is analyzed. According to the model test, the failure modes quite depend on the joint distribution, connectivity and stress states. According to the contrastive analysis of complete stress strain curve, different failure developing stages are found in the intact rock, across jointed rock mass and intermittent jointed rock mass. There are four typical stages in the stress strain curve of intact rock, namely shear contraction stage, linear elastic stage, failure stage and residual strength stage. There are three typical stages in the across jointed rock mass, namely linear elastic stage, transition zone and sliding failure stage. Correspondingly, five typical stages are found in the intermittent jointed rock mass, namely linear elastic stage, sliding of joint, steady growth of post-crack, joint coalescence failure, and residual strength. According to strength analysis, the failure envelopes of intact rock and across jointed rock mass are the upper bound and lower bound separately. The strength of intermittent jointed rock mass can be evaluated by reducing the bandwidth of the failure envelope with geo-mechanics analysis. (4) Numerical test of rock mass Two sets of methods, i.e. the distinct element method (DEC) based on in-situ geology mapping and the realistic failure process analysis (RFPA) based on high-definition digital imaging, are developed and introduced. The operation process and analysis results are demonstrated detailedly from the research on parameters of rock mass based on numerical test in the Jinping First Stage Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Then the applicable fields are figured out respectively. (5) Intelligent evaluation based on artificial neural network (ANN) The characters of both ANN and parameter evaluation of rock mass are discussed and summarized. According to the investigations, ANN has a bright application future in the field of parameter evaluation of rock mass. Intelligent evaluation of mechanical parameters in the Jinping First Stage Hydropower Station is taken as an example to demonstrate the analysis process. The problems in five aspects, i. e. sample selection, network design, initial value selection, learning rate and expected error, are discussed detailedly.

中国大中型气田天然气成藏物理化学模拟研究

Natural gas pays more important role in the society as clean fuel. Natural gas exploration has been enhanced in recent years in many countries. It also has prospective future in our country through "85" and "95" national research. Many big size gas fields have been discovered in different formations in different basins such as lower and upper Paleozoic in Erdos basin, Tertiary system in Kuche depression in Tarim basin, Triassic system in east of Sichuan basin. Because gas bearing basins had been experienced multiple tectogenesis. The characteristics of natural gases usually in one gas field are that they have multiple source rocks and are multiple maturities and formed in different ages. There has most difficult to research on the gas-rock correlation and mechanism of gas formation. Develop advanced techniques and methods and apply them to solve above problems is necessary. The research is focused on the critical techniques of geochemistry and physical simulation of gas-rock correlation and gas formation. The lists in the following are conclusions through research and lots of experiments. I 8 advanced techniques have been developed or improved about gas-rock correlation and gas migration, accumulation and formation. A series of geochemistry techniques has been developed about analyzing inclusion enclave. They are analyzing gas and liquid composition and biomarker and on-line individual carbon isotope composition in inclusion enclave. These techniques combing the inclusion homogeneous temperature can be applied to study on gas-rock correlation directly and gas migration, filling and formation ages. Technique of on-line determination individual gas carbon isotope composition in kerogen and bitumen thermal pyrolysis is developed. It is applied to determine the source of natural is kerogen thermal degradation or oil pyrolysis. Method of on-line determination individual gas carbon isotope composition in rock thermal simulation has being improved. Based on the "95"former research, on-line determination individual gas carbon isotope composition in different type of maceral and rocks thermal pyrolys is has been determined. The conclusion is that carbon isotope composition of benzene and toluene in homogenous texture kerogen thermal degradation is almost same at different maturity. By comparison, that in mixture type kerogen thermal pyrolysis jumps from step to step with the changes of maturity. This conclusion is a good proof of gas-rock dynamic correlation. 3. Biomarker of rock can be determined directly through research. It solves the problems such as long period preparing sample, light composition losing and sample contamination etc. It can be applied to research the character of source rock and mechanism of source rock expulsion and the path of hydrocarbon migration etc. 4. The process of hydrocarbon dynamic generation in source rock can be seen at every stage applying locating observation and thermal simulation of ESEM. The mechanism of hydrocarbon generation and expulsion in source rock is discussed according to the experiments. This technique is advanced in the world. 5. A sample injection system whose character is higher vacuum, lower leaks and lower blank has been built up to analyze inert gas. He,Ar,Kr and Xe can be determined continuously on one instrument and one injection. This is advanced in domestic. 7. Quality and quantity analysis of benzene ring compounds and phenolic compounds and determination of organic acid and aqueous gas analysis are applied to research the relationship between compounds in formation water and gas formation. This is another new idea to study the gas-rock correlation and gas formation. 8. Inclusion analysis data can be used to calculate the Paleo-fluid density, Paleo-geothermal gradient and Paleo-geopressure gradient and then to calculate the Paleo-fluid potential. It's also a new method to research the direction of hydrocarbon migration and accumulation. 9. Equipment of natural gas formation simulation is produced during the research to probe how the physical properties of rock affect the gas migration and accumulation and what efficiency of gas migrate and factors of gas formation and the models of different type of migration are. II study is focused on that if the source rocks of lower Paleozoic generated hydrocarbon and what the source rocks of weathered formation gas pool and the mechanism of gas formation are though many advanced techniques application. There are four conclusions. 1.The maturity of Majiagou formation source rocks is higher in south than that in north. There also have parts of the higher maturity in middle and east. Anomalous thermal pays important role in big size field formation in middle of basin. 2. The amount of gas generation in high-over maturity source rocks in lower Paleozoic is lager than that of most absorption of source rocks. Lower Paleozoic source rocks are effective source rocks. Universal bitumen exists in Ordovician source rocks to prove that Ordovician source rocks had generated hydrocarbon. Bitumen has some attribution to the middle gas pool formation. 3. Comprehensive gas-rock correlation says that natural gases of north, west, south of middle gas field of basin mainly come from lower Paleozoic source rocks. The attribution ratio of lower Paleozoic source rocks is 60%-70%. Natural gases of other areas mainly come from upper Paleozoic. The attribution ratio of upper Paleozoic source rocks is 70%. 4. Paleozoic gases migration phase of Erdos basin are also interesting. The relative abundance of gasoline aromatic is quite low especially toluene that of which is divided by that of methyl-cyclohexane is less than 0.2 in upper Paleozoic gas pool. The migration phase of upper Paleozoic gas may be aqueous phase. By comparison, the relative abundance of gasoline aromatic is higher in lower Paleozoic gas. The distribution character of gasoline gas is similar with that in source rock thermal simulation. The migration phase of it may be free phase. IH Comprehensive gas-rock correlation is also processed in Kuche depression Tarim basin. The mechanism of gas formation is probed and the gas formation model has been built up. Four conclusions list below. 1. Gases in Kuche depression come from Triassic-Jurassic coal-measure source rocks. They are high-over maturity. Comparatively, the highest maturity area is Kelasu, next is Dabei area, Yinan area. 2. Kerogen thermal degradation is main reason of the dry gas in Kuche depression. Small part of dry gas comes from oil pyrolysis. VI 3.The K12 natural gas lays out some of hydro-gas character. Oil dissolved in the gas. Hydro-gas is also a factor making the gas drier and carbon isotope composition heavier. 4. The mechanism and genesis of KL2 gas pool list as below. Overpressure has being existed in Triassic-Jurassic source rocks since Keche period. Natural gases were expulsed by episode style from overpressure source rocks. Hetero-face was main migration style of gas, oil and water at that time. The fluids transferred the pressure of source rocks when they migrated and then separated when they got in reservoir. After that, natural gas migrated up and accumulated and formed with the techno-genesis. Tectonic extrusion made the natural gas overpressure continuously. When the pressure was up to the critical pressure, the C6-C7 composition in natural gas changed. The results were that relative abundance of alkane and aromatic decreased while cycloalkane and isoparaffin increased. There was lots of natural gas filling during every tectonic. The main factors of overpressure of natural gas were tectonic extrusion and fluid transferring pressure of source rocks. Well preservation was also important in the KL2 gas pool formation. The reserves of gas can satisfy the need of pipeline where is from west to east. IV A good idea of natural gas migration and accumulation modeling whose apparent character is real core and formation condition is suggested to model the physical process of gas formation. Following is the modeling results. 1. Modeling results prove that the gas accumulation rule under cap layer and gas fraction on migration path. 2. Natural gas migration as free phase is difficult in dense rock. 3. Natural gases accumulated easily in good physical properties reservoirs where are under the plugging layer. Under the condition of that permeability of rock is more than 1 * 10~(-3)μm~(-1), the more better the physical properties and the more bigger pore of rock, the more easier the gas accumulation in there. On the contrary, natural gas canonly migrate further to accumulate in good physical properties of rock. 4. Natural gas migrate up is different from that down. Under the same situation, the amount of gas migration up is lager than that of gas migration down and the distance of migration up is 3 times as that of migration down. 5. After gas leaks from dense confining layer, the ability of its dynamic plug-back decreased apparently. Gas lost from these arils easily. These confining layer can confine again only after geology condition changes. 6. Water-wetted and capillary-blocking rocks can't block water but gases generally. The result is that water can migrate continuously through blocking rocks but the gases stay under the blocking rocks then form in there. The experiments have proved the formation model of deep basin gas.

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

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

中国粮食生产的区域演进研究

Grain is one of the primary material conditions of the human survival and the grain production concerns the stability and development of the society directly. The regional patterns influence greatly on the grain production and the rational production distribution the regional comparative advantages and promotes grain production. This thesis starts with summarizing of the characteristics of changes and the overall trend of regional pattern of grain production of our country since 1949. Then it carries on network analyses to the factors, which influences the evolvement of regional grain production patterns of our country. And finally it gives some proposals to the grain production distribution in the future. The main content includes: Firstly, Reviewing the regional evolvement of grain production in our country, and analyzing the changes of the regional pattern of grain production of our country on the provincial scale and county scale separately, since 1949, especially since the reform and opening up policy. The main grain production areas are acting an important position in ensuring the national grain security, so this thesis analyses the main matter of the main grain production areas, forecasts the grain production situation in the future, and selects the Northeastern main grain production areas as the typical area to carry on the positive research. Secondly, this thesis analyzes the origin causes from two respects of natural and social economy of the regional evolvement pattern of grain production in China. Thirdly, based on the summarizing to the status of the regional pattern of the grain production, this thesis proposes the precept of the grain production distribution in the future in our country. Therefore, the areas of three major cereal crops, rice, wheat and corn, are confirmed on the basis of the comparative advantages. Finally, this thesis puts forward the security system of guaranteeing the grain production progressing steady in China. According to the above analysis, some conclusions have been achieved as follows: (1) The grain gross production gets on extricating itself from awkward position frequently while fluctuating greatly annually since 1949 in China. (2) Since the reform, its traditional regional pattern of grain production, the most of which was concentrated in the south area, has changed rapidly. China's center of gravity of grain production has shifted from the south to the north, and on the belts of latitude, the grain production has represented a trend of focusing to the middle area in China. (3) The main grain production areas play a very important role in ensuring China's food security. With their relative severe situation of the problems of agriculture, rural area and peasant, China has carried out a series of measures, which aim at improving the food-producing conditions of the main grain production areas, and enhancing the grain yields there. Under this condition, a forecast of the producing amount of the main grain production areas under the nation's self-supplying rate of over 95% shows that the increasing provision production in these areas can meet the demand of the country. (4) The natural and social economic factors influence together on the changes of the grain production regional pattern. Along with the state system transition and progress of agricultural science and technology, the regional pattern of grain production is affected heavier by the agricultural policy and technological elements. (5) The grain production will be concentrated to the middle province in the future, which economic development level being medium-sized; According to crop allocation, although the rice superiority production area located in the South, its comparative advantage index is little in some degree. Meanwhile, the wheat and corn superiority production areas are in the North mainly and its scale superiority and production level advantage are all comparatively obviously.