积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

山东莒南地区韧性剪切带的变形特征与形成演化

Ju Nan of Shandong province is located at southwest of Sulu UHP (ultrahigh-pressure) metamorphic terrane. It is composed of gneiss, paragneiss, eclogites, ultramafic rocks, marble and quartzite. A large ductile shear zone extends east-west has been found at the Zhubian, The south of Junan county. The Zhubian ductile shear zone is composed of high srain rock and mylonites. The mylonites fall into 3 types: Initial gneiss mylonite, mylonite and altramylonit.obvious lineation of penetration,foliation,S-Cfabrics,porphyroclasts,folds,irregularundulatory,extinction,subgrain boundary, dynamic recrystallization microstructure, core-mantle structure and are common in the ductile shear zone. Based on field work and microstructural analyse, a conclution is arrived: The ductile shear zone is an approximately SE trending faults. The Zhubian ductile shear zone formed at Ep ―Hb facies conditions which could be proved by deformaed and metamorphosed mineral aggregates, Deformation behavior, Ternary-feldspar geothermometry and so on. Zircon MC―ICP―MS U-Pb analysis is performed on the mylonite and have an average age ―835.9±13.9Ma, it’s the primary rocks formed age. The Zhubian ductile shear zone maybe formed at 224-242Ma.

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

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.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

K/E之交环境地球化学元素的演化与全球变化及恐龙绝灭的关系-以黑龙江中下游地区和俄罗斯远东地区为例

In this paper, We analyzed the geological and geographical settings of dinosaurs extinction at the end of Cretaceous, especially the effect of the change of the elements contents on dinosaurs extinction. We studied basis on the two typical sections-Cretaceous-Paleocene boundary (Baishantou section (in Jiayin, Heilongjiang province of China) and Arkhara-Boguchan Coal Mine section (in Far East of Russian)) and Longgushan section (in Jiayin, Heilongjiang province of China) mainly. This work provided some evidences for forecasting the effects of global environmental change on bio-circle. The followings are the primary gains: According to the paleo-climate indexes (CaO/MgO,Sr/Ba) and the results of Factor Analysis, we found that there were similar climate in Baishantou section and Arkhara-Boguchan Coal Mine section near the K/E boundary, and both of them took on the trend of temperature declining and precipitation heightening after transitory high-temperature and drought. There are similar change and evlution rule of the elements contents near the boundary in the both sections (Baishantou section and Arkhara-Boguchan Coal Mine section). Both iron group elements and chalcophile elements appeared obvious abnormity. There are not visible correlation between the change of elements contents and climate indexes. This shows that the elements abnormity maybe came from the factors excluding climate or the factors were too many to conceal the influence of climate. --The result of cluster analysis showed that the boundary between BST3-8 and BST3-9 may be the K/E boundary of Baishantou section, and the top of twofold coal were the K/E boundary of Arkhara-Boguchan Coal Mine section which was consistent with accepted conclusion formerly. By contrast of elements contents in dinosaur bones and general organism, in surrounding rock and general sand stone, the regulation of the change of elements contents in dinosaur bones and surrounding rock, we confirmed that dinosaur extinction in Jiayin were relative with the high abnormities of Sr, Ba, Pb, Cr and the low abnormity of Zn, at least, it was them which speeded up dinosaurs extinction. After a series of analysis, we concluded that dinosaurs extinction of this areas tied up with the relative high background values of geo-chemical elements , paleo-climate and disaster incidents. First of all, high background values provided the necessary condition for the accumulation of the elements. Secondly, the drought climate adverse to the survival of dinosaurs, and led them to extinct gradually. finally, disaster incidents, the eruption of volcano or the collision of aerolites, made them exit this planet.

元谋地区晚新生代河湖相沉积物磁性地层学及古环境研究

Yuanmou area lies on the southeastern edge of the Tibetan Plateau, the middlesegment of Yunnan-Sichuan North-South Extending Tectonic Belt and the upper reaches of the Yangztze River, which is renowned for its thick late Cenozoic fluvial-lacustrine sequences that yield rich mammalian fossils including hominoid and early human. The sediments provides great potentials for understanding the relationships between uplift of the Tibetan Plateau, evolution of hominoids and other mammalian and evolution and formation of basins in Southwest China since late Miocene. However, lithostratigrphic and chronologic views on them are controversial and hinder further discussion of the relationships of them. To this end, we selected the Baozidongqing section and the Dapoqing section to carry out systematic lithostratigraphic, magnetostratlgraphic and environmental magnetism researches in this area.The Baozidongqing section was dated to about 10.95-7.17 Ma. The age estimation of the topmost hominoid-bearing layer was about 7.43-7.17 Ma. Rock magnetic results show that the dominant magnetic carrier is hematite, with minor amount of magnetite. Both the composition and concentration of magnetic minerals strongly correlate with the lithostratigraphy, indicating that Yuanmou basin is characterized by alternating of long-term torrid-humid climate and short-term dry-hot climate. But the pattern of these short-term hot-dry events, including both the lasting time and the frequency of their occurrence dramatically changed since -8.1 Ma. Our results infer that the drying process of the Asian west interior and a significant uplift of the Tibetan Plateau would have probably caused jointly the extinction of hominoids, or the emigration of hominoids from Yuanmou to adjacent relatively torrid-humid areas.The strata between the upper of the Dapoqing section, the Niujianbao Hill and Shangnabang area can be linked by three mark layers of conglomerate, which is rather continuous and coherent than physical disturbance by new tectonic activities. Rock magnetic studies indicate that hematite is the main magnetic carriers. The section is dated back to about 2.8-1.37 Ma. Its paleocurrent flowed northeastward, which was a close and stagnant river and swamp environment about 2.2 Ma ago. Then it ran northwestward and turned into an open overflown and braid river sedimentary face during 2.2 to 1.57 Ma. Since 1.57 Ma, the paleocurrent flowed intensely northwestern and about 1.37 Ma ago, it ended the basically continuous fluvial-lacustrine deposition.

下扬子北缘前陆盆地形成演化与油气赋存条件

The foreland basin on the northern margin of the lower reach of the Yangtze river (the lower Yangtze foreland basin) is tectonically situated in the basin-mountain transitional area along the southeastern flank of the Dabie mountains. The early formation and development of the basin is closely related to the open-up of the Mian-Lue paleo-oceanic basin on the southern margin of the Central Orogenic System represented by Qinling-Dabei orogenic belt, while the tectonic evolution of the middle-late stage of the basin is mainly related to development of the Mian-Lue tectonic zone that occurred on the basis of the previous Mian-Lue paleo-suture. The foreland basin of the northern rim of the lower reach of the Yangtze river was formed during the middle-Triassic collision between the Yangtze and North China plates and experienced an evolution of occuirence-development-extinction characterized by marine facies to continental facies and continental margin to intracontinent in terms of tectonic setting.The foreland basin (T2-J2) was developed on the basis of the passive continental marginal basin on the south side of the Mian-Lue paleo-ocean and superimposed by late Jurassic-Tertiary fault basin. The tectonic setting underwent a multiple transformation of rifting-collisional clososing-tensional faulting and depression, which resulted in changes of the property for the basin and the final formation of the superposed compose basin in a fashion of 3-story-building. According to the tectonic position and evolution stages of plate collision happening on the southeastern margin of the Dabie mountains, and tectono-tratigraphic features shown by the foreland basin in its main formational period, the evolution of the foreland basin can be divided into four stages: 1) pre-orogenic passive margin (P2-Ti). As the Mian-Lue ocean commenced subduction in the late-Permian, the approaching of the Yangtze and North China plates to each other led to long-periodical and large-scale marine regression in early Triassic which was 22 Ma earlier than the global one and generated I-type mixed strata of the clastic rocks and carbonate, and I-type carbonate platform. These represent the passive stratigraphy formed before formation of the foreland basin. 2) Foreland basin on continental margin during main orogenic episode (T2.3). The stage includes the sub-stage of marine foreland basin (T2X remain basin), which formed I-type stratigrphy of carbonate tidal flat-lagoon, the sub-stage of marine-continental transition-molasse showing II-type stratigraphy of marine-continental facies lake - continental facies lake. 3) Intracontinental foreland basin during intracontinental orogeny (Ji-2)- It is characterized by continental facies coal-bearing molasses. 4) Tensional fault and depression during post-orogeny (J3-E). It formed tectono-stratigraphy post formation of the foreland basin, marking the end of the foreland evolution. Fold-thrust deformation of the lower Yangtze foreland basin mainly happened in late middle-Jurassic, forming ramp structures along the Yangtze river that display thrusting, with deformation strength weakening toward the river from both the Dabie mountains and the Jiangnan rise. This exhibits as three zones in a pattern of thick-skinned structure involved the basement of the orogenic belt to decollement thin-skinned structure of fold-thrust from north to south: thrust zone of foreland basin on northern rim of the lower reach of the Yangtze river, foreland basin zone and Jiannan compose uplift zone. Due to the superposed tensional deformation on the earlier compressional deformation, the structural geometric stratification has occurred vertically: the upper part exhibits late tensional deformation, the middle portion is characterized by ramp fault -fold deformation on the base of the Silurian decollement and weak deformation in the lower portion consisting of Silurian and Neo-Proterozoic separated by the two decollements. These portions constitutes a three-layered structural assemblage in a 3-D geometric model.From the succession of the lower reach of the Yangtze river and combined with characteristics of hydrocarbon-bearing rocks and oil-gas system, it can be seen that the succession of the continental facies foreland basin overlies the marine facies stratigraphy on the passive continental margin, which formed upper continental facies and lower marine facies hydrocarbon-bearing rock system and oil-gas forming system possessing the basic conditions for oil-gas occurrence. Among the conditions, the key for oil-gas accumulation is development and preservation of the marine hydrocarbon-bearing rocks underlying the foreland basin. The synthetic study that in the lower Yangtze foreland basin (including the Wangjiang-Qianshan basin), the generation-reservoir-cover association with the Permian marine facies hydrocarbon-bearing rocks as the critical portion can be a prospective oil-gas accumulation.Therefore, it should aim at the upper Paleozoic marine hydrocarbon-bearing rock system and oil-gas forming system in oil-gas evaluation and exploration. Also, fining excellent reservoir phase and well-preserved oil-gas accumulation units is extremely important for a breakthrough in oil-gas exploration.

条件反射性免疫抑制及其脑机制的c-fos研究

To explore the neural mechanisms underlying conditioned immunomodulation, this study employed the classical taste aversion (CTA) behavioral paradigm to establish the conditioned humoral and cellular immunosuppression (CIS) in Wistar rats, by paring saccharin (CS) with intraperitoneal (i.p.) injection of an immunosuppressive drug cyclophophamide (UCS). C-fos immunohistochemistry method was used to observe the changes of the neuronal activities in the rat brain during the acquisition, expression and extinction of the conditioned immunosuppression (CIS). The followings are the main results: 1. Five days after one trial of CS-UCS paring, reexposure to CS alone significantly decreased the level of the anti-ovalbumin (OVA) IgG in the peripheral serum. Two trials of CS-UCS paring and three reexposures to CS not only resulted in further suppression of the primary immune response, but also reduced the numbers of peripheral lymphocytes and white blood cells. This finding indicates that CS can induce suppression of the immune function, and the magnitude of the effects is dependent on the intensity of training. 2. On day 5 following two trials of CS-UCS pairing, CS suppressed the spleen lymphocytes responsiveness to mitogens ConA, PHA and PWM, and decreased the numbers of peripheral lymphocytes and white blood cells. On day 15, only PHA induced lymphocyte proliferation was suppressed by CS. On day 30, presentation of CS did not have any effect on these immune parameters. These results suggest that the conditioned suppression of the cellular immune function can retain 5-15 days, and extinct after 30 days. 3. CTA was easily induced by one or two CS-UCS parings, and remained robust even after 30 days. These data demonstrate that CIS can be dissociated from CTA, and they may be mediated by different neural mechanisms. 4. Immunohistochemistry assays revealed a broad pattern of c-fos expression throughout the rat brain following the CS-UCS pairing and reexposure to CS, suggesting that many brain regions are involved in CIS. Some brain areas including the solitary tract nucleus (Sol), lateral parabrachial nucleus (LPB) and insular cortex (IC), showed high level c-fos expressions in response to both CS and UCS, suggesting that they may be involved in the transmission and integration of the CS and UCS signals in the brain. There were dense c-FOS positive neurons in the paraverntricular nucleus (PVN) and supraoptic nucleus (SO) of hypothalamus, subfornical organ (SFO) and area postrema (AP) etc. after two trials of CS-UCS paring and after the reexposure to CS 5 days later, but not in the first training and after the extinction of CIS (30 days later). The results reflect that these nuclei may have an important role in CIS expression, and may also response to the immunosuppression of UCS. The conditioned training and reexposure to CS 5 days later induced high level c-fos expression in the cingulate cortex (Cg), central amygdaloid nucleus (Ce), intermediate part of lateral septal nucleus (LSI) and ventrolateral parabrachial nucleus (VLPB) etc. But c-fos induction was not apparent when presenting CS 30 days later. These brain regions are mainly involved in CIS, and may be critical structures in the acquisition and expression of CIS. Some brain regions, including the frontal cortex (Fr), ventral orbital cortex (VO), IC, perirhinal cortex (PRh), LPB and the medial part of solitary nucleus (SolM), showed robust c-FOS expression following the conditioning training and reexposure to CS both on day 5 and day 30, suggesting that they are critically involved in CTA.