多代理技术在船体装配CAPP系统的应用

本文介绍了多代理技术在船体装配 CAPP系统中应用 .包括 :船体装配 CAPP系统的结构 ,Agent之间的通讯 ,Agent表达方法、提供的函数功能描述 ,以及各工艺代理的功能描述和封装实例 ;采用多代理技术的船体装配 CAPP系统可以生成满足工厂要求的装配工艺文件、工艺定额和材料定额

CIMS设计过程中存在的一些问题及解决方案

提出了一种新的集成的I2 DEF方法 ,并介绍了与之配套的设计开发工具 ,用以支持大型复杂信息系统的设计与开发 ,它可以成功地解决计算机集成制造系统设计开发过程中遇到的许多问题。本文指出了我国CIMS工程存在的问题 ,分析了这些问题产生的原因 ,并结合企业实际给出了应用I2 DEF方法的解决方案。

对象互操作的行为模型

对象互操作表达了一组对象在完成某一任务时的动态协作关系，对象互操作的行为描述与抽象是支持面向应用对象互操作的基础．对此，提出一种活动模型作为描述对象互操作行为的方法．该方法以一阶时态逻辑为基础，表达了互操作对象之间交换消息的时态顺序和不同活动之间的行为关系．在该方法中，提出了活动特化和活动聚合两种行为抽象机制，实现了对象互操作行为的复用．最后讨论了给定论域的类模式和活动模式的一致性集成问题。

地质导向钻井关键技术及方法

On the basis of analyzing the principle and realization of geo-steering drilling system, the key technologies and methods in it are systematically studied in this paper. In order to recognize lithology, distinguish stratum and track reservoirs, the techniques of MWD and data process about natural gamma, resistivity, inductive density and porosity are researched. The methods for pre-processing and standardizing MWD data and for converting geological data in directional and horizontal drilling are discussed, consequently the methods of data conversion between MD and TVD and those of formation description and adjacent well contrast are proposed. Researching the method of identifying sub-layer yields the techniques of single well explanation, multi-well evaluation and oil reservoir description. Using the extremum and variance clustering analysis realizes logging phase analysis and stratum subdivision and explanation, which provides a theoretical method and lays a technical basis for tracing oil reservoirs and achieving geo-steering drilling. Researching the technique for exploring the reservoir top with a holdup section provides a planning method of wellpath control scheme to trace oil and gas reservoir dynamically, which solves the problem of how to control well trajectory on condition that the layer’s TVD is uncertain. The control scheme and planning method of well path for meeting the demands of target hitting, soft landing and continuous steering respectively provide the technological guarantee to land safely and drill successfully for horizontal, extended-reach and multi-target wells. The integrative design and control technologies are researched based on geology, reservoir and drilling considering reservoir disclosing ratio as a primary index, and the methods for planning and control optimum wellpath under multi-target restriction, thus which lets the target wellpath lie the favorite position in oil reservoir during the process of geo-steering drilling. The BHA (bottomhole assembly) mechanical model is discussed using the finite element method, and the BHA design methods are given on the basis of mechanical analyses according to the shape of well trajectory and the characteristics of BHA’s structure and deformation. The methods for predicting the deflection rate of bent housing motors and designing their assemblies are proposed based on the principle of minimum potential energy, which can clearly show the relation between the BHA’s structure parameters and deflection rate, especially the key factors’ effect to the deflection rate. Moreover, the interaction model between bit and formation is discussed through the process of equivalent formation and equivalent bit considering the formation anisotropy and bit anisotropy on the basis of analyzing the influence factors of well trajectory. Accordingly, the inherence relationship among well trajectory, formation, bit and drilling direction is revealed, which lays the theory basis and technique for predicting and controlling well trajectory.

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

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

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

In 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

海拉尔盆地贝尔凹陷基岩潜山特征及储层预测

The lithology of the buried hill of Triassic Budate group in Beier depression is epimetamorphic clastic rock and volcanic clastic rock stratum. Recently the favorable hydrocarbon show was discovered in buried hill of base rock, and large-duty industrial oil stream was obtained in some wells in Beier depression. Based on the information of seismos and wells, the tectonic framework, tectonic deformation times and faulted system of the Beier depression are comprehensively studied, then configuration, evolutional history, genetic type and distributed regularity of buried hill are defined. According to observing description and analysis of core sample, well logging and interpretive result of FMI, the lithological component, diagenetic type and diagenetic sequence of buried hill reservoir are confirmed, then reservoir space system of buried hill is distinguished, and vegetal feature, genetic mechanism and distributed regularity of buried hill fissure are researched, at the same time the quantitative relationship is build up between core fissures and fissures interpreted by FMI. After that fundamental supervisory action of fault is defined to the vegetal degree of fissure, and the fissure beneficial places are forecasted using fractal theory and approach. At last the beneficial areas of Budate group reservoir are forecasted by reservoir appraisal parameters optimization such as multivariate gradually regression analysis et. al. and reservoir comprehensive appraisal method such as weighing analyze and clustering procedure et. al. which can provide foundation for the next exploratory disposition. Such production and knowledge are obtained in this text as those: 1. Four structural layers and two faulting systems are developed, and four structural layers are carved up by three bed succession boundary surfaces which creates three tectonic distortional times homology. Three types of buried hill are divided, they are ancient physiognomy buried hill, epigenetic buried hill, and contemporaneous buried hill. 2. Reservoir space of Budate buried hill is mainly secondary pore space and fissure, which distributes near the unconformity and/or inside buried hill in sections. The buried hill reservoir experienced multi-type and multi-stage diagenetic reconstruction, which led to the original porosity disappeared, and multi secondary porosity was created by dissolution, superficial clastation and cataclasis et. al. in diagenetic stage, which including middle crystal pore, inter crystal pore, moldic pore, inter particle emposieu, corrosion pore space and fissure et. al. which improved distinctly reservoir capability of buried hill. 3. The inner reservoir of buried hill in Beier depression is not stratigraphic bedded construction, but is fissure developing place formed by inner fault and broken lithogenetic belt. The fissures in inner reservoir of buried hill are developed unequally with many fissure types, which mainly are high angle fissure and dictyonal fissures and its developing degree and distribution is chiefly controlled by faulting. 4. The results of reservoir comprehensive evaluate and reservoir predicting indicates that advantageous areas of reservoir of buried hill chiefly distributes in Sudeerte, Beixi and Huoduomoer, which comprehensive evaluate mainly Ⅱand Ⅲ type reservoir. The clues and results of this text have directive significance for understanding the hydrocarbon reservoir condition of buried hill in Beier depression, for studying hydrocarbon accumulated mechanism and distributed regularity, and for guiding oil and gas exploration. The results of this text also can enrich and improve nonmarine hydrocarbon accumulated theory.

渤海海域河流相油藏精细描述技术及应用研究

Aiming at the character of Bohaii Sea area and the heterogeneity of fluvial facies reservoir, litho-geophysics experiments and integrated research of geophysical technologies are carried out. To deal with practical problems in oil fields of Bohai area, such as QHD32-6, Southern BZ25-1 and NP35-2 et al., technology of reservoir description based on seismic data and reservoir geophysical methods is built. In this dissertation, three points are emphasized: ①the integration of multidiscipline; ②the application of new methods and technologies; ③the integration of quiescent and dynamic data. At last, research of geology modeling and reservoir numerical simulation based on geophysical data are integrated. There are several innovative results and conclusion in this dissertation: (1)To deal with problems in shallow sea area where seismic data is the key data, a set of technologies for fine reservoir description based on seismic data in Bohai Sea area are built. All these technologies, including technologies of stratigraphic classification, sedimentary facies identification, structure fine characterization, reservoir description, fluid recognition and integration of geological modeling& reservoir numerical simulation, play an important role in the hydrocarbon exploration and development. In the research of lithology and hydrocarbon-bearing condition, petrophysical experiment is carried out. Outdoors inspection and experiment test data are integrated in seismic forward modeling& inversion research. Through the research, the seismic reflection rules of fluid in porosity are generated. Based on all the above research, seismic data is used to classify rock association, identify sedimentary facies belts and recognition hydrocarbon-bearing condition of reservoir. In this research, the geological meaning of geophysical information is more clear and the ambiguity of geophysical information is efficiently reduced, so the reliability in hydrocarbon forecasting is improved. The methods of multi-scales are developed in microfacies research aiming at the condition of shallow sea area in Bohai Sea: ① make the transformation from seismic information to sedimentary facies reality by discriminant analysis; ②in research of planar sedimentary facies, make microfacies research on seismic scale by technologies integration of seismic multi-attributes analysis& optimization, strata slicing and seismic waveform classification; ③descript the sedimentary facies distribution on scales below seismic resolution with the method of stochastic modeling. In the research of geological modeling and reservoir numerical simulation, the way of bilateral iteration between modeling and numerical simulation is carried out in the geological model correction. This process include several steps: ①make seismic forward modeling based on the reservoir numerical simulation results and geological models; ②get trend residual of forward modeling and real seismic data; ③make dynamic correction of the model according to the above trend residual. The modern integration technology of reservoir fine description research in Bohai Sea area, which is developed in this dissertation, is successfully used in (1)the reserve volume evaluation and development research in BZ25-1 oil field and (2)the tracing while drilling research in QHD32-6 oil field. These application researches show wide application potential in hydrocarbon exploration and development research in other oil fields.

叠前叠后联合反演技术及其在QHD32-6油田的应用研究

What geophysical inversion studied includes the common mathematics physical property of inversion and the constitution and appraisal method of solution in geophysics domain, i.e. using observed physical phenomenon from the earth surface to infer space changing and physical property structure of medium within the earth. Seismic inversion is a branch of geophysical inversion. The basic purpose of seismic inversion is to utilizing seismic wave propagating law in the medium underground to infer stratum structure and space distribution of physical property according to data acquisition, processing and interpretation, and then offer the vital foundation for exploratory development. Poststack inversion is convenient and swift, its acoustic impedance inversion product can reflect reservoir interior changing rule to a certain degree, but poststack data lack abundant amplitude and travel time information included in prestack data because of multiple superimpose and weaken the sensitiveness which reflecting reservoir property. Compared with poststack seismic inversion, prestack seismic inversion has better fidelity and more adequate information. Prestack seismic inversion, including waveform inversion, not only suitable for thin strata physical property inversion, it can also inverse reservoir oil-bearing ability. Prestack seismic inversion and prestack elastic impedance inversion maintain avo information, sufficiently applying seismic gathering data with different incident angle, partial angle stack, gradient and intercept seismic data cube. Prestack inversion and poststack inversion technology were studied in this dissertation. A joint inversion method which synthesize prestack elastic wave waveform inversion, prestack elastic impedance inversion and poststack inversion was proposed by making fully use of prestack inversion multiple information and relatively fast and steady characteristic of poststack inversion. Using the proposed method to extract rock physics attribute cube with clear physical significance and reflecting reservoir characterization, such as P-wave and S-wave impedance, P-wave and S-wave velocity, velocity ratio, density, Poisson ratio and Lame’s constant. Regarding loose sand reservoir in lower member of Minghuazhen formation, 32-6 south districts in Qinhuangdao,as the research object, be aimed at the different between shallow layer loose sand and deep layer tight sand, first of all, acquire physical property parameters suitable for this kind of heavy oil pool according to experimental study, establishing initial pressure and shear wave relational model; Afterwards, performing prestack elastic wave forward and inversion research, summarizing rules under the guidance of theoretical research and numerical simulation, performing elastic impedance inversion, calculating rock physics attributes; Finally, predicting sand body distribution according to rock physics parameters, and predicting favorable oil area combine well-logging materials and made good results.

松辽盆地长岭断陷火山岩储层识别与勘探目标评价

Changling fault depression is the biggest fault subsidence in south of Songliao Basin. In its Lower Cretaceous Yingcheng and Shahezi formations developed thick source rocks of deep lake facies and developed poly-phase volcanic rock reservoirs as well. In recent years, significant breakthroughs have been obtained in hydrocarbon exploration of volcanic rock reservoir in the different fault depressions in Songliao basin. Lately, I have been involved in hydrocarbon exploration in the Changling rift depression, especially volcanic rock reservoirs and exploration targets research, participating in the deployment of well Yaoshen 1 which gained over 40 × 104m3 natural gas flow. As quick changes of lithology and facies in Changling area in the south of Songliao basin, and the volcanic rock interludes distribution in continental clastic rock and shale in 3D space, so the identification of volcanic rock types and distribution become a difficult problem. Thus, based on the integrated research of the wild outcrop observation, gravity, magnetic and seismic data, geophysical logging, drilling and coring, laboratory test, this paper carried out the reservoir identification, description and prediction of volcanic rocks in Changling fault depression. In this area, this paper analyzed the volcanic rocks litho-facies, the eruption period, and characteristics of cycles. At the same time, tried to know how to use logging, seismic data to separate volcanic rocks from sandstone and shale, distinguish between volcanic reservoir and non-reservoir, distinguish between intermediate-basic and acidic volcanic rocks, and how to identify traps of volcanic rocks and its gas-bearing properties, etc. Also it is summarized forming conditions and distribution of traps, and possible gas-bearing traps were optimized queuing management. Conclusions as follows: There are two faulted basements in Changling fault depression, granite basement in the southeast and upper paleozoic epimetamorphic basement in the northwest. The main volcanic reservoirs developed in Yingcheng period, which was the intermediate-basic and acidic volcanic eruptions, from the south to north by the intermediate-basic to acid conversion. The volcanic vents are gradually young from south to north. According to information of the re-processing 3D seismic data and gravity-magnetic data, the large volcanic vent or conduit was mainly beaded-distributed along the main fault. The volcanic rocks thickness in Yingcheng formation was changed by the deep faults and basement boundary line. Compared with the clastic rocks, volcanic rocks in Changling area are with high resistance and velocity (4900-5800), abnormal Gamma. All kinds of volcanic rocks are with abnormal strong amplitude reflection on the seismic stacked section except tuff. By analyzing the seismic facies characteristics of volcanic rocks, optimizing seismic attributes constrained by logging, using seismic amplitude and waveforms and other attributes divided volcanic rocks of Yingcheng formation into four seismic zones in map. Currently, most volcanic gas reservoirs are fault-anticline and fault-nose structure. But the volcanic dome lithologic gas reservoirs with large quantity and size are the main gas reservoir types to be found.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

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

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

孤北洼陷沙三-沙四段隐蔽油藏形成机制及其预测

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in Gubei depression, this paper builds high resolution sequence stratigraphic structure, sedimentary system, sandbody distribution, the effect of tectonic in sequence and sedimentary system evolution and model of tectonic-lithofacies. The pool formation mechanism of subtle trap is developed. There are some conclusions and views as follows. 1.With the synthetic sequence analysis of drilling, seismic, and well log, the highly resolution sequence structure is build in Gubei depression. They are divided two secondary sequences and seven three-order sequences in Shahejie formation. They are include 4 kinds of system traces and 7 kinds of sedimentary systems which are alluvial fan, under water fan, alluvial fan and fan-delta, fan-delta, lacustrine-fan, fluvial-delta-turbidite, lakeshore beach and bar, and deep lake system. Sandbody distribution is show base on third order sequence. 2.Based on a lot of experiment and well log, it is point out that there are many types of pore in reservoir with the styles of corrosion pore, weak cementing, matrix cementing, impure filling, and 7 kinds of diagenetic facies. These reservoirs are evaluated by lateral and profile characteristics of diagenetic facies and reservoir properties. 3.The effect of simultaneous faulting on sediment process is analyzed from abrupt slope, gentle slope, and hollow zone. The 4 kinds of tectonic lithofacies models are developed in several periods in Gubei depression; the regional distribution of subtle trap is predicted by hydro accumulation characteristics of different tectonic lithofacies. 4.There are 4 types of compacting process, which are normal compaction, abnormal high pressure, abnormal low pressure and complex abnormal pressure. The domain type is normal compaction that locates any area of depression, but normal high pressure is located only deep hollow zone (depth more than 3000m), abnormal low pressures are located gentle slope and faulted abrupt slope (depth between 1200~2500m). 5.Two types dynamic systems of pool formation (enclosed and partly enclosed system) are recognized. They are composed by which source rocks are from Es3 and Es4, cap rocks are deep lacustrine shale of Esl and Es3, and sandstone reservoirs are 7 kinds of sedimentary system in Es3 and Es4. According to theory of petroleum system, two petroleum systems are divided in Es3 and Es4 of Gubei depression, which are high or normal pressure self-source system and normal or low pressure external-source system. 6.There are 3 kinds of combination model of pool formation, the first is litholgical pool of inner depression (high or normal pressure self-source type), the second is fault block or fault nose pool in marginal of depression (normal type), the third is fault block-lithological pool of central low lifted block (high or normal pressure type). The lithological pool is located central of depression, other pool are located gentle or abrupt slope that are controlled by lithological, faulting, unconfirmed. 7.This paper raise a new technique and process of exploration subtle trap which include geological modeling, coring description and logging recognition, and well log constrained inversion. These are composed to method and theory of predicting subtle trap. Application these methods and techniques, 6 hydro objects are predicted in three zone of depression.