准噶尔盆地热历史与油气

Based on the temperature data from 196 wells and thermal conductivity measurements of 90 rock samples, altogether 35 heat flow data are obtained. The results show that the Junggar basin is a relatively "cold basin" at present. The thermal gradients vary between 11.6 and 26.5 ℃/km, and the thermal conductivity change from 0.17 to 3.6 W/mK. Heat flow ranges from 23.4 to 53.7 mW/m~2 with a mean of 42.3 ± 7.7 mW/m~2. The heat flow pattern shows that heat flow is higher on the uplifts and lower on the depressions. The overall low present-day heat flow in the Junggar Basin reflects its stable cratonic basement and Cenozoic tectonothermal evolution characterized by lithospheric thickening, thrust and fault at shallow crust as well as consequently quick subsidence during the Late Cenozoic. The study of the basin thermal history, which is one of the important content of the basin analysis, reveals not only the process of the basin's tectonothermal evolution, but also the thermal evolution of the source rocks based on the hydrocarbon generation models. The latter is very helpful for petroleum exploration. The thermal history of the Junggar basin has been reconstructed through the heat flow based method using the VR and Fission track data. The thermal evolutions of main source rocks (Permian and Jurassic) and the formations of the Permian and the Jurassic petroleum systems as well as the influences of thermal fields to petroleum system also have been discussed in this paper. Thermal history reconstruction derived from vitrinite reflectance data indicates that the Paleozoic formations experienced their maximum paleotemperature during Permian to Triassic with the higher paleoheat flow of around 70-85 mW/m~2 and the basin cooled down to the present low heat flow. The thermal evolution put a quite important effect on the formation and evolution of the petroleum system. The Jurassic petroleum system in the Junggar basin is quite limited in space and the source rocks of Middle-Lower Jurassic entered oli-window only along the foreland region of the North Tianshan belt, where the Jurassic is buried to the depth of 5-7 km. By contrast, the Middle-Lower Permian source rocks have initiated oil and gas generation in latter Permian to Triassic, and the major petroleum systems, like Mahu-West Pen 1 Well, was formed prior to Triassic when later Paleozoic formation reached the maximum paleotemperature.

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

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.

库车坳陷-南天山盆山系统第三纪脆性构造研究

Study on the structural coupling relationship between basin and range is not only helpful to recognize the basin formation and evolution systematically, but also to guide petroleum exploration in the basin. As a late Paleozoic Orogen, the South Tianshan Mountains reactivated and uplifted rapidly during the Cenozoic, and led to the Mesozoic-Cenozoic considerable thick deposits in the Kuqa Depression. The researches of the dissertation were carried out in the Kuqa depression-South Tianshan M ountain s ystem, a nd t he b rittle m icrotectonics w. ere c hosen as t he m ost important object. Based on observations and measurements of the field, we made detailed investigations on the geometry and kinematics of this area, and analyzed the abutting and cutting relationships and relative sequence of many brittle structures, such as joint, shear fractures, faults and some small-scale structures related to them closely. According to those brittle fractures' relationships with stress, the nature and variation of regional palaeostress field during the Cenozoic were studied through inversion of fault slip data and inferring stress state from joint sequences. And the deformation time was estimated primarily via ESR dating of faulting. Results show that the stress field varies as well in times as in space. The maximal principal stress direction shifted from the vertical to the horizontal, and stress regime from weak extension to strong compression from the Paleogene to the Neogene regionally. During the late Neogene, the structural deformation of the South Tianshan and the basin-range boundary was dominated by near N-S extension, while near N-S compressive deformation in the interior of the Kuqa Depression. There exits obvious differential stress state from the north to the south. ESR dating of the faulting during the Cenozoic indicates that, the normal faulting in the north edge of the Kuqa Depression have been active all along from the Miocene to the early Pleistocene, but the thrusting and reverse faulting in the interior only been active from the Pliocene to the early Pleistocene. On the base of those geological data and some geophysical information and theoretical calculation results, we infer that, the different stress regime the basin-range system is ascribed to the vertical uplift of the Tianshan Mountain. It was the vertical uplift that lead to the gravity-driven gliding of thick layers lying on the faulted basement from the South Tianshan Mountain to the Kuqa depression, and to folding and thrusting in the interior and frontal of the Kuqa depression. Combining the structural evolution with petroleum geological conditions of the Kuqa Depression, we think that the strong compressive deformation of the Kuqa Depression during rapid uplifting of the Tianshan Mountains from the Pliocene to the early Pleistocene play crucial role in the structural trap formation and proliferous gas accumulation.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

中国东海陆架盆地构造-热演化研究

East China Sea Shelf Basin (ECSSB), as a basin with prospect of oil & gas resource and due to its special geological location on the west margin of the west Pacific, attracts a lot of attention from many geologists in the world.Based on systematic temperature measurements, bottom hole temperature (BHT) and the oil temperature data, the geothermal gradients in the ECS SB are calculated and vary from 25 to 43°C/km, with a mean of 32.7°C/km. The geothermal gradient in Fuzhou Sag has the higher value(40.6°C/km) in Taibei Depression than that in others. The lower value (27.2 °C/km) occurs in in Xihu Depression. The middle values occurs in Jiaojiang and Lishui sags in Taibei Depression with a mean value of 34.8 °C/km. Incorporated with the measured thermal conductivity, heat flow values show that the ECSSB is characterized by present-day heat flow around 70.6mW/m2, varying between 55 and 88 mW/m2. No significant difference in heat flow is observed between the Xihu and the Taibei Depressions. These heat flow data suggest that the ECSSB is geothermally not a modem back-arc basin.Applying the paleogeothermal gradient based method, thermal history is reconstructed using vitrinite reflectance (VR) and apatite fission track (AFT) data. The results suggest that the thermal history was different in the Taibei and the Xihu depressions. Paleo-heat flow values when the pre-Tertiary formations experienced their maximum temperature at the end of the Paleocene reached a mean of 81 mW/m2 in the Taibei Depression, much higher than the present-day value. The lower Tertiary sediments in the Xihu Depression experienced maximum temperatures at the end of Oligocene and reached a mean paleo-heat flow value of 83.4 mW/m2. The time, when the paleo-heat flow reached the maximum value, suggests that the ECSSB rifted eastward.Tectonic subsidence analysis shows that the timing of the major rifting episode was different across the ECSSB. The rifting occurred from the Late Cretaceous to the early Eocene in the Taibei Depression, followed by thermal subsidence from the late Eocene to the end of Miocene. In contrast, in the Xihu Depression the initial subsidence lasted until the early Miocene and thermal subsidence to the end of Miocene. From Pliocene to the present, an accelerated subsidence took place all along the West Pacific margin of the east Asia.The thermal lithosphere thickness is determined by temperature profile in the lithosphere, the mantle adiabat or the dry basalt solidus. It indicates that the thermal lithosphere reached the thinnest thickness at the end of Eocene in the Taibei Depression and the end of Oligocene in the Xihu Depression, respectively, corresponding with a value of 57-66km and 56-64km. In Taibei Depression, the lithosphere thickness decreased 16-22km from the end of Mesozoic to Paleocene. After Paleocene, the thickness increased 13-16km and reached 71-79 km at present-day. In Xihu Depression, From the end of Oligocene to present-day, the thickness increased 10-13km and reached 69-76km at present-day. The evolution of the lithosphere thickness is associated closely with the lithosphere stretching.Combining the reconstructed thermal history and the burial history, the maturation of the Jurassic oil-source rock shows that the main hydrocarbon generation phase was in the mid-Jurassic and a secondary hydrocarbon generation occurred at the end of Paleocene. The secondary generation was controlled mainly by the tectono-thermal background during the Paleocene.

句容海安区块热史、成烃史与二次生烃

Jurong -Hai'an block of lower Yangtz area is one of the important petroleum exploration area among the residual marine basins in the south China. In the history of the basin's evolution, the strongly compressing, napping, folding and deforming during Indosinian to Yanshan epoch resulted in destruction of the early formed petroleum pool. Therefore, the strategy exploration of the secondary hydrocarbon generation and later formation of petroleum pool was brought forward for petroleum exploration and planning in the study area..On the basis of tectonic and sedimentary as well as present-day thermal regime evolution, using vitrinite reflectance and apatite fission track data and the paleo-heat flow based method, the basin's thermal history is reconstructed and hence the strata's temperature history are obtained. In addition, the maturation histories of the main four sets of marine hydrocarbon source rocks in the block are calculated. Furthermore, taking the maturity in the end of early Cretaceous as the original maturity and according to the formulas fitted by the secondary hydrocarbon generation model, the secondary hydrocarbon generation potential of the four sets of source rock is evaluated.The results of thermal history reconstruction show that Jurong-Hai'an block was under an uniform thermal setting during the Caledonian to Hercynian period and characterized by middle heat flow (52~57rnW/ m2). The uniform thermal setting was divided during and after Indosinian to Yanshan epoch. Wuwei area of southern Anhui province was under the high heat flow setting(~90 mW/m2)between 236Ma and 137Ma; Jurong area of southern Jiangsu ,Huangqiao area and Subei basin reached its maximum heat flow of 90,84 and 78-82 mW/m2 at 101Ma,157Ma and 56Ma respectivelyThe study of secondary hydrocarbon generation shows that the upper Paleozoic and Triassic source rocks have excellent secondary hydrocarbon generation potential. Silurian source rock posses some secondary hydrocarbon generation. Cambrian source rock, however, nearly has no secondary hydrocarbon generation. Overall there are no advantages of secondary hydrocarbon generation in the southern area of Jiansu. The intensity of secondary hydrocarbon generation in North Jiangsu basin is definitely better than that in the southern area of Jiangsu.

注意在视觉刺激位置编码中的作用

There is extensive agreement that attention may play a role in spatial stimmlus coding (Lu & Proctor, 1995). Some authors investigated the effects of spatial attention on the spatial coding by using spatial cueing procedure and spatial Stroop task. The finding was that the stroop effects were modulated by spatial cueing. Three hypotheses including attentional shift account, referential coding account, and event integration account were used to explain the modulation of spatial cueing over the spatial Stroop effects. In these previous studies, on validly cued trials, cue and target not only appeared at the same location, but also in the same object, which resulted in both location and object cued. Consequently, the modulation of spatial attentional cueing over spatial Stroop effects was confounded with the role of object-based attention. In the third chapter of this dissertation, using a modification of double rectangles cueing procedure developed by Egly, Driver and Rafal (1994) and spatial Stroop task employed by Lupiáñez and Funes (2005), separate effects of spatial attention and object-based attention on the location code of visual stimuli were investigated. Across four experiments, the combined results showed that spatial Stroop effects were modulate by object-based attention, but not by location-based attention. This pattern of results could be well explained by event integration account, but not by attentional shift account and referential coding account. In the fourth chapter, on the basis of the prior chapter, whether the modulation of attentional cueing on location code occurred at the stage of perceptual identification or response choice was investigated. The findings were that object-based attention modulated spatial Stroop effects and did not modulate the Simon effects, whereas spatial attention did not modulate Stroop and Simon effects. This pattern of results partially replicated the outcome of the previous chapter. The previous studies generally argued that the conflicts of spatial Stroop task and Simon task respectively occurred at at the stage of perceptual identification and response choice. Therefore, it is likely to conclude that the modulation of attention over spatial Stroop effect was mediated by object-based attention, and this modulation occurred at the stage perceptual identification. Considering that the previous studies mostly investigated the effects of attention captured by abrupt onset on the spatial Stroop effects, few studies investigated the effects of attention captured by offset cue on the spatial Stroop effects. The aim of the fifth chapter was to investigate the role of attention induced by offset and abrupt onset cue in the spatial Stroop task. These results showed that attention elicited by offset cue or abrupt onset cue modulated the spatial Stroop effects, which reconciled with event integration account.

Elliptic Object Detection Based on Shape Preserving and Active Contour

本文通过形状约束方程(组)与一般主动轮廓模型结合,将目标形状与主动轮廓模型融合到统一能量泛函模型中,提出了一种形状保持主动轮廓模型即曲线在演化过程中保持为某一类特定形状。模型通过参数化水平集函数的零水平集控制演化曲线形状,不仅达到了分割即目标的目的,而且能够给出特定目标的定量描述。根据形状保持主动轮廓模型,建立了一个用于椭圆状目标检测的统一能量泛函模型,导出了相应的Euler-Lagrange常微分方程并用水平集方法实现了椭圆状目标检测。此模型可以应用于眼底乳头分割,虹膜检测及相机标定。实验结果表明,此模型不仅能够准确的检测出给定图像中的椭圆状目标,而且有很强的抗噪、抗变形及遮挡性能。

A Method Based on Digital Image Analysis for Measuring the Strain Field of Deformed Sheet

对薄板成形应变场传统的测量方法进行了研究,指出了其不足和误差的来源,提出了数字图像分析法测量薄板成形中的应变场,对测量原理、新的测量方法对传统方法的改进,以及如何降低误差进行了介绍,指出数字图像分析法的前景,提出了改进意见。

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

Search For Critical Slip Surface In Slope Stability Analysis By Spline-Based Ga Method

The stability of a soil slope is usually analyzed by limit equilibrium methods, in which the identification of the critical slip surface is of principal importance. In this study the spline curve in conjunction with a genetic algorithm is used to search the critical slip surface, and Spencer's method is employed to calculate the factor of safety. Three examples are presented to illustrate the reliability and efficiency of the method. Slip surfaces defined by a series of straight lines are compared with those defined by spline curves, and the results indicate that use of spline curves renders better results for a given number of slip surface nodal points comparing with the approximation using straight line segments.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

A two-step phase-retrieval method in Fourier-transform ghost imaging

A two-step phase-retrieval method, based on Fourier-transform ghost imaging, was demonstrated. For the complex objects, the phase-retrieval process was divided into two steps: first got the complex object's amplitude from the Fourier-transform patterns of the squared object function, then combining with the Fourier-transform patterns of the object function to get the phase. The theoretical basis of this technique is outlined, and the experimental results are presented. (C) 2008 Elsevier B.V. All rights reserved.