Boundary Characterization Experiment Series Overview

Ocean acoustic propagation and reverberation in continental shelf regions is often controlled by the seabed and sea surface boundaries. A series of three multi-national and multi-disciplinary experiments was conducted between 2000-2002 to identify and measure key ocean boundary characteristics. The frequency range of interest was nominally 500-5000 Hz with the main focus on the seabed, which is generally considered as the boundary of greatest importance and least understood. Two of the experiments were conducted in the Mediterranean in the Strait of Sicily and one experiment in the North Atlantic with sites on the outer New Jersey Shelf (STRATAFORM area) and on the Scotian Shelf. Measurements included seabed reflection, seabed, surface, and biologic scattering, propagation, reverberation, and ambient noise along with supporting oceanographic, geologic, and geophysical data. This paper is primarily intended to provide an overview of the experiments and the strategies that linked the various measurements together, with detailed experiment results contained in various papers in this volume and other sources

可变误差多面体法同时测定镧及其它轻稀土合量

电子计算机应用于分光光度法为同时分析相互干扰的多组分体系开辟了一个新的研究领域。文献报道的线性规划、因素轮换优选、因子分析、正交分解及可变误差多面体法已显示了计算分光光度法的优越性。在所报道的工作中,测定某单组分的吸光常数时,大都没有其它组分共存,依据α=A/bc来计算其吸光常数。但结果产生较大的误差。本文以化学性质极其相似的La、Ce、Pr、Nd、Sm为分析对象,采用五个已知组成的标样,通过全选主元的Gauss���去法求得单组分的吸光常数。据此再用可变误差多面体法求解样品中各组分的浓度。实践证明,此法不仅减小了实验误差和手工计算时间,而且考虑了组分间的相互作用,所以获得比较准确的结果。可变误差多面体法也是解决这类问题更为合适的方法。

基于测角的自主移动机器人定位算法

给出了自主移动机器人定位的两种算法：解析算法和数值算法。解析法公式较以往的简洁。数值算法结合解析法和高斯－牛顿算法，不仅能避免因初值选取不合理而导致求解过程发散的问题，而且能提高运算精度和速度，通过对两种算法的计算机仿真，表明了解析算法具有运算速度快，而数值算法精度高的特点。其结果已用于自主移动机器人的研制中。

基于弹性参数分析的碳酸盐岩储层预测方法研究

Prediction of Carbonate Reservoir Based on the Elastic Parameter Analysis Zhang Guangzhi (Solid Geophysics) Directed by Professor Liu Hong Abstract With the exploration and development of Puguang Oilfield, oil-gas exploration of carbonate rock in China has shown good prospects. Research on earthquake prediction methods for carbonate reservoir becomes the key of oil and gas exploration. Starting with analysis of geological characteristics of carbonate rock, prestack AVO inversion method, prestack elastic impedance inversion and parameter calculation method and seismic attribute extraction and optimization method were studied based on the analysis of rock physics in this work. First, variation characteristic and law of carbonate rock reservoir parameters were studied based on experimental data of rock physics, log data, analysis assay data, mud logging data and seismic data, so as to lay a foundation for the further reservoir identification and description. Then, the structure, type and propagation law of seismic wave field were analyzed through seismic forward modeling of the reservoir, and contact between information from log and geology data with elastic parameters, such as compressional wave and shear wave velocity and density were established, so as to provide a standard for reservoir identification and hydrocarbon detection using seismic reflection characteristics of the research area. Starting with the general concept of inverse problem, through analysis of Zoeppritz equation, three kinds of pre-stack inversion methods were derived and analyzed in detail, the AVO 3-parameter inversion based on Bayesian theory, the prestack AVO waveform inversion method and the simultaneous inversion method, based on the statistical hypothesis of inversion parameters and observation data and the Gauss distribution assumption of noise. The three methods were validated by model data and real data. Then, the elastic wave impedance inversion method of carbonate reservoir was investigated and the method of elastic parameter extraction from elastic impedance data was put forward. Based on the analysis of conventional methods of seismic attribute extraction and optimization, the time-frequency attributes and the wavelet attributes with time and amplitude feature were presented, and the prestack seismic attribute calculation method which can characterize the reservoir rock and fluid characteristic was presented. And the optimization of seismic attribute using the nonlinear KPCA method was also put forward. A series of seismic prediction technologies for carbonate reservoir were presented based on analysis of rock physics and seismic forward simulation technology. Practical application of these technologies was implemented in A oil field of Southern China and good effect has been achieved. Key words: carbonate rock; reservoir prediction; rock physics, prestack seismic inversion; seismic attribute

甘肃临夏盆地龙担剖面晚新生代磁性地层与古环境

Linxia Basin, situated in the northeast belt of the Tibetan Plateau, is a late Cenozoic depression basin bounded by the Tibetan Plateau and the Chinese Loess Plateau. The Cenozoic deposition, spanning over 30Ma, in which very abundant mammal fossils were discovered, is very suitable for study of uplift processes and geo-morphological evolution of the Tibetan Plateau. The Longdan section (35°31′31.6″N，103°29′0.6″E) is famous for the middle Miocene Platybelodon fauna and the late Miocene Hipparion fauna for a long time and is also one of the earliest known places for wooly rhino, which lies on the east slope of Longdan, a small village of township Nalesi in the south of the Dongxiang Autonomous County, Linxia Hui Nationallity Autonomous Prefecture. The Longdan mammal fauna was discovered at the base of the Early Pleistocene loess deposits at Dongxiang, where the lithology is different from the typical Wucheng Loess on the Chinese Loess Plateau. The rich fossils contain many new species and the major two layers of fossils are in the loess beds. Geologically the fossiliferous area is located in the central part of the Linxia Cenozoic sedimentary basin. Tectonically the Linxia Basin is an intermountain fault basin, bordered by the Leijishan major fault in the south and the north Qinling and Qilianshan major faults in the north. The section is 51.6m thick above the gravel layer, including the 1.6m Late Pleistocene Malan Loess on the top and the other loess-paleosol sequences in the middle of the section. The base of the section is the Jishi Formation, consisting of gravel layer of 13 ~ 17m thick. In this study, 972 bulk samples were collected with an interval of 5cm and other 401 orientied samples were taken with a magnetic compass. In the laboratory, the paleomagnetism, medium grain size, susceptibility, color, micromorphology, anisotropy of magnetic susceptibility were analyzed. From the stratigraphic analysis, the Longdan section from the top 0.3m to the bottom 51.6m, containing 5 normal polarities (N1-N5) and 5 reversal polarities (R1-R5). The paleomagnetic results show N3 is the Olduvai subchron in the middle of the Matuyama chron, and then the chronology of the Longdan mammal fauna is constructed along the section. The Matuyama-Gauss boundary is 45m and N5 enters Gauss chron. The Olduvai subchron with the age of 1.77 ~ 1.95Ma is found just in the upper fossiliferous level of Longdan mammal fauna. Taking the deposit rate of the section into account, the geological age of the upper fossiliferous level of Longdan mammal fauna is estimated to be about 1.9Ma. The lower fossiliferous level is just below the Reunion subchron and its age is estimated to be 2.25Ma. In addition, anisotropy of magnetic susceptibility of the loess-paleosol and other climatic indexes were used for discussing the late Cenozoic paleoenvironmental changes at Longdan, from which the Longdan area should have been an area of predominantly steppe the same as the Longdan mammal fauna.

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

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.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.