水利水电岩溶发育深度的定量化理论研究

The distribution in outcrop of carbonate rocks in China is about 3.44×106km2 and accounts for over one third of the total area of the country, mainly in the Yangtze valley and the Yungui tableland ,and there is a great plenty of Water and hydropower resources in these carbonate area. A large number of projects will be built in order to develop the Water and hydropower resources. They are facing amount of complicated problems of karsts, specially the problem of the depths of karsts below water, which is a key problem that hasn’t been solved well theoretically. So, systematically research in quantitative analysis of the problem is necessary, which is very important to solve the inconsistency between current theories and facts, and is helped to foresee the depth of karsts below water and decrease the expense of the engineering. For the problem mentioned above, the thesis makes a detailed research on the length of corrosion of karsts water, based on the kinetics of corrosion, and founds the primary theory about the processes of four typical geologic elements, which are single fracture, single conduit, porous limestone and fractured limestone. And it has done a deep research upon the characterties and kinetics of corrosion of the four typical geologic elements by simulation using the programme edited by the author. The thesis also makes a discussion on the general model and process of the form of caves by using the founded theory. According to the characteristic of Water Resource and Hydropower engineering, the thesis creates three representative flow-dissolution models for three types of common geologic conditions under the induction of the theory kinetics of corrosion of karsts water. It has done quantitative research upon the process of dissolution and brought out primary theory about quantitative analysis of the depth of karsts below water as well. It found that the depth of karsts below water relates to the specific geologic conditions and the time of dissolution by simulation based on typical parameters and deep analysis of the result. That is to say there are karsts caves in any depths in flow areas of groundwater in specific geologic conditions and appropriate time. The thesis also discourses three basic problems that frequently encountered in the reconnaissance of Water Resource and Hydropower engineering, which are karsts base level, essential conditions of karsts, the comparability of caves and terrace, and indicates the importance of time in karsts. Finally, as an application, the karsts of the scheming Qianzhong water resource engineering is analyzed by the primary theory founded of quantitative analysis. The result of the application comes to a conclusion that the theory accords with the facts properly.

黄海海域前新生代残留盆地综合地质地球物理研究

China’s annual oil import volume has been increasing in recent years, but the oil price in the international market fluctuates and poses a severe threat to China’s economic development and national security. Therefore, it is of great importance to study the gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea. Yellow Sea has widespread and thick Mesozoic and Paleozoic strata that contain multilayer source rock. Hence, Yellow Sea Mesozoic and Paleozoic strata have good conditions of forming Pre-Cenozoic hydrocarbon reservoirs. Pre-Cenozoic Residual Basins are usually buried deep and then transformed many times in its long evolutional history. These characteristics make it difficult to apply a single method in exploring Pre-Cenozoic Residual Basins. On the other hand, it is highly effective to solve key problems of gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea by using integrated geological and geophysical methods which make full use of the advantages of various exploring techniques. Based on the principle of “the region controls the local; the deep restricts the shallow,” this study focuses on Pre-Cenozoic Residual Basins in Yellow Sea to describe the structure frame of its distribution, with gravity, magnetic, seismic, drill-hole and geological data and previous research findings. In addition, the distribution characteristics of Pre-Cenozoic Residual Basins in Yellow Sea are also analyzed. This paper explores the characteristics of error between gravity forward with constant density and gravity forward with variable density through the study on 2-D and 3-D gravity forward in frequency domain. The result shows that there is a linear relationship between error and depth of 2-D geological model but there is a nonlinear relationship between error and depth of 3-D geological model. The error can be removed according to its linear characteristics or statistical nature of nonlinear characteristics. There is also error between gravity inversion with constant density and gravity inversion with variable density due to variable density and edge-effect. Since there are not noticeable rules between the error and the two causes as variable density and edge-effect, this study adopts gravity inversion with variable density and methods to eliminate the edge-effect in basement inversion to improve inversion accuracy. Based on the study on the rock physical properties and strata distribution of Yellow Sea and adjacent regions, this study finds that there is a big density contrast between Cretaceous-Jurassic strata and their substratum. The magnetic basement of south Yellow Sea is regarded as top of Archeozoic-Proterozoic early strata, and there are double magnetic basements in north Yellow Sea. Gravity and magnetic data are used to inverse the gravity basement and magnetic basement of Yellow Sea, with seismic and drill-hole data as constrains. According to data of gravity and magnetic basement distribution, the depth of Cenozoic strata and previous research findings, this paper calculates the thickness of the Mesozoic and Pre-Mesozoic Residual Basins, draws the distribution outline of Pre-Cenozoic Residual Basins in Yellow Sea, and analyzes its macro-distribution characteristics. Gravity inversion is applied on a typical geological profile in Yellow Sea to analyze the characteristics of its fractures and magnetic basements. The characteristics of Pre-Cenozoic Residual Basins distribution outline in Yellow Sea and the fractures and magnetic basements of its typical profile shown by profile inversion provides new geophysical evidence for these structure views such as “the South Yellow Sea and the North Yellow Sea belong to different structural units” and “Sino-Korea and Yangtze blocks combine along Yellow Sea East Fractured Zone in Yellow Sea”.

AVO分析与地震资料叠前储层预测技术

In this paper we base on the anisotropic theory and Zoeppritz function of the transmission theory and the law of amplitude versus offset simplify seismic reflection coefficient of different media, analyze the characteristic of the gas or oil saturated stratum or the VTI and HTI models. Discuss the P wave reflection relationship and the meanings of the different parameters. We use measured parameters of a reservoir to simulate the characteristic of the reservoir, study the different effects of stratum saturated with gas or oil and analyze the characteristic of the seismic response of different models which change with different incident angles and different azimuths. Using the field data of logs ,analyze the rock property parameters, build the relationship of logs and parameters by Gassmann theory or empirical function. Calculate the density and the shear modulus and bulk modulus, reconstruct the log curves, calculate shear wave logs and correlate the logs affected by mud and other environmental factors. Finally perform the relationship of the seismic data log of saturated stratum and enhance the ability and reliability in reservoir prediction. Our aim is by the prestack seismic processing to get high solution and amplitude preserved seismic data. Because in incident angle gathers or azimuthal gathers, the low signal to noise ratio and low different covers affect the result of the prestack reservoir prediction. We apply prestack noise erase, cell regularization process and relatively amplitude preservation in the high solution seismic process routine to preserve the characteristic of stratum response, and erase the effects of the noise. In this paper we finished prestack invertion in the BYT survey and fractured reservoir depiction in MB survey. By the invertion and multiple attributes crossplot. we can get the stratum profiles and oil indicator profiles which can predict the distribution of the reservoir and oil. In the MB survey, we get orientation and density of fractured reservoir by the azimuthal seismic amplitude and depict the potential oil and gas reservoir. Prestak invertion works better in distinguishing oil and reservoir.

束鹿凹陷泥灰岩储层综合评价与加砂压裂技术研究与应用

Luo Ning ( Mineralogy, Petrology, Deposit Mineralogy) Directed by Fu Liyun With the increase of the level of exploration and development, North China field, as one of the maturing fields in the east, has gradually turned their prospecting targets to frontiers such as deep zones, lithologic hydrocarbon reservoirs, low permeable layers, special lithostromes, etc, which propose new challenges to mating technique of exploration engineering. In it, the special lithostrome of clay carbonate in Shu-Lu cave in Middle Flank exploration area locates in Es_3 generating rock. The area distribution is large, formation thickness is over 100 meters, the oil accumulation condition is excellent, prognostic reserves is over 80,000,000 tons, but how to effectively stimulate the special low permeable and fractured reservoir has become the bottle neck problem of stimulation and stable yields. In this thesis, through comprehensive evaluation and analysis of lithology, lithomechanics, hydrocarbon reservoir characteristics, the characteristics of fluid flow through porous medium and the stimulation measures in the past, we acquire new cognition of clay carbonate reservoirs, in addition, the research and application of first hydraulic fracturing has gained positive effect and formed commensurable comprehensive reservoir evaluation technique and mating engineering technique of hydraulic fracturing. The main cognitions and achievements are as follows: 1.Study of geological information such as lithololy analysis and nuclear magnetic logging, etc, indicates that clay carbonate formation of Shu-Lu cave is anisotropic, low permeable with high shale content, whose accumulation space gives priority to microcracks. 2.The analysis of lithomechanics of clay carbonate indicates that the hardness is moderate, Young’s modulus is between that of sandstone and limestone, clay carbonate presents plastic property and its breakdown pressure is high because of the deep buried depth. 3.The analysis of the drillstem test curves indicates that the flow and build-up pressure curve of clay carbonate of Shu-Lu cave mainly has three types: formation contamination block-up type, low permeable type, formation energy accumulation slowness type; the reservoir characteristics presents double porosity media, radial compounding, uniform flow vertical fracture, isotropy, moniliform reservoir type. The target well Jingu 3 belongs to moniliform reservoir type. 4.Through recognition and re-evaluation of the treatment effect and technologic limitations of acidizing, acid fracturing and gelled acidizing in the past, based on the sufficient survey and study of hydraulic fracturing home and abroad, combined with comprehensive formation study of target well, we launched the study of the optimization of hydraulic fracturing technique, forming the principal clue and commensurable mating technology aimed at clay carbonate formation, whose targets are preventing leak off, preventing sand bridge, preventing embedment, controlling fracture height, forming long fracture. 5. Recognition of stimulation effect evaluation.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

滑坡型泥石流转化机理研究—以四川省达县青宁乡滑坡为例

Slide-debris flow is debris flow which is transformed from landslide consecutively in a short time, it comprises of two phases: First, Landslide starts to slide; Second, landslide changes to debris flow. Slide-debris flow which brings great property and life loss happens frequently at home and abroad. In order to forecast the happening possibility and scope of slide-debris flow, transfromation mechanism of Slide-debris flow must be studied. Research on transformation mechanism of slide-debris flow is intersectant science of landslide kinetics and debris flow starting theory, It is a fringe problem as well as front problem of geological hazard. This paper takes Qingning slide-debris flow in Da County, Sichuan Province for example and has studied the mechanism of its instability and transfromation into debris flow through indoor test (including usual soil test and ring shear test) and digital modeling method.The research gets the following conclusions. Qingning Landslide took place mainly because of confined water head arising from rainfall infiltration. Before Landslide occurring, it rained continuously for 22 days, accumulated precipitation arrived at 521.6mm.Investigation shows that strata of Qingning Landslide contains quaternary loose accumulation, slip soil and highly weathered bedrock, which is a good condition for formation of confined water in the slope. Further more, groundwater seepage in the slope body and corresponding slope safety factor before landslide occurring have been computed through finite element method. The result shows that because of infiltration of rainfall, confined water head in the slope arose sharply, accordingly, the safety factor of the slope declined quickly. The result also shows that force put on the slide body by the rock mass detached from Dazhaiyan mountain was the direct factor for landslide occurring. Qingning slide-debris transformation mode has been summarized, the process the landslide changed into debris flow is divided into three phases in the prospective of macroscopic geological condition: landslide occurring, transformation and debris flow. Landslide occurring phase is from slope’ local creeping slide to Landslide occurring; transformation phase contains slide body sliding on the slide bed after slide occurring and sliding on the slope after shearing opening; debris flow phase is that slide body breaks up completely and flows downward into the ditches. The transformation mechanism of Qingning slide-debris flow has been studied through indoor ring shear test of slip soil. The result shows that transformation mechanism contains two points: first, during slide body sliding on the slide bed and slope after shearing opening, shearing shrinkage, grain crushing and grain layering brought about declining of its volume and produced excess pore water pressure, and because producing velocity of excess pore water pressure is much greater than its dissipating velocity, shear strength of slide body decreased sharply because of accumulated pore water pressure. Second, grains crushing and grains layering during slide body sliding brought about thick liquefied layer at the bottom of the slidebody, liquefied layer contained high water content and its shear strength was very low, its thickness increased as the sliding displacement increasing. Liquefied layer makes slide body sliding fast and easily break down to debris flow. Excess pore water pressure and liquefied layer made shear strength of slidebody became very low, furthermore, water in the pit of slope joining in the slidebody was also a facter that made slidebody accelerate the transformation. Influence of slide body thickness and fine grains content to transformation of slide-debris flow has been studied through ring shear test. The result reaches two conclusions. First, thickness of slide body affects transformation of slide-debris flow by two ways, porewater pressure and effect of “soft base” increases as thickness of slide body increasing.so the thicker slide body is ,the easier transformation is. Second, actual dissipating velocity of porewater pressure should be considered when studying the influence of fine grains content to tranformation of slide-debris flow. There should be a critical content of fine grains which makes the difference of producing and dissipating velocity of water pore pressre greatest, this value is the best for slide-debris transformation. The whole process of slide-debris flow transformation is reproduced through discrete element method. Transformation mechanism of slide-debris flow is studied through monitoring various parameters including pore water pressure, grain crushing and grain layering in the slide body during the transformation. The result confirms and supplements the transformation mechanism of slide-debris flow got from ring shear test well.

裂缝型储层检测方法——P波方位各向异性技术

Nowadays, the exploration of fractured reservoir plays a vital role in the further development of petroleum industry through out the world. Fractured hydrocarbon reservoirs are widely distributed in China. Usually, S-wave technique prevails, but it also has its disadvantage, prohibitive expense in S-wave data acquisition and processing. So directly utilizing P-wave data to detect fractures, comes to our mind. We briefly introduce theoretical model (HTI) for fractured reservoir. Then study Ruger’s reflectivity method to recognize reflection P-wave reflection coefficient of the top and bottom interface of HTI layer respectively, and its azimuth anisotropy character. Base on that study, we gives a review and comparison of two seismic exploration technologies for fractures available in the industry-- P-wave AVO and AVA. They has shown great potential for application to the oil and gas prediction of fractured reservoir and the reservoir fine description．Every technique has its disadvantage, AVO limited to small reflection angle; and AVA just offering relatively results. So that, We can draw a conclusion that a better way to any particular field is using synthesis of multiple data sources including core、outcrop、well-test、image logs、3D VSPs, generally to improve the accuracy.

块体理论赤平投影解析法与关键块体滑落概率研究—以锦屏二级水电站皮带机隧道为例

Block theory is an effective method on stability analysis of fractured rigid rock mass. There are a lot of discontinuous planes developed in rock mass of Jinping II hydropower station conveyor tunnel, so the stability of conveyor tunnel is related with whether there are unstable blocks on excavation planes. This paper deals with the stability of conveyor tunnel with stereo-analytical method for block theory on the basis of detailed investigation of rock mass data, and makes judgements on the movable blocks sliding types which are induced by all rock discontinuous planes and every excavation plane of conveyor tunnel. A conclusion is obtained that the sliding type of blocks is mainly single sliding, and a relatively few sliding types of double-sided sliding and vertical block falling; Also, the obvious statistical distribution result on movable blocks in conveyor tunnel indicates that there are a bit more instability blocks in left wall, left and right arches than right wall. In this paper, the stochastic probability model is drawn into block theory to study the sliding probability of key block on the basis of detailed investigation of its rock mass data and the development of the discontinuous planes in rock mass of Jinping II hydropower station conveyor tunnel. And some following conclusions are obtained. The relationship between trace length and the probability of instability of key block is inverse ratio. The probability of 1-3m primary joints are relatively higher. Key block containing joints J2 is relatively stable and the reinforcement of the arch would be crucial in the conveyor tunnel. They are all useful to offer effective reinforcement design and have important engineering values.

纵波方位属性裂缝检测技术综合研究

Fractured oil and gas reservoir is an important type of oil and gas reservoir, which is taking a growing part of current oil and gas production in the whole world. Thus these technologies targeted at exploration of fractured oil and gas reservoirs are drawing vast attentions. It is difficult to accurately predict the fracture development orientation and intensity in oil and gas exploration. Focused on this problem, this paper systematically conducted series study of seismic data processing and P-wave attributes fracture detection based on the structure of ZX buried mountain, and obtained good results. This paper firstly stimulated the propagation of P-wave in weak anisotropic media caused by vertical aligned cracks, and analyzed the rule of P-wave attributes’ variation associated with observed azimuth, such as travel-time, amplitude and AVO gradient and so on, and quantitatively described the sensitive degree of these attributes to anisotropy of fracture medium. In order to further study the sensitive degree of these attributes to anisotropy of fractures, meanwhile, this paper stimulated P-wave propagation through different types and different intensity anisotropic medium respectively and summarized the rule of these attributes’ variation associated with observed azimuth in different anisotropic medium. The results of these studies provided reliable references for predicting orientation, extensity and size of actual complicated cracked medium by P-wave azimuth attributes responses. In the paper, amounts of seismic data processing methods are used to keep and recover all kinds of attributes applied for fracture detection, which guarantee the high accurate of these attributes, thus then improve the accurate of fracture detection. During seismic data processing, the paper adopted the three dimensional F-Kx-Ky field cone filter technique to attenuate ground roll waves and multiple waves, then enhances the S/N ratio of pre-stack seismic data; comprehensively applying geometrical spread compensation, surface consistent amplitude compensation, residual amplitude compensation to recover amplitude; common azimuth processing method effectively preserves the azimuthal characteristics of P-wave attributes; the technique of bend ray adaptive aperture pre-stack time migration insures to obtain the best image in each azimuth. Application of these processing methods guaranteed these attributes’ accuracy, and then improved the accuracy of fracture detection. After comparing and analyzing a variety of attributes, relative wave impedance (relative amplitude) attribute is selected to inverse the orientation of fracture medium; attenuation gradient and corresponding frequency of 85% energy are selected to inverse the intensity of fracture medium; then obtained the fracture distribution characteristics of lower Paleozoic and Precambrian in ZX ancient buried mountains. The results are good accord with the characteristics of faults system and well information in this area.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.

复杂碎屑岩储层动、静态参数计算理论与方法研究

Reservoir characterization and reservoir modeling are two key techniques in petroleum exploration and development. They all are based on the reliable static and dynamic properties of the reservoirs, especially the static and dynamic properties of the reservoirs at each borehole. Without the static and dynamic properties of the reservoir, reservoir characterization and reservoir modeling will pass into nothingness. In fact, the static and dynamic properties of the reservoir are needed in every domain and stage of petroleum exploration & development Today, petroleum industry has reached a stage worldwide that most of the simple & large massive reservoirs have been well explored and developed. As a result, oil companies are paying more and more attention to the exploration & development of the complex & middle to small clastic reservoirs (such as low resistivity sandstone reservoirs, low or no resistivity contrast sandstone reservoirs, conglomerate reservoirs, volcanoclastic reservoirs). In the recent years, oil companies inside and outside China are focusing on the exploration and development elastic reservoirs. Most of the theories & methods being applicable for simple clastic reservoirs can not be used in complex clasic reservoirs. Some theories & methods that are not resolved in the case of simple clasic reservoirs become more impossible to be resolved in the case of complex elastic reservoirs. A set of theories & methods being applicable for computing the static and dynamic properties of the complex elastic reservoirs are developed in this paper and they have been put into practice successfully. These theories & methods are developed by integrating multi-subjects such as geology, well logging and reservoir engineering, in which geology is used as direction and modern well logging technology is used as basis and reservoir engineering is used as assistance and computer technology is used as tool. There are three outstanding breakthroughs in this paper: of the low porosity fractured and/or vuggy carbonate/igneous reservoirs too. A set of practical theories and methods of computing the static properties (such as porosity, saturation, lithology and fluid type) & dynamic properties (such as permeability and production rate) of simple clastic reservoirs have been developed with the hard efforts of many petroleum engineers and scientists in the past 70 years. However, only some of the theories & methods being applicable for simple clastic reservoirs can be used in complex clastic reservoirs after little modification because of the complexity of the complex clastic reservoirs. Most of the theories & methods being applicable for simple clastic reservoirs can not be used in complex clasic reservoirs. Some theories & methods that are not resolved in the case of simple clasic reservoirs become more impossible to be resolved in the case of complex clastic reservoirs. A set of theories & methods being applicable for computing the static and dynamic properties of the complex clastic reservoirs are developed in this paper and they have been put into practice successfully. These theories & methods are developed by integrating multi-subjects such as geology, well logging and reservoir engineering, in which geology is used as direction and modern well logging technology is used as basis and reservoir engineering is used as assistance and computer technology is used as tool. There are three outstanding breakthroughs in this paper:

泥岩裂缝储层地震识别的理论和方法

During the exploration of fractured reservoirs, worldwide difficult problems will be encountered: how to locate the fractured zones, how to quantitatively determine the azimuth, density, and distribution of the fractures, and how to compute the permeability and porosity of the fractures. In an endeavor to solve these problems, the fractured shale reservoir in SiKou area of ShengLi oil field was chosen as a study area. A study of seismic predictive theory and methods for solving problems encountered in fractured reservoir exploration are examined herein. Building on widely used current fractured reservoir exploration techniques, new seismic theories and methods focusing on wave propagation principles in anisotropic medium are proposed. Additionally, integrated new seismic data acquisition and processing methods are proposed. Based on research and application of RVA and WA methods from earlier research, a new method of acoustic impedance varying with azimuth (IPVA) creatively is put forth. Lastly combining drilling data, well log data, and geologic data, an integrated seismic predictive method for cracked reservoir bed was formed. A summary of the six parts of research work of this paper is outlined below. In part one, conventional geologic and geophysical prediction methods etc. for cracked reservoir exploration are examined, and the weaknesses of these approaches discussed. In part two, seismic wave propagation principles in cracked reservoirs are studied. The wave equation of seismic velocity and attenuation factor in three kinds of fracture mediums is induced, and the azimuth anisotropy of velocity and attenuation in fracture mediums is determined. In part three, building on the research and application of AVA and WA methods by a former researcher, a new method of acoustic impedance creatively varying with azimuth (IPVA) is introduced. A practical software package utilizing this technique is also introduced. In part four, Base on previously discussed theory, first a large full azimuth 3d seismic data (70km~2) was designed and acquired. Next, the volume was processed with conventional processing sequence. Then AVA, WA, and IPVA processing was applied, and finally the azimuth and density of the fractures were quantitatively determined by an integrated method. Predictions were supported by well data that indicate the approach is highly reliable. in part five, geological conditions contributing to cracked reservoir bed formation are analyzed in the LuoJia area resulting in the discovery that the main fractured zones are related to fault distribution in the basin, that also control the accumulation of the oil and gas, the generation mechanisms and types of fractured shale reservoirs are studied. Lastly, by using full 3D seismic attributes, azimuth and density of cracked reservoir zones are successfully quantitative predicted. Using an integrated approach that incorporates seismic, geologic and well log data, the best two fractured oil prospects in LouJia area are proposed. These results herein represent a break through in seismic technology, integrated seismic predictive theory, and production technology for fractured reservoirs. The approach fills a void that can be applied both inside China, and internationally. Importantly, this technique opens a new exploration play in the ShengLi oil field that while difficult has substantial potential. Properly applied, this approach could play an important role toward stabilizing the oil field' production. In addition, this technique could be extended fracture exploration in other oil fields producing substantial economic reward.

火成岩储层裂缝识别、预测理论与方法研究－以济阳坳陷火成岩油藏为例

Since 1970s, igneous reservoirs such as Shang741, Bin674 and Luol51 have been found in Jiyang depression, which are enrichment and heavy-producing. Showing good prospect of exploration and development, igneous reservoirs have been the main part of increasing reserves and production in Shengli oilfield. As fracture igneous reservoir being an extraordinary complex concealed reservoir and showing heavy heterogeneity in spatial distribution, the study of recognition, prediction, formation mechanism and the law of distribution of fracture is essential to develop the reservoir. Guided by multiple discipline theory such as sedimentology, geophysics, mineralogy, petroleum geology, structural geology and reservoir engineering, a set of theories and methods of recognition and prediction of fractured igneous rock reservoir are formed in this paper. Rock data, three-dimensional seismic data, log data, borehole log data, testing data and production data are combined in these methods by the means of computer. Based on the research of igneous rock petrography and reservoir formation mechanism, emphasized on the assessment and forecast of igneous rock reservoir, aimed at establishing a nonhomogeneity quantification model of fractured igneous rock reservoir, the creativity on the fracture recognition, prediction and formation mechanism are achieved. The research result is applied to Jiyang depression, suggestion of exploration and development for fractured igneous rock reservoir is supplied and some great achievement and favourable economic effect are achieved. The main achievements are gained as follows: 1. The main facies models of igneous rock reservoir in JiYang depression are summarized. Based on data and techniques of seism, well log and logging,started from the research of single well rock facies, proceeded by seismic and log facies research, from point to line and line to face, the regional igneous facies models are established. And hypabyssal intrusion allgovite facies model, explosion volcaniclastic rock facies model and overfall basaltic rocks facies model are the main facies models of igneous rock reservoir in JiYang depression. 2. Four nonhomogenous reservoir models of igneous reservoirs are established, which is the base of fracture prediction and recognition. According to characteristics of igneous petrology and spatial types of reservoir, igneous reservoirs of Jiyang depression are divided into four categories: fractured irruptive rock reservoir, fracture-pore thermocontact metamorphic rock and irruptive rock compound reservoir, pore volcanic debris cone reservoir and fracture-pore overfall basaltic rock reservoir. The spatial distribution of each model's reservoir has its features. And reservoirs can be divided into primary ones and secondary ones, whose mechanism of formation and laws of distribution are studied in this paper. 3. Eight geologic factors which dominate igneous reservoirs are presented. The eight geologic factors which dominates igneous reservoirs are igneous facies, epigenetic tectonics deformation, fracture motion, intensity of intrusive effect and adjoining-rock characters, thermo-contact metamorphic rock facies, specific volcano-tectonic position, magmatic cyclicity and epigenetic diagenetic evolution. The interaction of the eight factors forms the four types nonhomogenous reservoir models of igneous reservoirs in Jiyang depression. And igneous facies and fracture motion are the most important and primary factors. 4. Identification patterns of seismic, well log and logging facies of igneous rocks are established. Igneous rocks of Jiyang depression show typical reflecting features on seismic profile. Tabular reflection seismic facies, arc reflection seismic facies and hummocky or mushroom reflection seismic facies are the three main facies. Logging response features of basic basalt and diabase are shown as typical "three low and two high", which means low natural gamma value, low interval transit-time, low neutron porosity, high resistivity and high density. Volcaniclastic rocks show "two high and three low"-high neutron porosity, high interval transit-time, low density, low-resistance and low natural gamma value. Thermo-contact metamorphic rocks surrounding to diabase show "four high and two low" on log data, which is high natural gamma value, high self-potential anomaly, high neutron porosity, high interval transit-time and low density and low-resistance. Based on seismic, well log and logging data, spatial shape of Shang 741 igneous rock is described. 5. The methods of fracture prediction and recognition for fractured igneous reservoir are summarized. Adopting FMI image log and nuclear magnetic resonance log to quantitative analysis of fractured igneous reservoir and according to formation mechanism and shape of fracture, various fractures are recognized, such as high-angle fracture, low-angle fracture, vertical fracture, reticulated fracture, induced fracture, infilling fracture and corrosion vug. Shang 741 intrusive rock reservoir can be divided into pore-vug compound type, pore fracture type, micro-pore and micro-fracture type. Physical properties parameters of the reservoir are computed and single-well fracture model and reservoir parameters model are established. 6. Various comprehensive methods of fracture prediction and recognition for fractured igneous reservoir are put forward. Adopting three-element (igneous facies, fracture motion and rock bending) geologic comprehensive reservoir evaluation technique and deep-shallow unconventional laterolog constrained inversion technique, lateral prediction of fractured reservoir such as Shang 741 is taken and nonhomogeneity quantification models of reservoirs are established.