Tectonic and geodynamic pattern of marginal seas on the West Pacific inferred from multi-satellite altewetry geoid

With the rapid development of satellite observations, we can use the altimetry geoid to study submarine tectonics and geodynamics. On the basis of the 4' x 4' geoid undulation calculated from altimeter data of Geosat, ERS-1/2 and Topex/Poseidon on the West Pacific, located between 0degreesN similar to 45degreesN, 100degreesC similar to 150degreesE, Bouguer, Glenni and isostatic geoid undulation are obtained from correction of gravitational potential of the global topography and isostacy. Moho discontinuity depth is inversed by the Glenni geoid undulation, and the stress field from small-scale mantle convection is reasonably calculated from the isostatic geoid undulation. The results show that within the Philippine Sea and the South China Sea, short-wavelength lineations of the geoid undulation are parallel or cross to magnetic lineations and rifting ridges. The Moho depth of marginal sea basins becomes shallow southward, and its values are similar to that of the Philippine Sea. These facts show that strength of tectonic activities are almost the same on the both sides of the Ryukyu-Taiwan-Philippine are. Various kinds of tectonic features with different driving mechanisms of small-middle and large-scale of mantle convection, however, display a special pattern of tectonics and geodynamics of the continental marginal seas distinguished from oceans and continents.

Periodical mixing of MORB magmas near East Pacific Rise 13N: evidence from modeling and zoned plagioclase phenocrysts

Thirty-six basalt samples from near East Pacific Rise 13N are analyzed for major and trace elements. Different types of zoned plagioclase phenocrysts in basalts are also backscatter imaged, and major element profiles scanned and analyzed for microprobe. Basalts dredged from a restricted area have evolved to different extents (MgO=9.38wt%-6.76wt%). High MgO basalts are modeled for crystallization to MgO of about 7wt%, and resulted in the Ni contents (28 ppm) that are generally lower than that in observed basalts (> 60 ppm). It suggests that low MgO basalts may have experienced more intensive magma mixing. High MgO (9.38wt%) basalt is modeled for self-"mixing-crystallization", and the high Ni contents in low MgO basalts can be generated in small scale and periodical self-mixing of new magma (high MgO). "Mixing-crystallization" processes that low MgO magmas experienced accord with recent 226Ra/230Th disequilibria studies for magma residence time, in which low MgO magmas have experienced more circles of "mixing-crystallization" in relatively longer residence time. Magma mixing is not homogeneous in magma chamber, however, low MgO magmas are closer to stable composition produced by periodical "mixing-crystallization", which is also an important reason for magma diversity in East Pacific Rise. Zoned plagioclase phenocrysts can be divided into two types: with and without high An# cores, both of which have multiple reversed An# zones, suggesting periodical mixing of their host magmas. Cores of zoned plagioclase in low MgO (7.45wt%) basalt differ significantly with their mantle in An#, but are similar in An# with microlite cores (products of equilibrium crystallization) in high MgO (9.38wt%) basalt, which further shows that plagioclase phenocryst cores in low MgO basalts may have formed in their parental magmas before entering into the magma chamber.

COASTAL SEA-LEVEL AND THE LARGE-SCALE CLIMATE STATE - A DOWNSCALING EXERCISE FOR THE JAPANESE ISLANDS

A major problem which is envisaged in the course of man-made climate change is sea-level rise. The global aspect of the thermal expansion of the sea water likely is reasonably well simulated by present day climate models; the variation of sea level, due to variations of the regional atmospheric forcing and of the large-scale oceanic circulation, is not adequately simulated by a global climate model because of insufficient spatial resolution. A method to infer the coastal aspects of sea level change is to use a statistical ''downscaling'' strategy: a linear statistical model is built upon a multi-year data set of local sea level data and of large-scale oceanic and/or atmospheric data such as sea-surface temperature or sea-level air-pressure. We apply this idea to sea level along the Japanese coast. The sea level is related to regional and North Pacific sea-surface temperature and sea-level air pressure. Two relevant processes are identified. One process is the local wind set-up of water due to regional low-frequency wind anomalies; the other is a planetary scale atmosphere-ocean interaction which takes place in the eastern North Pacific.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

基于Multi-Agent的生产过程建模与仿真技术研究

某些流程行业由于采用按配方进行分组加工的模式组织生产,在排产时存在多条路径调度优化的问题,应用一般的优化算法对于现场在线调度难以给出满意结果,而基于Agent的过程仿真在解决离散、非线性系统模拟方面有显著的优势,本文采用Agent的方法对生产过程建模,然后对方案组内的备选方案进行仿真,通过对比各方案的仿真结果找到最优的方案作为执行方案,为现场的优化排产提供决策支持。

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

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

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

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

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

弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

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

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.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

Automatically Recovering Geometry and Texture from Large Sets of Calibrated Images

Three-dimensional models which contain both geometry and texture have numerous applications such as urban planning, physical simulation, and virtual environments. A major focus of computer vision (and recently graphics) research is the automatic recovery of three-dimensional models from two-dimensional images. After many years of research this goal is yet to be achieved. Most practical modeling systems require substantial human input and unlike automatic systems are not scalable. This thesis presents a novel method for automatically recovering dense surface patches using large sets (1000's) of calibrated images taken from arbitrary positions within the scene. Physical instruments, such as Global Positioning System (GPS), inertial sensors, and inclinometers, are used to estimate the position and orientation of each image. Essentially, the problem is to find corresponding points in each of the images. Once a correspondence has been established, calculating its three-dimensional position is simply a matter of geometry. Long baseline images improve the accuracy. Short baseline images and the large number of images greatly simplifies the correspondence problem. The initial stage of the algorithm is completely local and scales linearly with the number of images. Subsequent stages are global in nature, exploit geometric constraints, and scale quadratically with the complexity of the underlying scene. We describe techniques for: 1) detecting and localizing surface patches; 2) refining camera calibration estimates and rejecting false positive surfels; and 3) grouping surface patches into surfaces and growing the surface along a two-dimensional manifold. We also discuss a method for producing high quality, textured three-dimensional models from these surfaces. Some of the most important characteristics of this approach are that it: 1) uses and refines noisy calibration estimates; 2) compensates for large variations in illumination; 3) tolerates significant soft occlusion (e.g. tree branches); and 4) associates, at a fundamental level, an estimated normal (i.e. no frontal-planar assumption) and texture with each surface patch.

A Neural Model of 3D Shape-From-Texture: Multiple-Scale Filtering, Boundary Grouping, and Surface Filling-In

A neural model is presented of how cortical areas V1, V2, and V4 interact to convert a textured 2D image into a representation of curved 3D shape. Two basic problems are solved to achieve this: (1) Patterns of spatially discrete 2D texture elements are transformed into a spatially smooth surface representation of 3D shape. (2) Changes in the statistical properties of texture elements across space induce the perceived 3D shape of this surface representation. This is achieved in the model through multiple-scale filtering of a 2D image, followed by a cooperative-competitive grouping network that coherently binds texture elements into boundary webs at the appropriate depths using a scale-to-depth map and a subsequent depth competition stage. These boundary webs then gate filling-in of surface lightness signals in order to form a smooth 3D surface percept. The model quantitatively simulates challenging psychophysical data about perception of prolate ellipsoids (Todd and Akerstrom, 1987, J. Exp. Psych., 13, 242). In particular, the model represents a high degree of 3D curvature for a certain class of images, all of whose texture elements have the same degree of optical compression, in accordance with percepts of human observers. Simulations of 3D percepts of an elliptical cylinder, a slanted plane, and a photo of a golf ball are also presented.

Multi-Area Visuotopic Map Complexes in Macaque Striate and Extra-striate Cortex

We propose that a simple, closed-form mathematical expression--the Wedge-Dipole mapping--provides a concise approximation to the full-field, two-dimensional topographic structure of macaque V1, V2, and V3. A single map function, which we term a map complex, acts as a simultaneous descriptor of all three areas. Quantitative estimation of the Wedge-Dipole parameters is provided via 2DG data of central-field V1 topography and a publicly available data set of full-field macaque V1 and V2 topography. Good quantitative agreement is obtained between the data and the model presented here. The increasing importance of fMRI-based brain imaging motivates the development of more sophisticated two-dimensional models of cortical visuotopy, in contrast to the one-dimensional approximations that have been in common use. One reason is that topography has traditionally supplied an important aspect of "ground truth", or validation, for brain imaging, suggesting that further development of high-resolution fMRI will be facilitated by this data analysis. In addition, several important insights into the nature of cortical topography follows from this work. The presence of anisotropy in cortical magnification factor is shown to follow mathematically from the shared boundary conditions at the V1-V2 and V2-V3 borders, and therefore may not causally follow from the existence of columnar systems in these areas, as is widely assumed. An application of the Wedge-Dipole model to localizing aspects of visual processing to specific cortical areas--extending previous work in correlating V1 cortical magnification factor to retinal anatomy or visual psychophysics data--is briefly discussed.