一种基于模糊petri网专家系统的面向对象方法设计

简要介绍了模糊petri网以及模糊产生式规则,给出了基于模糊petri网的专家系统的框架设计,并提出了模糊产生式规则和模糊petri网的详细设计,根据本设计方案开发了汽车变速箱故障诊断专家系统,证明设计方案简洁高效,扩充性和实用性好。

构造复杂区变速成图方法研究

In the complex structure areas, velocity field building and structure mapping are important for seismic exploration. With the development of seismic exploration, the methods of structure mapping, reservoir prediction and reservoir description all require high precious velocity field. And more accurate depth-structure maps are required for well site design. Aiming at the problems and defects in velocity analysis and structure mapping in oil seismic exploration, the paper which is based on the studies of real data in several areas combines the theories with practical application, and analyzes the precision and applicability of several methods of velocity model building. After that, the following methods are mainly studied: the coherence inversion methods based on the pre-stack CMP gathers or stacking velocity; the interval velocity inversion methods constrained by multi-well; the Random Simulation method; 3D Image Ray Map Migration method and the structure mapping in floating datum and in fixed datum, and then we conclude the method of building high precious seismic velocity field and structure mapping with variable velocity. Firstly, the paper analyses the distributing rule of the velocity variation in the areas with complex structures in the northwest of China, then points out that velocity is a crucial factor which influences the precision of structure mapping, and the velocity variations have something to do with the shapes of the structures, the variety of lithology and so on. The key point of improving the precision of seismic velocity field is to obtain a structure mapping with high precision. We also describe the range and conditions of these methods. Secondly, by comparing many popular methods of velocity model building, we propose a new method in the use of velocity model building. The new method is more effective in velocity model building under every kind of complex condition and is worthy of spreading. At last, the paper fingers out that it is a system engineering to study variable velocity mapping in every kind of complex structure areas. Every step of the work can affect the final results. So it is important to build high efficient and practical velocity model and the flows of mapping processing. The paper builds the flows and gives some examples. The method has been applied in more than ten exploring surveys. The application proves that this method could bring good effect on researching on low-amplitude trap, reservoir prediction, reservoir description and the integrated research of oil&gas geology. Keywords: structure mapping velocity model building complex structure variable velocity media

全信息精细储层预测方法研究与实践-以周清庄油田沙三段油藏为例

With the continuously proceeding of petroleum exploratory development in China, exploratory development becomes more and more difficult. For increasing reserve volume and production, lithologic hydrocarbon reservoir has been the most workable, potential and universality exploration targets. In the past, Dagang Oil Field use the complicated fault reservoir theory as the guide, develop and form a suit of matching construction and instrument in prospecting complicated fault reservoir that reach top of exploration industry in China. But the research of lithologic hydrocarbon reservoir is not much, which affects the exploitation progress of lithologic hydrocarbon reservoir. In this thesis, is object, through the depth study of lithologic deposition in Shasan segment of Zhouqingzhuang Oil Field, a suit of holographic fine reservoir bed forecasting techniques is built up and finally gets following main results: 1. Applying geology, seism, drilling, logging and other information to sensitivity preferences, geological model, inversion and integrated stratum evaluation, realizing the method and flow of refined multi-information stratum forecast. 2. Built up a full three dimensional fine structural interpretation method: in view of r problem of accurately demarcating 90% inclined well, propose a inclined well air space demarcating method, make bed demarcating more exactly; in view of problem of faults demarcating and combination in seismic interpretation, propose a computational method of seismic interference based on wavelet translation, make identify the fault in different level more dependable and reasonable; for exactly identifying structural attitude, propose a velocity modeling method under multi-well restriction, make structural attitude closer to the facts. 3. Built up a high accuracy reservoir bed inversion method: in view of problem in exactly identifying reservoir and nonreservoir with conventional wave impedance inversion method in this place, propose a reservoir log response characteristic analysis and sensible log parameter inversion method. ①analysis log response of reservoir and nonreservoir in region of interest, make definite the most sensible log parameter in identifying reservoir and nonreservoir in this region; ②make sensible log parameter inversion based on wave impedance inversion, to improve inversion accuracy, the thickness of recognizable reservoir bed reach 4-5m. 4. Built up a 4-D reservoir forcasting circuit: in view of difficulty that in lithologic hydrocarbon reservoir making reservoir space characteristic clear by using structural map and reservoir forecasting techniques once only, propose a 4-D reservoir forcasting circuit. In other words, based on development conceptual design, forcast reservoir of different time, namely multiple 3D reservoir forcasting in time queue, each time the accuracy degree of reservoir forcasting is improved since apply the new well material, thereby achieve high quality and highly efficient in exploratory development. During exploratory development lithologic depositin in Shasan segment of Zhouqingzhuang Oil Field, there are thirteen wells get 100% success rate, which sufficiently proves that this suit of method is scientific and effective.

东海陆架盆地复杂天然气藏地震勘探判别技术研究

Exploration study proves that East sea shelf basin embeds abundant hydrocarbon resources. However, the exploration knowledge of this area is very low. Many problems in exploration are encountered here. One of them is that the gas reservoir of this area, with rapid lateral variation, is deeply buried. Correlation of Impendence between sandstone, gas sand and shale is very poor. Another problem is that the S/N ratio of the seismic data is very low and multiples are relatively productive which seriously affect reservoir identification. Resolution of the seismic data reflected from 2500-3000 meter is rather low, which seriously affects the application of hydrocarbon direct identification (HDI) technology. This research established a fine geological & geophysical model based on drilling、well logging、geology&seismic data of East sea Lishui area. A Q value extraction method from seismic data is proposed. With this method, Q value inversion from VSP data and seismic data is performed to determine the subsurface absorption of this area. Then wave propagation and absorption rule are in control. Field acquisition design can be directed. And at the same time, with the optimization of source system, the performance of high resolution seismic acquisition layout system is enhanced. So the firm foundation is ensured for east sea gas reservoir exploration. For solving the multiple and amplitude preserving problems during the seismic data processing, wave equation pre-stack amplitude preservation migration and wave equation feedback iteratively multiple attenuation technologies are developed. Amplitude preservation migration technology can preserve the amplitude of imaging condition and wave-field extrapolation. Multiple removing technology is independent of seismic source wavelet and velocity model, which avoiding the weakness of Delft method. Aiming at the complicated formation condition of the gas reservoir in this area, with dissecting typical hydrocarbon reservoir, a series of pertinent advanced gas reservoir seismic identification technologies such as petrophysical properties analyzing and seismic modeling technology、pre-stack/post-stack joint elastic inversion, attribute extraction technology based on seismic non-stationary signal theory and formation absorption characteristic and so on are studied and developed. Integrated analysis of pre-stack/post-stack seismic data, reservoir information, rock physics and attribute information is performed. And finally, a suit of gas reservoir identification technology is built, according to the geological and geophysical characteristics of this area. With developed innovative technologies, practical application and intergrated interpretation appraisal researches are carried out in Lishui 36-1.The validity of these technologies is tested and verified. Also the hydrocarbon charging possibility and position of those three east sea gas exploration targets are clearly pointed out.

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

西峰油田特低渗透油藏渗流机理及开发对策研究

Micro-pore-throat, micro-fracture and low permeability are the most obvious characters of Xifeng ultra-low permeability reservoir, and threshold pressure gradient and medium deformation during the period of oilfield developing results non-linear seepage feature of the formation liquid flowing in the porous medium underground. It is impossible to solve some problems in the ultra-low permeability reservoir development by current Darcy filtration theory and development techniques. In the view of the characters of ultra-low permeability and powerful-diagenesis and fracture up-growth, the paper quantitatively characterizes of through-going scope for reservoir parameters together with some materials such as similarity field outcrop, rock core, drilling, well logging and production dynamic, which provides geological base for further development adjustment. Based on the displacement experiment of different kinds of seepage fluid and oil-water two phases, this paper proves the relationship between threshold pressure gradient and formation permeability in experiment and theory, which is power function and its index is about -1. The variation rule and the mechanism of oil-water two phases threshold pressure gradient are studied. At the same time, based on the experiment of medium deformation, the variation rule of formation physical property parameters and the deformation mechanism are researched, and the influential factors on the medium deformation are analyzed systematically. With elastic unsteady filtration theory, nonlinear mathematical models of the steady and unsteady flow of single phase as well as horizontal well flow and oil-water two phases flow are deduced with the influence of nonlinear factors including threshold pressure gradient and media deformation. The influences of nonlinear factors upon well deliverability and reservoir pressure distribution as well as the saturation variation pattern of oil-water front are analyzed. By means of the researches such as reasonable well pattern, reasonable well array ration, artificial fracture length optimization advisable water flood timing and feasibility of advanced water flooding, it is necessary to find out effective techniques in order to improve development result of this kind of reservoir. This research result develops and improves on low-velocity nonlinear seepage theory, and offers ways to study similar kind of reservoir; it is meaningful to the development of the ultra-low permeability oil and gas reservoir.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

滑坡泥石流信息遥感获取与空间预测

On the issue of geological hazard evaluation（GHE）, taking remote sensing and GIS systems as experimental environment, assisting with some programming development, this thesis combines multi-knowledges of geo-hazard mechanism, statistic learning, remote sensing (RS), high-spectral recognition, spatial analysis, digital photogrammetry as well as mineralogy, and selects geo-hazard samples from Hong Kong and Three Parallel River region as experimental data, to study two kinds of core questions of GHE, geo-hazard information acquiring and evaluation model. In the aspect of landslide information acquiring by RS, three detailed topics are presented, image enhance for visual interpretation, automatic recognition of landslide as well as quantitative mineral mapping. As to the evaluation model, the latest and powerful data mining method, support vector machine (SVM), is introduced to GHE field, and a serious of comparing experiments are carried out to verify its feasibility and efficiency. Furthermore, this paper proposes a method to forecast the distribution of landslides if rainfall in future is known baseing on historical rainfall and corresponding landslide susceptibility map. The details are as following: (a) Remote sensing image enhancing methods for geo-hazard visual interpretation. The effect of visual interpretation is determined by RS data and image enhancing method, for which the most effective and regular technique is image merge between high-spatial image and multi-spectral image, but there are few researches concerning the merging methods of geo-hazard recognition. By the comparing experimental of six mainstream merging methods and combination of different remote sensing data source, this thesis presents merits of each method ,and qualitatively analyzes the effect of spatial resolution, spectral resolution and time phase on merging image. (b) Automatic recognition of shallow landslide by RS image. The inventory of landslide is the base of landslide forecast and landslide study. If persistent collecting of landslide events, updating the geo-hazard inventory in time, and promoting prediction model incessantly, the accuracy of forecast would be boosted step by step. RS technique is a feasible method to obtain landslide information, which is determined by the feature of geo-hazard distribution. An automatic hierarchical approach is proposed to identify shallow landslides in vegetable region by the combination of multi-spectral RS imagery and DEM derivatives, and the experiment is also drilled to inspect its efficiency. (c) Hazard-causing factors obtaining. Accurate environmental factors are the key to analyze and predict the risk of regional geological hazard. As to predict huge debris flow, the main challenge is still to determine the startup material and its volume in debris flow source region. Exerting the merits of various RS technique, this thesis presents the methods to obtain two important hazard-causing factors, DEM and alteration mineral, and through spatial analysis, finds the relationship between hydrothermal clay alteration minerals and geo-hazards in the arid-hot valleys of Three Parallel Rivers region. (d) Applying support vector machine (SVM) to landslide susceptibility mapping. Introduce the latest and powerful statistical learning theory, SVM, to RGHE. SVM that proved an efficient statistic learning method can deal with two-class and one-class samples, with feature avoiding produce ‘pseudo’ samples. 55 years historical samples in a natural terrain of Hong Kong are used to assess this method, whose susceptibility maps obtained by one-class SVM and two-class SVM are compared to that obtained by logistic regression method. It can conclude that two-class SVM possesses better prediction efficiency than logistic regression and one-class SVM. However, one-class SVM, only requires failed cases, has an advantage over the other two methods as only "failed" case information is usually available in landslide susceptibility mapping. (e) Predicting the distribution of rainfall-induced landslides by time-series analysis. Rainfall is the most dominating factor to bring in landslides. More than 90% losing and casualty by landslides is introduced by rainfall, so predicting landslide sites under certain rainfall is an important geological evaluating issue. With full considering the contribution of stable factors (landslide susceptibility map) and dynamic factors (rainfall), the time-series linear regression analysis between rainfall and landslide risk mapis presented, and experiments based on true samples prove that this method is perfect in natural region of Hong Kong. The following 4 practicable or original findings are obtained: 1) The RS ways to enhance geo-hazards image, automatic recognize shallow landslides, obtain DEM and mineral are studied, and the detailed operating steps are given through examples. The conclusion is practical strongly. 2) The explorative researching about relationship between geo-hazards and alteration mineral in arid-hot valley of Jinshajiang river is presented. Based on standard USGS mineral spectrum, the distribution of hydrothermal alteration mineral is mapped by SAM method. Through statistic analysis between debris flows and hazard-causing factors, the strong correlation between debris flows and clay minerals is found and validated. 3) Applying SVM theory (especially one-class SVM theory) to the landslide susceptibility mapping and system evaluation for its performance is also carried out, which proves that advantages of SVM in this field. 4) Establishing time-serial prediction method for rainfall induced landslide distribution. In a natural study area, the distribution of landslides induced by a storm is predicted successfully under a real maximum 24h rainfall based on the regression between 4 historical storms and corresponding landslides.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

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

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

地震波散射非线性反演问题不动点理论及方法的研究

How to create a new method to solve the problem or reduce the influence of that the result of the seismic waves scattering nonlinear inversion is not uniqueness is a main purpose of this research work in the paper. On the background of research into the seismic inversion, new progress of the nonlinear inversion is introduced at the first chapter in this paper. Especially, the development, basic theories and assumptions on some major theories of seismic inversion are analyzed, discussed and summarized in mathematics and physics. Also, the problems faced by the mathematical basis of investigations of the seismic inversion are discussed, and inverse questions of strongly seismic scattering due to strong heterogeneous media in the Earth interior are analyzed and viewed. What the kernel of paper is that gathers all our attention making a new nonlinear inversion method of seismic scattering. The paper provides a theory and method of how to introduce the fixed-point theory into the nonlinear seismic scattering inversion and how to obtain the solution, and gives the actually method to create a serials of contractive mappings of velocity parameter's in the mapping space of wave. Therefore, the results testify the existence of fixed point of velocity parameter and give the method the find it. Further, the paper proves the conclusion that the value obtained by taking the fixed point of velocity parameter into wave equation is the fixed point of the wave of the contractive mapping. Thence, the fixed point is the global minima since the stabilities quality of the fixed point. Based on the new theory, in the chapter three, many inverse results are obtained in the numerical value test. By analysis the results one could find a basic facts that all the results, which are inversed by the different initial model, are tended to the true value in theoretical true model. In other words, the new method mostly eliminates the non-uniqueness that which is existed in seismic waves scattering nonlinear inversion in degree. But, since the test results are quite finite now, more test is need here to positive our theory. As a new theoretical method, it must be existed many weaken in it. The chapter four points out all the questions which is bother us. We hope more people to join us to solve the problem together.

中国南方热带地区田洋玛珥湖24万年以来的气候变化及湖泊演变

A continuous long (224m) and high-resolution core TY2 was recovered from paleo-maar-lake Tianyang, tropical South China. Based on the diatom records of the upper 130-m core, this paper focuses on exploring climate change and the lake evolution history in tropical South China during the past 240ka. The most typical and unique characteristics of the diatom assemblages is that, Aulacoseira granulata was dominant or absolutely dominant species (80-90%) during most parts of the 130-m core, while Cyclotella stelligera var. tenuis and Fragilaria construens var. venter were subdominant species in only limited parts of the lower and upper core, respectively. Time scale is always the biggest problem for the study of TY2 core, so although diatom is seldom used for establishing time scale, here we attempt this by correlating the diatom-reconstructed temperature sequence with the time scale of ODP core 806B from Equatorial Western Pacific. Verified by the few most reliable ages from TY2 core and the parallel core TYl, a rather reasonable and reliable time scale was established. 01S 7/6 falls at the depth of 100m (ca. 194kaBP), OIS 6/5 at 75m (ca. 132kaBP), OIS 5/4 at 46m (ca. 75kaBP), OIS 4-3 at 35m (ca. 60kaBP). Qualitative and quantitative environmental reconstructions are made on the basis of diatom assemblage ecotype and EDDI dataset. Correlation of diatom-reconstructed temperature and moisture changes of Core TY2 with pollen-reconstructed temperature and rainfall sequence of Core TYl proves that the results are quite consistent in most periods. Thus the reconstruction results from diatom are quite reliable, and probably have a much higher resolution than pollen results. Combined with lithological and magnetic susceptibility variations, the diatom analysis reveals that, the general climate in tropical South China during the past 240ka was warm and wet. On the time scale of glacial-interglacial, warm and wet, cool and dry are not always synchronous. It was relatively warm-wet during the penultimate interglacial, cool-dry during the penultimate glacial, warm-dry during the last interglacial, and cooler-drier during the last glacial. In contrast, on the time scale of subglacial-subinterglacial scales, warm and dry, cool and wet corresponds very obviously, showing very clear 21-23 ka precession cycle. Analysis also shows that, the water of Tianyang paleo-maar-lake was generally warm, turbulent, turbid, meso-trophic, slightly alkaline, low conductivity and fresh during the past 240 kaBP, with small variations in some parts. Tianyang paleolake experienced shallow to semi-deep lake in OIS7d, open shallow lake in OIS7c-OIS5b, shallow coastal lake in OIS5a-OIS4c, swamp in OIS4b, and then completely dried up in OIS3c. The lake evolution was mainly controlled by temperature and precipitation changes in tropical China. While temperature and precipitation changes were probably controlled by the migrations of monsoon rainband and the evaporation rate, which was in turn controlled by the evolution of East Asian monsoon. Therefore, when the summer monsoon was strongest the climate was warm-dry, when stronger the climate was warm-wet; when the winter monsoon was strongest the climate was cool-dry, stronger cool-wet. This mechanism caused the warm-dry sub interglacial and cool-wet subglacial climate in the tropical South China.

水库滑坡复活机理及其发展趋势预测研究

Various numbers of ancient landslides of various scales are frequently distributed on both banks of reservoirs, especially large reservoirs, both in China and abroad. During inundation and operation of theses reservoirs, some of the landslides are reactivated, which caused losses of people's lives and properties to various extents, some even disasters. Systematic studies are, however, very few on the reservoir-induced reactivation mechanism and development tendency prediction. Based on investigation of reservoir-induced reactivation phenomena of ancient landslides and relevant existing research problems, a systematic study is carried out on the field identification, induced reactivation mechanism, development tendency prediction, risk decision-making and treatment of reservoir-related ancient landslides, through analysis of large numbers of engineering geological investigation results, scientific experimental and research results, in combination with prevention and treatment practices of reservoir-related landslides both in China and abroad, and a series of research results have been obtained. 1. On the basis of study of the distribution features, genesis mechanism of ancient landslides on river banks, a set of scientific methods are summarized on field identification of ancient landslides, and a significant method named "lithologic sequence method" or "indicator layer method", is proposed, which is proved to be very useful. 2. A detail study is made on the reservoir-induced hydraulic effects and material mechanic effects (or softening effects) on the ancient landslide through model and case studies, which concludes that the magnitude and properties of reservoir-induced hydraulic effects are related to the shapes of sliding planes, water content and permeability of landslide materials and variation rate and magnitude of reservoir levels; the magnitude of material mechanic effects are related to the material composition (including mineral composition and grain size), natural water content and saturation state of sliding zones. Also a sensitive analysis is made on the factors that are related to the stabilities of the landslides, which indicate that the stability of a landslide is more sensitive to the groundwater head h_w in the slides and the inner friction angleψof sliding zones than others. 3. The joint inducing mechanism of rainfall and reservoir is also discussed in the paper through model analysis and case study, which proves that reservoir inundation increases firstly the sensitivity of a landslides to rainfall through reduction of its stability or cracking deformation which will increase the rainfall infiltration to the slide body, and then rainfall triggers reactivation or intensifies the reservoir-induced deformation of a landslide. 4. Based on rheologic test results of sliding zones of several reservoir-related ancient landslides, the rheologic characteristics of sliding zones have been discussed in detail and several typical rheologic models have been set up, which well explains the dynamic process of slide deformation. The response types to reservoir inundation and development tendency of reservoir -related ancient landslide are discussed in the paper based on field investigation results. And prediction methods for reservoir-related landslides have been studied based on the Mate-Synthetic principle of quantitative and qualitative analysis, as well as combination of computation and internal mechanism analysis, and a rheologic analytical method is proposed which is proved very useful for prediction of the landslide development tendency. 6. In disaster-prevention and treatment of reservoir-related landslides, risk decision-making has been proved very significant both in engineering and economics. Based on the practices in disaster-prevention and treatment of reservoir-related landslides both in China and abroad, the disaster-prevention risk decision-making for reservoir-related landslides has been proposed in terms of philosophy, methods and procedures, and well put into practice. A summary is also made through case study of the experiences of treatment of reservoir-related landslides both in China and abroad in terms of principle, methods and technical lines. 7 A detail study is made as a case study of the reactivated Maoping ancient landslide on the left bank of the Geheyan Reservoir on Qingjiang river in Hubei province, China, including its field identification features, reservoir-induced reactivation characteristics and mechanism, development tendency prediction and proposed counter measures based on risk analysis.

东营凹陷沙三中亚段湖相浊积岩储层特征及油气成藏规律研究

Turbidity sandstone reservoirs have been an important field of hydrocarbon exploration and development in the basins all over the world, as well as in China. Lithologic pools are composed of turbidity sandstones and other sandstones are frequently found in the Jiyang Depression that is a Mesozoic-Cenozoic non-marine oil-bearing basin. The Dongying Sag lies in the sedimentary center of the basin. The subtle traps with turbidity reservoirs are generally difficult to be predicted and described by using current techniques. The studies on turbidity reservoirs plays thus an important theoretical and theoretical practical role in exploration and development in the Jiyang Depression. The attention is, in this thesis, focused on the petrologic properties and oil accumulating behaviors in lake turbidity sedimentary systems in the middle part of the third section of Shahejie Formation in the Dongying Sag, especially in Dongxin area, which lies on the central uplift of the Sag. The paper has disclosed the origin types of turbidity sandstones, distribution pattern and controlling factors of turbidity sandstones, and set up hydrocarbon accumulation patterns of the middle part of the third section of Shahejie Formation in Dongxin, based on nonmarine high resolution sequence stratigraphy, event sedimentology and new theories of hydrocarbon forming. By studying prediction method and technology of turbidity sandstone reservoirs, using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, the paper has forecast low permeability turbidity sandstone reservoirs and pointed out advantage exploration aims to progressive exploration and development. The paper has obtained mainly many productions and acknowledges as follows: 1.Turbidity sandstone reservoirs of the third section of Shahejie Formationin Dongying Sag are formed in such specifical geological background as rift and extension of basin. The inherited Dongying delta and transgression make up many turbidity distribution areas by overlaying and joining together. The hydrocarbon migrates from depression area to adjacent turbidity sandstone continuously. Accumulation area which is sufficient in oil is formed. 2.The paper has confirmed distinguishable sign of sequence boundary , established stratigraphic framework of Dongying Sag and realized isotime stratigraphic correlation. Es3 of Dongying delta is divided into eleven stages. Among them, the second period of the lower section in Es3, the sixth period of the middle section in Es3, the third period of the upper section in Es3 correspond to eleven sedimentary isotime surface in seismic profile, namely Es3 is classified into eleven Formations. 3.According to such the features of turbidity sandstone as deep in burial, small in area, strong in subtle property, overlaying and joining together and occurring in groups, management through fault and space variations of restriction quantum are realized and the forecast precision of turbidity sandstone by using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, based on the analysis of all kinds of interwell seismic inversion techniques. 4.According to the features of low permeable turbidity sandstone reservoirs, new method of log interpretation model is put forward. At the same time, distinguish technology of familiar low resistivity oil layer in the turbidity sandstone reservoirs is studied based on petrophysical laboratory work and "four properties" interrelationship between lithological physical Jogging and bearing hydrocarbon properties. Log interpretation model and reservoir index interpretation model of low resistivity oil layer are set up. So the log interpretation precision is improved. 5.The evolution law and its difference of the turbidity sandstone are embodies as follows: the source of sediments come from the south and east of the study area in the middle period of Es3. East source of sediments is pushed from west to east. However, the south source supply of sediments in the early and middle period of Es3 is in full, especially in Es3. subsequently, the supply is decreased gradually. Turbidity fan moves back toward the south and the size of fan is minished accordingly. The characteristic of turbidity sandstone in Dongying Sag is different in different structural positions. Dongxin in the middle-east of the central lift and Niuzhuang Sag He in Dongying delta front and prodelta deep lake subfacies. Although the turbidity sandstone of the two areas root in the Dongying delta sedimentary system, the sand body has different remarkably characteristic. 6.The sedimentary model of the turbiditys in study area have three types as follows: (1) collapse turbidity fan in respect of delta; (2) fault trench turbidity fan; (3) other types of microturbidity sandstone. Middle fan and outer fan, can be found mainly in sublacustrine fan. Middle fan includes braided channel microfacies, central microfacies and braided interchannel microfacies, which is main prospecting oil-bearing subfacies. The middle section of the third section of Shahejie Formation in study area (for example the central lift) can be divided into middle-lower and upper part. The middle-lower part is characteristic of turbidity fan. The upper part is sedimented mainly by delta-collapse fan. 7.The turbidity reservoirs of the middle part of the third section of Shahejie Formation in study area characterize by low maturity both in component and texture, strong in diagenesis and low in permeability. The reservoir can be classified into four types. Type III is the body of reservoir and comprises two types of H a and HI b. M a belongs to middle porosity - low permeability reservoir and distributes in the central lift. Hlb belongs to low porosity - low permeability and distributes in Haojia region. 8.A11 single sand body of lens turbidity reservoir of the middle part of the third section of Shahejie Formation in study area are surrounded by thick dark source rocks. The oil-water system is complex and behaves that every sandstone is single seal unit. The water body is 1/3-1-5 of the sand body. The edge water is not active. The gas exists in the top of reservoir in the form of mixed gas. For far-range turbidity fan with big scale channel, the area and volume of sand body is large and the gap is big in oil packing degree. There are lots of edge water and bottom water, and the latter increases rapidly during the course of development. 9.By exerting the modern hydrocarbon forming theories, the third section of Shahejie Formation in study area belongs to abnormally pressured fluid compartment. The lithological reservoir of the third section of Shahejie Formation is formed in the compartment. The reservoir-formed dynamic system belongs to lower self-source enclosed type. The result and the practice indicate that the form and accumulation of lithological oil reservoirs are controlled by the temperature and pressure of stratum, microfacies, thickness of sand body, fault and reservoir heterogeneity. 10. Based on studies above, the emphases focus on in south and north part of Dongying structure, west Dongxin region and south part Xinzhen structure in the application of production. The practice proves that the turbidity sandstone reservoirs in Ying 11 block and the fault-lithological reservoirs in Xin 133 block have been obtained significant breakthrough. The next target is still sandstone groups of the third section of Shahejie Formation in the bordering areas of Dongxin region for instance Xin 149 area, He 89 area, Ying 8 area etc.

利用分频与反演提高地震资料信噪比与分辨率方法研究

This thesis mainly talks about the wavelet transfrom and the frequency division method. It describes the frequency division processing on prestack or post-stack seismic data and application of inversion noise attenuation, frequency division residual static correction and high resolution data in reservoir inversion. This thesis not only describes the frequency division and inversion in theory, but also proves it by model calculation. All the methods are integrated together. The actual data processing demonstrates the applying results. This thesis analyzes the differences and limitation between t-x prediction filter and f-x prediction filter noise attenuation from wavelet transform theory. It considers that we can do the frequency division attenuation process of noise and signal by wavelet frequency division theory according to the differences of noise and signal in phase, amplitude and frequency. By comparison with the f-x coherence noise, removal method, it approves the effects and practicability of frequency division in coherence and random noise isolation. In order to solve the side effects in non-noise area, we: take the area constraint method and only apply the frequency division processing in the noise area. So it can solve the problem of low frequency loss in non-noise area. The residual moveout differences in seismic data processing have a great effect on stack image and resolutions. Different frequency components have different residual moveout differences. The frequency division residual static correction realizes the frequency division and the calculation of residual correction magnitude. It also solves the problems of different residual correction magnitude in different frequency and protects the high frequency information in data. By actual data processing, we can get good results in phase residual moveout differences elimination of pre-stack data, stack image quality and improvement of data resolution. This thesis analyses the characters of the random noises and its descriptions in time domain and frequency domain. Furthermore it gives the inversion prediction solution methods and realizes the frequency division inversion attenuation of the random noise. By the analysis of results of the actual data processing, we show that the noise removed by inversion has its own advantages. By analyzing parameter's about resolution and technology of high resolution data processing, this thesis describes the relations between frequency domain and resolution, parameters about resolution and methods to increase resolution. It also gives the processing flows of the high resolution data; the effect and influence of reservoir inversion caused by high resolution data. Finally it proves the accuracy and precision of the reservoir inversion results. The research results of this thesis reveal that frequency division noise attenuation, frequency residual correction and inversion noise attenuation are effective methods to increase the SNR and resolution of seismic data.