地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

非均匀性对岩石破裂过程影响的研究

Rock heterogeneity plays an important role in rock fracturing processes. However, because fracturing is a dynamic process and it is very difficult to quantify materials' heterogeneity, most of the theories dealing with local failure were based on the homogeneity assumption, very few involving stress distribution heterogeneity and successive local failure due to rock heterogeneity. Therefore, based on various references, the author studied the laws and mechanism of influences of heterogeneity on rock fracturing processes, under the frame of the project "Study on Associate Mechanism between Rock Mass Fracture and Strength Failure", funded by Nation Natural Science Fund. the research consists of such aspects as size effect correction to rock fracture parameters, SEM (Scanning Electron Microscope) real-time observation on rock samples under different loads, micro-hardness testing, and numerical simulating based on microstructure. There are some important research results as followed: 1. Unifying formula for nonlinear and non-singularity correction, simplifying the complex process of correcting size effect on rock fracture toughness. 2. Using the methods of micro-hardness testing mineral grain and random jointing micrograph digitizing mineral slice, preliminarily solving the problems of numerical simulating and quantitatively describing the heterogeneous strength and its distribution rules, which has certain innovation and better practicability. 3. Based on SEM real-time observation, studying the micro-process of fracturing in marble, sandstone, granite, and mushroom stone samples with premanufactured cracks under tension, pure-shear and compression-shear conditions. Strength Failure was observed: there was some kind failure occurred before Fracture Failure in marble and sandstone samples with double cracks under pure-shearing. It is believed that the reason of strength failure developing is that stress concentrations is some locations are larger than that near the end of pre-manufactured cracks. 4. Based on the idea that rock macro-constitute is composed of complex microstructure, the promising method used to handle heterogeneity considers not only the heterogeneity of the rock medium, but also the heterogeneity of the rock structure. 5. Putting forward two types of rock strength failure: medium strength failure induced by heterogeneity of rock medium and structure strength failure induced by heterogeneity rock structure. 6. By evaluating potential fracture cell with proper failure priority, the numerical simulating method solved the problem of simulating the coextensive strength failure and fracture failure with convention strength failure rules. The result of numerical analysis shows that the influence of heterogeneity on rock fracturing processes is evident. The sinuosity of the rock fracture-propagation path, and the irregular fluctuation of loading displacement curve, is mainly controlled by the heterogeneity of rock medium.

河流成因储层沉积微相与剩余油分布及形成机理研究

Based on the study of fluvial sandstone reservoir in upper of Guantao group in Gudao and Gudong oilfields, this paper first introduces A.D.Miall's(1996a) architectural-element analysis method that was summarized from ground outcrop scale into the reservoir formation research of the study area, more subtly divides sedimentary microfacies and establishes sedimentary model of research area.on this base, this paper summarizes the laws of residual oil distribution of fluvial formation and the control effect of sedimentary microfacies to residual oil distribution, and reveals residual oil formation mechanism. These results have been applied to residual oil production, and the economic effect is good. This paper will be useful for residual oil research and production and enhancement of oil recovery in similar reservoir. The major conclusions of this paper are as follows. 1. Using the architectural-element analysis method to the core data, a interfacial division scheme of the first to the dixth scale is established for the studied fluvial formation. 2.Seven architectural-elements are divided in upper of Guantao group of study area. The sandstone group 5～1+2 of Neogene upper Gutao group belongs to high sinuous fine grain meandering river, and the sandstone group 6 is sandy braided river. 3. Inter layer, the residual oil saturation of "non-main layer" is higher than "main layer", but the residual recoverable reserve of former is larger. Therefore, "main layer" is the main body of residual oil distribution. The upper and middle part of inner layer has lower permeability and strong seeping resistance. Addition to gravity effect in process of driving, its driving efficiency is low; residual oil saturation is high. Because of controlling of inside non-permeable interlayer or sedimentary construction, the residual oil saturation of non-driving or lower driving efficiency position also is high. On plane, the position of high residual oil saturation mostly is at element LV, CS, CH (FF), FF etc, Which has lower porosity and permeability, as well as lens sand-body and sand-body edge that is not controlled by well-net, non-perfect area of injection and production, lower press difference resort area of inter-well diffiuent-line and shelter from fault, local high position of small structure. 4.Microscopic residual oil mainly includes the non-moved oil in the structure of fine pore network, oil in fine pore and path, oil segment in pore and path vertical to flow direction, oil spot or oil film in big pore, residual oil in non-connective pore. 5.The most essential and internal controlling factor of fluvial formation residual oil distribution is sedimentary microfacies. Status of injection and production is the exterior controlling factor of residual oil distribution. 6. The controlling effect of formation sedimentary microfacies to residual oil distribution indicates inter-layer vertical sedimentary facies change in scale of injection and production layer-series, planar sedimentary face change and inner-layer vertical sedimentary rhythm and interbed in single layer to residual oil distribution. 7. It is difficult to clear up the inter-layer difference in scale of injection and production layer-series. The using status of minor layer is not good and its residual oil saturation is high relatively. It is obvious that inter-layer vertical sedimentary facies changes control inter-layer residual oil distribution at the same or similar conditions of injection and production. For fluvial formation, this vertical sedimentary facies change mainly is positive gyration. Namely, from down to top, channel sediment (element CHL, LA) changes into over-bank sediment (element LV, CR, CS). 8. In water-injection developing process of transverse connecting fluvial sandstone oil formation, injection water always comes into channel nearby, and breaks through along channel and orientation of high pressure gradient, does not expand into side of channel until pressure gradient of channel orientation changes into low. It brings about that water-driving status of over-bank sedimentary element formation (LV, CR, CS) is not good, residual oil saturation is high. In non-connective abandoned channel element (CH) formation with channel, because this reverse is difficult to control by injection and production well-series, its using status is not good, even terribly not good, residual oil is enrichment. 9. The rhythm and sedimentary structure, sedimentary facies change in single sand body brings about vertical changes of formation character, growth character of inner layer interbed. These are important factor of controlling and affecting vertical water spread volume and inner layer residual oil forming and distribution in single sand body. Positive rhythm, is the principal part of fluvial sandstone inner layer sedimentary rhythm. Namely, from down to upside, rock grain granularity changes from coarse to fine, seeping ability changes from strong to feebleness. It brings about that water-driving status of inner layer upside is not good, residual oil saturation is high. Inner layer interbed has different degree affecting and controlling effect to seeping of oil and water. Its affecting degree lies on interbed thickness, extending scale, position, and jeted segment of production or injection well. The effect of interbed at upside of oil formation to oil and water seeping is less; the effect of interbed at middle of oil formation to oil and water seeping is more. 10. Indoor experiment and research indicate that wettability, permeability step, vertical permeability, position of Kmax and ratio of oil viscousity and water viscousity all have great effect on the water-driving recovery ratio. 11. Microscopic residual oil distribution is affected and controlled by formation pore network structure, pressure field distribution, and oil characteristic. 12.The residual oil forming mechanism: the over-bank sedimentary element and upper part of a positive rhythm sandstone have fine pore and throat network, permeability is low, displacement pressure of pore and throat is high. The water-driving power usually falls short of displacement pressure that brings about injection water does not spread into these pore and throat network, thereby immovable oil area, namely residual oil, is formed. At underside of channel sedimentary element and positive rhythm sandstone, porosity and permeability is relatively high, connecting degree of pore and throat is high, displacement pressure of pore and throat is low. Thereby injection water is easy to enter into pore and throat, driving oil in them. Because the pore space is irregular, the surface of pore wall is coarse and non-flat. That the oil locate on concave hole of pore wall and the dead angle of pore, and the oil attaches on surface of pore wall by surface tension, are difficult to be peeled off, becoming water-driving residual oil (remaining oil). On the other hand, Because flowing section lessens, flowing resistance increase, action of capillary fore, or seeping speed decreases at process of transfer at pass narrow throat path in the course carried by driving water. The "oil drop", "oil bead", or "oil segment" peeled off by driving water is difficult to carry and to drive out by water at less pressure difference. Thereby they are enclosed in pore to form discontinuous residual oil. 13.This results described above have been applied in nine develop blocks of Gudao and Gudong oilfield. Its applying effect is marked through local injection production adjustment, deploying replacement well, repair hole, replacement envelop, block off water and profile control etc. Relative method and technology can be applied to other oil production area of Shengli oilfield, and obtain better economic and societal effect.

沾化凹陷泥质岩油气藏成藏特征研究

Mudstone reservoir is a subtle reservoir with extremely inhomogeneous, whose formation is greatly related to the existence of fracture. For this kind of reservoir, mudstone is oil source rock, cover rock and reservoir strata, reservoir type is various, attitude of oil layer changes greatly, and the distribution of oil and gas is different from igneous or clastic rock reservoir as well as from carbonate reservoir of self-producing and self-containing of oil and gas. No mature experience has been obtained in the description, exploration and development of the reservoir by far. Taking Zhanhua depression as an example, we studied in this thesis the tectonic evolution, deposit characteristics, diagenesis, hydrocarbon formation, abnormal formation pressure, forming of fissure in mudstone reservoir, etc. on the basis of core analysis, physical simulation, numerical simulation, integrated study of well logging and geophysical data, and systematically analyzed the developing and distributing of mudstone fissure reservoir and set up a geological model for the formation of mudstone fissure reservoir, and predicted possible fractural zone in studied area. Mudstone reservoir mainly distributed on the thrown side of sedimentary fault along the sloping area of the petroleum generatiion depression in Zhanhua depression. Growing fault controlled subsidence and sedimentation. Both the rate of subsidence and thickness of mudstone are great on the thrown side of growing fault, which result in the formation of surpressure in the area. The unlocking of fault which leads to the pressure discharges and the upward conduct of below stratum, also makes for the surpressure in mudstone. In Zhanhua depression, mudstone reservior mainly developed in sub-compacted stratum in the third segment of Shahejie formation, which is the best oil source rock because of its wide spread in distribution, great in thickness, and rich in organic matter, and rock types of which are oil source mudstone and shale of deep water or semi-deep water sediment in lacustrine facies. It revealed from core analysis that the stratum is rich in limestone, and consists of lamina of dark mudstone and that of light grey limestone alternately, such rock assemblage is in favor of high pressure and fracture in the process of hydrocarbon generation. Fracture of mudstone in the third segment of Shahejie formation was divided into structure fracture, hydrocarbon generation fracture and compound fracture and six secondary types of fracture for the fist time according to the cause of their formation in the thesis. Structural fracture is formed by tectonic movement such as fold or fault, which develops mainly near the faults, especially in the protrude area and the edge of faults, such fracture has obvious directivity, and tend to have more width and extension in length and obvious direction, and was developed periodically, discontinuously in time and successively as the result of multi-tectonic movement in studied area. Hydrocarbon generation fracture was formed in the process of hydrocarbon generation, the fracture is numerous in number and extensively in distribution, but the scale of it is always small and belongs to microfracture. The compound fracture is the result of both tectonic movement and hydrocarbon forming process. The combination of above fractures in time and space forms the three dimension reservoir space network of mudstone, which satellites with abnormal pressure zone in plane distribution and relates to sedimentary faces, rock combination, organic content, structural evolution, and high pressure, etc.. In Zhanhua depression, the mudstone of third segment in shahejie formation corresponds with a set of seismic reflection with better continuous. When mudstone containing oil and gas of abnormal high pressure, the seismic waveform would change as a result of absorb of oil and gas to the high-frequency composition of seismic reflection, and decrease of seismic reflection frequency resulted from the breakage of mudstone structure. The author solved the problem of mudstone reservoir predicting to some degree through the use of coherent data analysis in Zhanhua depression. Numerical modeling of basin has been used to simulate the ancient liquid pressure field in Zhanhua depression, to quantitative analysis the main controlling factor (such as uncompaction, tectonic movement, hydrocarbon generation) to surpressure in mudstone. Combined with factual geologic information and references, we analyzed the characteristic of basin evolution and factors influence the pressure field, and employed numerical modeling of liquid pressure evolution in 1-D and 2-D section, modeled and analyzed the forming and evolution of pressure in plane for main position in different periods, and made a conclusion that the main factors for surpressure in studied area are tectonic movement, uncompaction and hydrocarbon generation process. In Zhanhua depression, the valid fracture zone in mudstone was mainly formed in the last stage of Dongying movement, the mudstone in the third segment of Shahejie formation turn into fastigium for oil generation and migration in Guantao stage, and oil and gas were preserved since the end of the stage. Tectonic movement was weak after oil and gas to be preserved, and such made for the preserve of oil and gas. The forming of fractured mudstone reservoir can be divided into four different stages, i.e. deposition of muddy oil source rock, draining off water by compacting to producing hydrocarbon, forming of valid fracture and collecting of oil, forming of fracture reservoir. Combined with other regional geologic information, we predicted four prior mudstone fracture reservoirs, which measured 18km2 in area and 1200 X 104t in geological reserves.

地震反射系数成像方法研究

Reflectivity sequences extraction is a key part of impedance inversion in seismic exploration. Although many valid inversion methods exist, with crosswell seismic data, the frequency brand of seismic data can not be broadened to satisfy the practical need. It is an urgent problem to be solved. Pre-stack depth migration which developed in these years becomes more and more robust in the exploration. It is a powerful technology of imaging to the geological object with complex structure and its final result is reflectivity imaging. Based on the reflectivity imaging of crosswell seismic data and wave equation, this paper completed such works as follows: Completes the workflow of blind deconvolution, Cauchy criteria is used to regulate the inversion(sparse inversion). Also the precondition conjugate gradient(PCG) based on Krylov subspace is combined with to decrease the computation, improves the speed, and the transition matrix is not necessary anymore be positive and symmetric. This method is used to the high frequency recovery of crosswell seismic section and the result is satisfactory. Application of rotation transform and viterbi algorithm in the preprocess of equation prestack depth migration. In equation prestack depth migration, the grid of seismic dataset is required to be regular. Due to the influence of complex terrain and fold, the acquisition geometry sometimes becomes irregular. At the same time, to avoid the aliasing produced by the sparse sample along the on-line, interpolation should be done between tracks. In this paper, I use the rotation transform to make on-line run parallel with the coordinate, and also use the viterbi algorithm to complete the automatic picking of events, the result is satisfactory. 1. Imaging is a key part of pre-stack depth migration besides extrapolation. Imaging condition can influence the final result of reflectivity sequences imaging greatly however accurate the extrapolation operator is. The author does migration of Marmousi under different imaging conditions. And analyzes these methods according to the results. The results of computation show that imaging condition which stabilize source wave field and the least-squares estimation imaging condition in this paper are better than the conventional correlation imaging condition. The traditional pattern of "distributed computing and mass decision" is wisely adopted in the field of seismic data processing and becoming an obstacle of the promoting of the enterprise management level. Thus at the end of this paper, a systemic solution scheme, which employs the mode of "distributed computing - centralized storage - instant release", is brought forward, based on the combination of C/S and B/S release models. The architecture of the solution, the corresponding web technology and the client software are introduced. The application shows that the validity of this scheme.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

中新世以来中国北方风尘堆积的粘土矿物组成及其古环境意义

In northern China, the loess-soil sequence of the last 2.6 Ma, the Hipparion Red-Earth of eolian origin and recently reported Pliocene-Miocene loess-soil sequence provide a near continuous continental eolian record of climatic history for the past 22.0 Ma. This work aims to investigate the composition and structure of clay minerals contained in deposits, and to explore their implications for environmental evolutions over the last 22.0 Ma. Clay minerals, which were extracted from eolian samples collected at Xifeng (0-6.2 MaBP) and Qinan (6.2-22.0 MaBP) sections, were analyzed qualitatively and semi-quantitatively by using X-ray diffraction (XRD), differential thermal analysis (DTA), thermogravimetric analysis (TG) and chemical analysis. The main conclusions are as follows: Over the last 22.0 Ma, the clay mineral assemblage among Quaternary loess-soils, Hipparion Red-Earth, and Miocene loess-soils shows similar components, mainly consisting of illite (55-80%), kaolinite (7-20%), chlorite (4-13%), smectite (2-23%) as results calculated by comparing major peak areas. There are no obvious differences in both types and amounts of clay minerals between loess and interbedded soils, suggesting that overwhelming part of the clay minerals is derived from the source. According to the components of clay minerals, the whole sequence of eolian deposits in the Loess Plateau can be divided into ten clay mineral assemblage zones over the last 22.0 Ma, whose corresponding ages are: 22.0-21.0 MaBP, 21.0-18.0 MaBP, 18-16.2 MaBP, 16.2-13.0 MaBP, 13.0-10.0 MaBP, 10.0-5.5 MaBP, 5.5-4.4 MaBP, 4.4-2.8 MaBP, 2.8-1.0 MaBP, 1.0-0 MaBP, respectively. This may imply that dust supply changed at least nine times over the past 22.0 Ma. The loess illite has a better crystaliinity, higher value of the FWHM and IC, than the interbedd soils. Previous studies indicated that irregular mixed layer minerals could form under relatively warm and humid conditions (Han, 1982). According to the general distribution of clay minerals of zonal soil (Chamley, 1989), the clay mineral assemblage of eolian deposits in Xifeng and Qinan sections is typical of temperature-humid and warm-subarid environment. Therefore, our results indicate climatic environment in Loess Plateau did not change remarkably since 22.0 Ma, and fluctuated between temperature-humid and warm-subarid climate. 4. The illite generally presents poorer crystaliinity during the period of 22.0 to 2.8 MaBP than in the last 2.8 Ma BP, especially at the intervals of 3.5-4.5 Ma BP, 14.0-17.0 MaBP and 20.0-22.0 Ma BP, which indicates that the weathering intensity was stronger in Neogene than in Quaternary. 5. The relatively low ice volume and high global temperature may be responsible for the strange weathering intensity during the interval of the 3.5-4.5 Ma BP, 14.0-17.0 Ma BP and 20.0-22.0 Ma BP.

东南极格罗夫山土壤研究与古环境分析

Three soil spots were found in Grove Mountains, east Antarctica during 1999-2000, when the Chinare 16th Antarctic expedition teams entered the inland Antarctica. The characteristics of soils in Grove Mountains are desert pavement coating the surface, abundant water soluble salt, negligible organ matter, and severe rubification and salinization, scarces of liquid water, partly with dry permafrost, corresponding with the soils of McMurdo, Transantarctic. The soils age of Grove Mountains is 0.5-3.5Ma. Podzolization and redoximorphism are the main features in coastal Wilks region, in addition, there is strong enrichment of organic matter in many soils of this region. The main soil processes of Fildes Peninsula of King George Island include the intense physical weathering, decalcification and weakly biochemical processes. Peat accumulation is the main processes in Arctic because of humid and cold environment.Based on synthesis of heavy minerals, particle size, quartz grain surface textures, as well as pollen in soils, the soils parent materials of Grove Mountains derived from alluvial sediment of the weathering bedrocks around soils, and formed during the warm period of Pliocene. The detailed information is followed .l)The results of heavy minerals particle size showed the parent minerals derived form the weathering bedrocks around soils. 2)The quartz sand surface textures include glacial crushing and abrasion such as abrasive conchoidal fractures and grain edges, abrasive subparallel linear fractures and angularity, subaqueous environments produce V-shaped and irregular impact pits, polished surface, and chemical textures, such as beehive solution pits, which showed the water is the main force during the sediment of the soil parent minerals. 3)The pollen consist of 40 plant species, of which at least 5 species including Ranunculaceae, Chenopodiaceae, Artemisia, Gramineae, Podocarpus belong to the Neogene vegetation except the species from the old continent. Compared with Neogene vegetation of Transantarctic Mountains, Antarctic, we concluded that they grow in warm Pliocene.

大方不规则I/S间层的铝层柱化研究

层柱粘土是一类新型的矿物功能材料,具有大孔径、大比表面积、孔径可调、热稳定性好、表面酸性强及含有金属活性组分等特点.它是当前国际粘土矿物学、化学及材料科学界研究的一个热点.该文以贵州大方猫场不规则伊/蒙（简称I/S）间层为研究对象,利用聚合羟基铝作为柱化剂对其进行了层柱化研究,系统地介绍了大方不规则I/S间层的矿物学特征、柱化剂Al<,13>的形态特征、铝层柱I/S间层的制备及表征等.

伊/蒙不规则间层粘土的矿物学特征及其层柱化研究

该文首先详细而系统地介绍了一类新型纳米矿物材料——层柱粘土.论文以贵州大方I/S不规则有序间层粘土为例,在采用多种测试手段了解该矿物的组成、结构及性能的基础上,首次对I/S不规则间层粘土进行了较系统的铝、锆、钛的层柱化研究.XRD、HRTEM、FTIR及DTA-TG共同揭示了大方I/S间层粘土是以Rl型为主的I/S不规则有序间层粘土,化学分析数据计算表明结构中的八面体为二八面体,蒙皂石为钠型蒙脱石.通过对大方I/S不规则间层粘土的铝层柱研究发现,较低的初始铝浓度有利于形成较多的Keggin离子,产物层间距可达3.2～3.5nm,N<,2-BET比表面积达110m<'2>/g.首次对I/S不规则间层粘土进行了钻、钛层柱化研究,所获产物层间距都大于铝层柱者,分别为3.3～4.0nm和3.5~4.7nm.文章首次对层柱I/S不规则间层粘土进行了除层间距以外的孔结构、表面酸性、(水)热稳定性和显微形貌的表征.探讨了层与柱之间、颗粒与颗粒之间的排列和连接对层柱产物性能的影响,比较成功地用IR证据解释了层在机理中四面体的反转问题.最后初步试验了层柱I/S不规则间层粘土对水中重金属pb<'2+>和垃圾淋滤液中COD的去除效果,结果表明卫Ti-I/S对Pb<'2+>的吸附容量可达8.39mg/g,对COD的去除率可高达76.71%.

运动与重定向时物-我、物-物空间表征的使用

Eight experiments tested how object array structure and learning location influenced the establishing and utilization of self-to-object and object-to-object spatial representations in locomotion and reorientation. In Experiment 1 to 4, participants learned either at the periphery of or amidst regular or irregular object array, and then pointed to objects while blindfolded in three conditions: before turning (baseline), after rotating 240 degrees (updating), and after disorientation (disorientation). In Experiment 5 to 8, participants received instruction to keep track of self-to-object or object-to-object spatial representations before rotation. In each condition, the configuration error, which means the standard deviation of the means per target object of the signed pointing errors, was calculated as the index of the fidelity of representation used in each condition. Results indicate that participants form both self-to-object and object-to-object spatial representations after learning an object-array. Object-array structure influences the selection of representation during updating. By default, object-to-object spatial representation is updated when people learned the regular object-array structure, and self-to-object spatial representation is updated when people learned the irregular object array. But people could also update the other representation when they are required to do so. The fidelity of representations will confine this kind of “switch”. People could only “switch” from a low fidelity representation to a high fidelity representation or between two representations of similar fidelity. They couldn’t “switch” from a high fidelity representation to a low fidelity representation. Leaning location might influence the fidelity of representations. When people learned at the periphery of object array, they could acquire both self-to-object and object-to-object spatial representations of high fidelity. But when people learned amidst the object array, they could only acquire self-to-object spatial representation of high fidelity, and the fidelity of object-to-object spatial representation was low.

汉字阅读中的字形、字音、字义及频率效应的研究

Unlike alphabetic languages, Chinese language is ideographical writing system. Each Chinese character is single-syllable and usually has a direct meaning. So Chinese characters are a kind of valuable experimental material used for research on reading and comparisons of the reading mechanism of different language. In this paper, the normal persons and the patients with semantic dementia were respectively scheduled for two parts of experimental studies on the orthographic, phonologic, semantic and frequency effects of reading of Chinese characters. The Stroop-like character-picture interference experimental paradigm was used to investigate the orthographic, phonologic, semantic and frequency effects of Chinese characters on picture naming when they were presented with pictures to normal persons. The results indicated that the orthographic facilitation effect, phonologic facilitation effect, and semantic interference effect occurred at different SOA values. The orthographic and phonologic facilitation effects were independent. It was for the first time shown that the interaction between orthographic variable and semantic variable occurred when the high-frequency Chinese characters were read. Phonologic representation was activated quicker than semantic representation, by comparison of their SOA. Generally, it means that there is reading without meaning in Chinese character among the normal persons. The orthographic, phonologic, semantic, frequency and concrete effects of Chinese characters were further investigated among the dementia patients with DAT(dementia of Alzheimer's type disease) or CVA or both. They all have an impaired semantic memory. The results showed that patients with dementia could read the names of the pictures aloud while they could not name them or match them with a right character correctly. This is reading impairment without meaning in Chinese among the dementia patients. Meanwhile, they had a selective reading impairment and more LARC(a legitimate alternative reading of components) mistakes especially when reading low-frequency irregular, low-frequency inconsistent and abstract Chinese characters. With the patients' semantic impairment developed, their ability to read the pictures names would remain whereas their ability to read low-frequency irregular and low-frequency inconsistency Chinese characters was reduced. These results indicated that low-frequency irregular Chinese characters can be read correctly only when it is supported by their semantic information. Based on the above results of reading without meaning and of reading of low-frequency irregular Chinese characters supported by their semantic information, it is reasonable to suggest that at least two routes are involved in the process of reading Chinese characters. They are direct phonologic route and indirect semantic route; moreover, the two routes are independent.