砂泥岩储层测井产能预测方法研究

Productivity prediction is very important in the exploration and development of oilfields. Using well log data to predict productivity is a front-line technology， which is key issue in petroleum exploration phase. The essential factors of productivity prediction is building practical models and correcting various causes to improve precision of prediction parameters. Any errors of parameters selections can affect the calculation of productivity prediction; therefore， how to improve research means and calculation accuracy is an important task of productivity prediction. Theory and case examples are deeply and comprehensively studied in the paper. Based on the theory of mud-filtrate invasion and experimental results， the damage of drilling， cementing， perforating，acidizing and fracturing were investigated. The damage depth was quantitatively evaluated by log data， based on this， the processing results of reservoir sensitivity were used to analysis quantitatively the damage of reservoir. The productivity prediction and reservoir damage were initiatively incorporated according to well logging， and the precision of productivity prediction was effectively improved. The method of NMR was explored to calculate the fluid viscosity on the basis of reservoir physical method， and the differences between the two methods were compared in the paper. From the theory fluid flow in porous media， various of theoretical models of production prediction were explored and several practical models were consided， such as productivity index method， improved productivity index method， improved Bearder method， SVM and so on. The characteristic and the application scope of these methods were studied. The inflow productivity and outflow productivity were incorporated and nodal analysis method was used to forecast wellhead yield， thus achieved scientifically production. On the applied background of conventional logging suite， the applying of special items or new logging method which is practical in the research area were studied， the logging suite was further optimized， and the precision of forecast was improved. On the basis of the modeling and the calculation of parameters， these methods were verified and analyzed， and the reconstruct principle was also built for block reservoir. The research block was processed by these methods and compared with testing data. Based on above the research， a technological system which is practical for shaly sand profiles in Shengli Oilfield was built. The system can reach commercialized degree，and satisfied the need of exploration and development of the oilfield.

小湾电站坝基岩体卸荷工程地质力学研究

Dam is the key main works in the construction of water power. The success or failure of the construction of the dam mainly depends on the stability of the dam foundation. The double curvature arch dam-XiaoWan Dam is the highest one among the dams with the same type in the world, and the water thrust acted on it reaches ton, so the rock bearing capacity of dam foundation becomes more important. Because of the high and steep valley-side slope, the large scale of excavation and the complex body type of excavation, it is prominent that the problem of stress release of the rock mass in dam foundation. More great attentions should be paid for the stability and the degraded of rock properties of rock mass induced by the stress release. In this paper, the phenomena of stress release of rock mass in XiaoWan Dam foundation and its mechanisms were analyzed based on the collection of data, the detailed field engineering investigations, measurement of the rock mass and the 2D numerical calculations. The rock mass under the foundation is weak-weathered to intact, the quality of which is good. After excavation of the foundation, the rock mass near the slope surface occurred extend, stretch and stick-slip along original textures till the new fracture surface formed. Then platy structure of the rock mass takes on. The rock mass in the dam foundation occur resilience due to stress release towards free faces with the characteristics of time effect and localized deformation. In-situ measurements show that the rock mass near the surface are degraded. The stress release induced by excavation is a process of the interaction between engineering structures and geologic body. The stress release of rock mass in dam foundation is related to the changed degree of geometrical conditions. The rock near excavation surface failed nearly under uniaxial stresses. The bending-breaking mechanism of plate girder can interpret the failure model of the rock mass with platy structure in dam foundation slope. In essence, the stress release is the change of stress field including the change of directions and magnitudes of stress induced by excavation, which can induce the comedown of the safety margin. In this paper, the inducing conditions of stress release were calculated by numerical analyses. Moreover, from the point of view that the change of stress field, the coefficient of K, i.e. the variable load coefficient was proposed. Then the law of the change of it is interpreted. The distributional characteristics of fracture zone were expressed by the coefficient. The stress release of hard rock has the characteristic of localization. The measuring technique of sound wave can not reflect the small cracks in this kind of rock mass due to stress release. So, the spectral analysis method was proposed. At the same time, the application foreground in engineering of the Stockwell Time-Frequency- Spectrum method was discussed with a view to the limitation of it.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

波动方程偏移速度分析

On the subject of oil and gas exploration, migration is an efficacious technique for imagining structures underground. Wave-equation migration (WEM) dominates over other migration methods in accuracy, despite of higher computational cost. However, the advantages of WEM will emerge as the progress of computer technology. WEM is sensitive to velocity model more than others. Small velocity perturbations result in grate divergence in the image pad. Currently, Kirrchhoff method is still very popular in the exploration industry for the reason of difficult to provide precise velocity model. It is very urgent to figure out a way to migration velocity modeling. This dissertation is mainly devoted to migration velocity analysis method for WEM: 1. In this dissertation, we cataloged wave equation prestack depth migration. The concept of migration is introduced. Then, the analysis is applied to different kinds of extrapolate operator to demonstrate their accuracy and applicability. We derived the DSR and SSR migration method and apply both to 2D model. 2. The output of prestack WEM is in form of common image gathers (CIGs). Angle domain common image gathers (ADCIGs) gained by wave equation are proved to be free of artifacts. They are also the most potential candidates for migration velocity analysis. We discussed how to get ADCIGs by DSR and SSR, and obtained ADCIGs before and after imaging separately. The quality of post stack image is affected by CIGs, only the focused or flattened CIGs generate the correct image. Based on wave equation migration, image could be enhanced by special measures. In this dissertation we use both prestack depth residual migration and time shift imaging condition to improve the image quality. 3. Inaccurate velocities lead to errors of imaging depth and curvature of coherent events in CIGs. The ultimate goal of migration velocity analysis (MVA) is to focus scattered event to correct depth and flatten curving event by updating velocities. The kinematic figures are implicitly presented by focus depth aberration and kinetic figure by amplitude. The initial model of Wave-equation migration velocity analysis (WEMVA) is the output of RMO velocity analysis. For integrity of MVA, we review RMO method in this dissertation. The dissertation discusses the general ideal of RMO velocity analysis for flat and dipping events and the corresponding velocity update formula. Migration velocity analysis is a very time consuming work. Respect to computational convenience, we discus how RMO works for synthetic source record migration. In some extremely situation, RMO method fails. Especially in the areas of poorly illuminated or steep structure, it is very difficult to obtain enough angle information for RMO. WEMVA based on wave extrapolate theory, which successfully overcome the drawback of ray based methods. WEMVA inverses residual velocities with residual images. Based on migration regression, we studied the linearized scattering operator and linearized residual image. The key to WEMVA is the linearized residual image. Residual image obtained by Prestack residual migration, which based on DSR is very inefficient. In this dissertation, we proposed obtaining residual migration by time shift image condition, so that, WEMVA could be implemented by SSR. It evidently reduce the computational cost for this method.

活动构造区区域地壳稳定性研究——以金沙江水电梯级工程规划区为例

The study of regional crustal stability of active tectonic region basically includes analysis of recent activity of Earth's crust, single factor assessment, study of complexity, and comprehensive assessment of crust stability. In this thesis, some work are made as follows: · Based on abundant data from gravity field, aeromagnetic survey, magnetism, magnetotelluric deep sounding, remote sensing and geotectonic as well as earthquakes observed in recent years around this region and adjacent zones, we can get a through understanding about the structural features and activity of the earth's crust in Chuan-Dian region. The results from explosion earthquake and telluric electromagnetic sounding are consistent with the structural features of the crust manifested by the geophysical field. The data of deep geologic structures are important for us to work out a vivid three-dimensional model of the earth's crustal structure of the Jinsha River region. According to a synthesis, the author of this thesis proposes some indicators for dividing the faulted blocks. It can also be inferred that the movement of the Chuan-Dian faulted block, which is the relatively active part of southwestern China, is controlled by the boundary faults, and the intensive activity and deformation are concentrated along the boundaries of the block. · Mainly discussing respectively the mechanism and laws of active faults, earthquakes, and geological hazards activity, and their influences on the stability and security of engineering, also trying to probe into the way to assess the risk of single factor in this section. Especially with the method of fractal geometry, the thesis has discussed how to study the complexity of each factor. These geologic hazards which are distributed at the uppermost part of the crust in this region form a typical mountainous set of the active tectonic areas. The results of survey show that some slopes are liable- to -sliding with a weak layer of low shear strength. Occurrences of landslides are to a great extent related to local geological structures, in particular active faults. This is why numerous landslides have occurred simultaneously around the epicenter of a strong earthquake or the center of a strong rainfall, which are related to active faults. · The analysis of the crustal stability is based on a regional grid division, and a fuzzy comprehensive analysis method is used to determine the grade of the quality in each grid. The evaluation factors and their weights are taken from the results of the hierarchical analysis. The evaluation indexes consist of qualitative and quantitative ones. The qualitative ones can be quantified through the experts weighing system, while the quantitative ones can be obtained from statistical analysis. For quality grades, four levels are used: stable, essentially stable, sub-stable, and unstable. The results of the evaluation on Jinshajiang region demonstrate that the crustal stability become distinctly worse in the areas controlled by active deep faults. Therefore, detailed investigations on the active faulting and geologic hazards, include earthquake activity are especially necessary for those areas adjacent to the deep fault belts. On the bases of the data available and the survey results, we have made a preliminary assessment for the construction conditions and adaptability of every planned site in the middle or lower reaches of Jinsha River. Finally, the thesis prospected the vista of the study of crustal stability.

扩散波勘探方法基础及其在石油油气藏开发中的应用研究

As we know, the essence of exploration is objective body determined by getting the information. Such as seismic、electrical and electromagnetic prospecting, they are the common methods of the exploration. Therefore, They have a complete set of theory now. In fact, the effective information can also be got by the diffusion way, it is called diffusion prospecting. The diffusion way prospecting is necessary and important. The way of diffusion prospecting is studied in the paper and main works include below: (1) On the basis of studying basic law of the diffusion, the paper gives the idea of diffusion wave and the formulas of computing diffusion wave function. (2) The paper studies the way of the diffusion prospecting and the methods of data processing. At the same time, it also expounds the characteristics and the applied foreground of the diffusion prospecting. (3) The paper gives the tomography idea and the basic method of diffusion CT. Meanwhile, it also expounds the foreground that the diffusion CT is applied in oil development prospecting. (4) As the inversion of the diffusion equation is a part of the diffusion prospecting way, the methods of diffusion equation inversion are studied and the two formulas are deduced --Laplace transform and polynomial fitting inversion formulas. As the other important result of diffusion equation inversion, the inversion can offer a new analysis method for well Testing in oil development. In order to show a set of methods in the paper feasible, forward、inversion and CT numerical simulation are done in the paper.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

均质弹性条件下多洞平面问题解析算法及其应用研究

The scholars in the world have been trying to find an effective analytic algorithm of multiple hole problems usually meet in engineering designs. Though some studies on circular or elliptic holes had been achieved under specific conditions, no efforts were made to any multiple hole problems that is most significant for engineering designs. The author has made further studies on any multiple hole problems, using complex variable function method and Schwarz alternating method. After solving a series of technological difficulties, the author obtains an effective analytic algorithm, and acquires stress field and displacement field with high accuracy, which can be conducted for arbitrary many iterations according to practical accuracy requirements. In addition, th solution of stress and displacement fields, even for multiple holes of complex shapes and smaller distances. Further, the author made preliminary studies on viscoelastic displacement solution for any double holes. In terms of the obtained displacement solution of any multiple holes, this paper studies displacement back-analysis for the excavations of two tunnels, and find that the back-analysis method is accurate. Additionally, the author presents the mathematical prove of inversion uniqueness for ground stresses, elastic modulus and Poisson ratio. The author believes that the accurate analytic algorithm provided in this paper will presents an effective way to stress and displacement analysis for any multiple hole problems, optimal arrangement of multiple holes, hole shape optimization of multiple holes, etc..

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

河流成因储层沉积微相与剩余油分布及形成机理研究

Based on the study of fluvial sandstone reservoir in upper of Guantao group in Gudao and Gudong oilfields, this paper first introduces A.D.Miall's(1996a) architectural-element analysis method that was summarized from ground outcrop scale into the reservoir formation research of the study area, more subtly divides sedimentary microfacies and establishes sedimentary model of research area.on this base, this paper summarizes the laws of residual oil distribution of fluvial formation and the control effect of sedimentary microfacies to residual oil distribution, and reveals residual oil formation mechanism. These results have been applied to residual oil production, and the economic effect is good. This paper will be useful for residual oil research and production and enhancement of oil recovery in similar reservoir. The major conclusions of this paper are as follows. 1. Using the architectural-element analysis method to the core data, a interfacial division scheme of the first to the dixth scale is established for the studied fluvial formation. 2.Seven architectural-elements are divided in upper of Guantao group of study area. The sandstone group 5～1+2 of Neogene upper Gutao group belongs to high sinuous fine grain meandering river, and the sandstone group 6 is sandy braided river. 3. Inter layer, the residual oil saturation of "non-main layer" is higher than "main layer", but the residual recoverable reserve of former is larger. Therefore, "main layer" is the main body of residual oil distribution. The upper and middle part of inner layer has lower permeability and strong seeping resistance. Addition to gravity effect in process of driving, its driving efficiency is low; residual oil saturation is high. Because of controlling of inside non-permeable interlayer or sedimentary construction, the residual oil saturation of non-driving or lower driving efficiency position also is high. On plane, the position of high residual oil saturation mostly is at element LV, CS, CH (FF), FF etc, Which has lower porosity and permeability, as well as lens sand-body and sand-body edge that is not controlled by well-net, non-perfect area of injection and production, lower press difference resort area of inter-well diffiuent-line and shelter from fault, local high position of small structure. 4.Microscopic residual oil mainly includes the non-moved oil in the structure of fine pore network, oil in fine pore and path, oil segment in pore and path vertical to flow direction, oil spot or oil film in big pore, residual oil in non-connective pore. 5.The most essential and internal controlling factor of fluvial formation residual oil distribution is sedimentary microfacies. Status of injection and production is the exterior controlling factor of residual oil distribution. 6. The controlling effect of formation sedimentary microfacies to residual oil distribution indicates inter-layer vertical sedimentary facies change in scale of injection and production layer-series, planar sedimentary face change and inner-layer vertical sedimentary rhythm and interbed in single layer to residual oil distribution. 7. It is difficult to clear up the inter-layer difference in scale of injection and production layer-series. The using status of minor layer is not good and its residual oil saturation is high relatively. It is obvious that inter-layer vertical sedimentary facies changes control inter-layer residual oil distribution at the same or similar conditions of injection and production. For fluvial formation, this vertical sedimentary facies change mainly is positive gyration. Namely, from down to top, channel sediment (element CHL, LA) changes into over-bank sediment (element LV, CR, CS). 8. In water-injection developing process of transverse connecting fluvial sandstone oil formation, injection water always comes into channel nearby, and breaks through along channel and orientation of high pressure gradient, does not expand into side of channel until pressure gradient of channel orientation changes into low. It brings about that water-driving status of over-bank sedimentary element formation (LV, CR, CS) is not good, residual oil saturation is high. In non-connective abandoned channel element (CH) formation with channel, because this reverse is difficult to control by injection and production well-series, its using status is not good, even terribly not good, residual oil is enrichment. 9. The rhythm and sedimentary structure, sedimentary facies change in single sand body brings about vertical changes of formation character, growth character of inner layer interbed. These are important factor of controlling and affecting vertical water spread volume and inner layer residual oil forming and distribution in single sand body. Positive rhythm, is the principal part of fluvial sandstone inner layer sedimentary rhythm. Namely, from down to upside, rock grain granularity changes from coarse to fine, seeping ability changes from strong to feebleness. It brings about that water-driving status of inner layer upside is not good, residual oil saturation is high. Inner layer interbed has different degree affecting and controlling effect to seeping of oil and water. Its affecting degree lies on interbed thickness, extending scale, position, and jeted segment of production or injection well. The effect of interbed at upside of oil formation to oil and water seeping is less; the effect of interbed at middle of oil formation to oil and water seeping is more. 10. Indoor experiment and research indicate that wettability, permeability step, vertical permeability, position of Kmax and ratio of oil viscousity and water viscousity all have great effect on the water-driving recovery ratio. 11. Microscopic residual oil distribution is affected and controlled by formation pore network structure, pressure field distribution, and oil characteristic. 12.The residual oil forming mechanism: the over-bank sedimentary element and upper part of a positive rhythm sandstone have fine pore and throat network, permeability is low, displacement pressure of pore and throat is high. The water-driving power usually falls short of displacement pressure that brings about injection water does not spread into these pore and throat network, thereby immovable oil area, namely residual oil, is formed. At underside of channel sedimentary element and positive rhythm sandstone, porosity and permeability is relatively high, connecting degree of pore and throat is high, displacement pressure of pore and throat is low. Thereby injection water is easy to enter into pore and throat, driving oil in them. Because the pore space is irregular, the surface of pore wall is coarse and non-flat. That the oil locate on concave hole of pore wall and the dead angle of pore, and the oil attaches on surface of pore wall by surface tension, are difficult to be peeled off, becoming water-driving residual oil (remaining oil). On the other hand, Because flowing section lessens, flowing resistance increase, action of capillary fore, or seeping speed decreases at process of transfer at pass narrow throat path in the course carried by driving water. The "oil drop", "oil bead", or "oil segment" peeled off by driving water is difficult to carry and to drive out by water at less pressure difference. Thereby they are enclosed in pore to form discontinuous residual oil. 13.This results described above have been applied in nine develop blocks of Gudao and Gudong oilfield. Its applying effect is marked through local injection production adjustment, deploying replacement well, repair hole, replacement envelop, block off water and profile control etc. Relative method and technology can be applied to other oil production area of Shengli oilfield, and obtain better economic and societal effect.

叠前深度偏移速度模型的建立方法研究

Seismic exploration is the main tools of exploration for petroleum. as the society needs more petroleum and the level of exploration is going up, the exploration in the area of complex geology construction is the main task in oil industry, so the seismic prestack depth migration appeared, it has good ability for complex construction imaging. Its result depends on the velocity model strongly. So for seismic prestack depth migration has become the main research area. In this thesis the difference in seismic prestack depth migration between our country and the abroad has been analyzed in system. the tomographical method with no layer velocity model, the residual curve velocity analysical method based on velocity model and the deleting method in pre-processing have been developed. In the thesis, the tomographysical method in velocity analysis is been analyzed at first. It characterized with perfection in theory and diffculity in application. This method use the picked first arrivial, compare the difference between the picked first arrival and the calculated arrival in theory velocity model, and then anti-projected the difference along the ray path to get the new velocity model. This method only has the hypothesis of high frequency, no other hypothesis. So it is very effective and has high efficiency. But this method has default still. The picking of first arrival is difficult in the prestack data. The reasons are the ratio of signal to noise is very low and many other event cross each other in prestack data. These phenomenon appear strongly in the complex geology construction area. Based on these a new tomophysical methos in velocity analysis with no layer velocity model is been developed. The aim is to solve the picking problem. It do not need picking the event time contiunely. You can picking in random depending on the reliability. This methos not only need the pick time as the routine tomographysical mehtod, but also the slope of event. In this methos we use the high slope analysis method to improve the precision of picking. In addition we also make research on the residual curve velocity analysis and find that its application is not good and the efficiency is low. The reasons is that the hypothesis is rigid and it is a local optimizing method, it can solve seismic velocity problem in the area with laterical strong velocity variation. A new method is developed to improve the precision of velocity model building . So far the pattern of seismic prestack depth migration is the same as it aborad. Before the work of velocity building the original seismic data must been corrected on a datum plane, and then to make the prestack depth migration work. As we know the successful example is in Mexico bay. It characterized with the simple surface layer construction, the pre-precessing is very simple and its precision is very high. But in our country the main seismic work is in land, the surface layer is very complex, in some area the error of pre-precessing is big, it affect the velocity building. So based on this a new method is developed to delete the per-precessing error and improve the precision of velocity model building. Our main work is, (1) developing a effective tomographical velocity building method with no layer velocity model. (2) a new high resolution slope analysis method is developed. (3) developing a global optimized residual curve velocity buliding method based on velocity model. (4) a effective method of deleting the pre-precessing error is developing. All the method as listed above has been ceritified by the theorical calculation and the actual seismic data.

陕甘宁盆地储层物性预测及油气综合判识方法研究

The thesis show you seven interpretation models of erosion of MAs1+2 in the west of prospect and eroded gully of middle and east, form the interpretation technique how to built up eroded gully of the Ordovician top, according of the practical demand of oil-gas exploration in the ShanGanNing basin, using seismic information, combining well logging and well drilling data, Carefully analyzing geologic deposition background and well logging data, through a great quantity forward and inversion for geologic model and combination geologic model with seismic section. Related to research of reservoir absorption in the ShanGanNing basin, it firstly introduces PRONY transformation multidimensional filter. It can simultaneously express relationship of frequency and absorption decay coefficient, better than FUSAIPU analysis method; PRONY filter have obtain the better effect in the gas field of ZhenChuanBao in the ShanBei area after adopting PRONY filtering method to predict reservoir absorption, by analyzing fixed well and prediction of non-well drilling. In the ShanGanNing basin, general seismic inversion method can produce evident different results or misunderstanding because wave impedance and lithology, physical property, gas property are not sole, especially while have little impedance contrast and even have contract direction; the author carefully analyzes multi-parameter inversion technique, add natural gamma ray and natural potential and other parameter combined making model inversion method according of theory of seismic inversion and applying reservoir velocity and wave impedance information at last, we get the more directly reservoir physical property parameter, judging reservoir physical property is more exact. In accordance with geologic, seismic feature of Shan basin, the thesis conclude Ordovician system top erosion interpretation technology with ChangQing character, and reservoir thickness prediction technique combining inversion technique with wave character analysis, Reservoir physical property that is mainly absorption factor analysis and multi-parameter inversion and oil-gas prediction technology. These technologies obtain the better result in the oil-gas field exploration and have formed comprehensive research method and technology series with ShanGanNing character.