Dielectric responses of anisotropic graded granular composites having arbitrary inclusion shapes

The transformation field method (TFM) originated from Eshelby's transformation field theory is developed to estimate the effective permittivity of an anisotropic graded granular composite having inclusions of arbitrary shape and arbitrary anisotropic grading profile. The complicated boundary-value problem of the anisotropic graded composite is solved by introducing an appropriate transformation field within the whole composite region. As an example, the effective dielectric response for an anisotropic graded composite with inclusions having arbitrary geometrical shape and arbitrary grading profile is formulated. The validity of TFM is tested by comparing our results with the exact solution of an isotropic graded composite having inclusions with a power-law dielectric grading profile and good agreement is achieved in the dilute limit. Furthermore, it is found that the inclusion shape and the parameters of the grading profile can have profound effect on the effective permittivity at high concentrations of the inclusions. It is pointed out that TFM used in this paper can be further extended to investigate the effective elastic, thermal, and electroelastic properties of anisotropic graded granular composite materials.

机器人标定技术研究

机器人标定是离线编程技术实用化的关键技术之一 ,测量手段是制约其发展的一个重要因素 .本文提出了一种简单实用的测量方法 .同时建立了具有完备性、比例性、连续性的 RH6 - A型机器人的真实运动学模型 .并根据测量结果 ,成功地辨识出了机器人的几何结构参数 ,提高了机器人的绝对精度 .

基于视觉伺服的输电线巡检机器人抓线控制

巡检机器人在越障时,需要完成机器人手臂的准确抓线控制.结合输电线的几何特征和摄像机成像原理,提出了一种基于单摄像机的立体视觉方法来确定输电线的位置和姿态.基于该定位方法及视觉伺服理论,建立机械手抓线伺服控制模型.利用自行研制的巡检机器人进行了视觉伺服抓线实验;实验结果验证了该方法的有效性.

基于TOA的无线传感器网络自定位技术的研究

提出一种基于TOA(tim e of arrival)测距、利用节点之间的几何关系实现传感器网络中未知节点自定位的算法.该方法计算量小、算法简单且定位精度高.仿真实验证明了该方法的有效性.*

一种基于线性模型的机器人手眼标定新方法

针对基于传统手眼标定模型AX=XB的标定问题,引入了由摄像机投影矩阵M构成的数学模型MY=M'YB,并针对该模型提出了一种新解法。该方法运用矩阵直积和特征向量理论,将关于机器人手眼关系矩阵X的问题描述转换成关于关联矩阵Y的线性方程,运用最小二乘法求得线性闭解。该方法消除了传统方法多次分解摄像机传感器系统的内外参数而引入的摄像机传感器系统的固有误差和随机误差的影响。实验表明该方法可以提高手眼系统标定精度。

可重构星球探测机器人的运动学建模及轨迹规划

提出了可重构星球探测机器人的概念,对系统中子机器人的研究进行了重点论述.通过设计恰当的子机器人连杆坐标系,利用Denavit Hartenberg方法完成了子机器人的运动学建模,并直接给出了子机器人的运动学正解模型.由于使用单一的求解算法不能求出工作空间的封闭解,因此综合利用代数法、几何法原理及空间投影关系,结合子机器人的结构特殊性推导出了运动学逆解,从而得到了工作空间内的所有解.在此基础上,考虑结构间的约束关系,给出了子机器人的工作空间及轨迹规划方法.最后,使用OpenGL对设计的子机器人系统进行了运动学仿真实验,实验以末端操作器的直线运动为例,充分考虑空间几何的关系,其结果有效地证明了建模及轨迹规划的正确性.

全方位六足步行机器人运动规划的相对运动算法

本文提出了步行机器人运动控制算法。该方法以相对运动学原理为基础,把机体的运动规划问题转化为腿的足端轨迹规划问题,从而使步行机器人运动控制问题得到大大简化.并应用该方法对全方位三角步态算法及稳定性进行分析求解.

小湾电站坝基岩体卸荷工程地质力学研究

Dam is the key main works in the construction of water power. The success or failure of the construction of the dam mainly depends on the stability of the dam foundation. The double curvature arch dam-XiaoWan Dam is the highest one among the dams with the same type in the world, and the water thrust acted on it reaches ton, so the rock bearing capacity of dam foundation becomes more important. Because of the high and steep valley-side slope, the large scale of excavation and the complex body type of excavation, it is prominent that the problem of stress release of the rock mass in dam foundation. More great attentions should be paid for the stability and the degraded of rock properties of rock mass induced by the stress release. In this paper, the phenomena of stress release of rock mass in XiaoWan Dam foundation and its mechanisms were analyzed based on the collection of data, the detailed field engineering investigations, measurement of the rock mass and the 2D numerical calculations. The rock mass under the foundation is weak-weathered to intact, the quality of which is good. After excavation of the foundation, the rock mass near the slope surface occurred extend, stretch and stick-slip along original textures till the new fracture surface formed. Then platy structure of the rock mass takes on. The rock mass in the dam foundation occur resilience due to stress release towards free faces with the characteristics of time effect and localized deformation. In-situ measurements show that the rock mass near the surface are degraded. The stress release induced by excavation is a process of the interaction between engineering structures and geologic body. The stress release of rock mass in dam foundation is related to the changed degree of geometrical conditions. The rock near excavation surface failed nearly under uniaxial stresses. The bending-breaking mechanism of plate girder can interpret the failure model of the rock mass with platy structure in dam foundation slope. In essence, the stress release is the change of stress field including the change of directions and magnitudes of stress induced by excavation, which can induce the comedown of the safety margin. In this paper, the inducing conditions of stress release were calculated by numerical analyses. Moreover, from the point of view that the change of stress field, the coefficient of K, i.e. the variable load coefficient was proposed. Then the law of the change of it is interpreted. The distributional characteristics of fracture zone were expressed by the coefficient. The stress release of hard rock has the characteristic of localization. The measuring technique of sound wave can not reflect the small cracks in this kind of rock mass due to stress release. So, the spectral analysis method was proposed. At the same time, the application foreground in engineering of the Stockwell Time-Frequency- Spectrum method was discussed with a view to the limitation of it.

悬索桥隧道式锚碇与围岩系统的破坏模式研究

The anchorages are unparalleled structures only in a suspension bridge, and as main bearing facilities, play an important role in connecting the superstructures and the ground. The tunnel anchorage, as one alternative type of the anchorages, has more advantages over its counterpart, the gravity anchorage. With the tunnel anchorages adopted, not only can surface excavation be reduced to protect the environment, and natural condition of the rock be utilized and potential bearing capacity of surrounding rock be mobilized to save engineering cost, but also the technological predominance of auxiliary engineering measures, such as prestressed concrete, anchoring piles, rock anchors and collar beam between the two separated anchorages, can be easily cooperated to work together harmoniously under the circumstances of poor rock quality. There are plentiful high mountains and deep canyons in west part of China, and long-span bridge construction is inevitably encountered in order to realize leapfrogging development of the transportation infrastructure. Western mountainous areas usually possess the conditions for constructing tunnel anchorages, and therefore, the tunnel anchorages, which are conformed to the conception of resource conservative and sustainable society, extremely have application and popularization value in western underdeveloped region. The scientific and technological problem about the design, construction and operation of tunnel anchorages should be further investigated. Combining the engineering of western tunnel anchorages for the Balinghe Suspension Bridge, this paper probed into the survey method and in-situ test method for tunnel anchorages, scientific rock quality evaluation of surrounding rock to provide reasonable physical and mechanical parameters for design, construction and operation of tunnel anchorages, bearing capacity estimation for tunnel anchorage, deformation prediction of the anchorage-rockmass system, tunnel-anchorage slope stability analysis and the evaluation of excavation stability and degree of safety of the anchorage tunnel. The following outcomes were obtained: 1. Materials of tunnel anchorages of suspension bridge built (and in progress) at home and abroad were systematically sorted out, with the engineering geological condition and geomechanical property of surrounding rock around the anchorage tunnel, the design size of anchorages and the construction method of anchorage tunnel paid more emphasis on, to unveil the internal relationship between the engineering geological conditions of surrounding rock and the design size and axis angle of anchorages and provide references for future design, construction and study of tunnel anchorages. 2. Physical and mechanical parameters were recommended based on three domestic and foreign methods of rock quality evaluation. 3. In-situ tests, adopting the back-thrust method, of two kinds of reduced scale model, 1/30 and 1/20, for the tunnel anchorages were conducted in the declining exploration drift with rock mass at the test depth being the same as surrounding rock around real anchorages, and reliable field rockmass displacement data were acquired. Attenuation relation between the increment of distance from the anchorage and the decrement of rockmass displacement under maximum test load, and influential scope suffered by anchorage load were obtained. 4. Using similarity theory, the magnitude of real anchorage and rockmass displacement under design load and degree of safety of the anchorage system were deduced. Furthermore, inversion analysis to deformation modulus of slightly weathered dolomite rock, the surrounding rock of anchorage tunnel, was performed by the means of numerical simulation. 5. The influential law of the geometrical size to the limit bearing capacity of tunnel anchorage was studied. 6. Based on engineering geological survey data, accounting for the combination of strata layer and adverse discontinuities, the failure patterns of tunnel anchorage slope were divided into three modes: sliding of splay saddle pier slope, superficial-layer slippage, and deep-layer slippage. Using virtual work principle and taking anchorage load in account, the stability of the three kinds of failure patterns were analyzed in detail. 7. The step-by-step excavation of anchorage tunnel, the numerical overload and the staged decrement of rock strength parameters were numerically simulated to evaluate the excavation stability of surrounding rock around anchorage tunnel, the overload performance of tunnel anchorage, and the safety margin of strength parameters of the surrounding rock.

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

纵波方位属性裂缝检测技术综合研究

Fractured oil and gas reservoir is an important type of oil and gas reservoir, which is taking a growing part of current oil and gas production in the whole world. Thus these technologies targeted at exploration of fractured oil and gas reservoirs are drawing vast attentions. It is difficult to accurately predict the fracture development orientation and intensity in oil and gas exploration. Focused on this problem, this paper systematically conducted series study of seismic data processing and P-wave attributes fracture detection based on the structure of ZX buried mountain, and obtained good results. This paper firstly stimulated the propagation of P-wave in weak anisotropic media caused by vertical aligned cracks, and analyzed the rule of P-wave attributes’ variation associated with observed azimuth, such as travel-time, amplitude and AVO gradient and so on, and quantitatively described the sensitive degree of these attributes to anisotropy of fracture medium. In order to further study the sensitive degree of these attributes to anisotropy of fractures, meanwhile, this paper stimulated P-wave propagation through different types and different intensity anisotropic medium respectively and summarized the rule of these attributes’ variation associated with observed azimuth in different anisotropic medium. The results of these studies provided reliable references for predicting orientation, extensity and size of actual complicated cracked medium by P-wave azimuth attributes responses. In the paper, amounts of seismic data processing methods are used to keep and recover all kinds of attributes applied for fracture detection, which guarantee the high accurate of these attributes, thus then improve the accurate of fracture detection. During seismic data processing, the paper adopted the three dimensional F-Kx-Ky field cone filter technique to attenuate ground roll waves and multiple waves, then enhances the S/N ratio of pre-stack seismic data; comprehensively applying geometrical spread compensation, surface consistent amplitude compensation, residual amplitude compensation to recover amplitude; common azimuth processing method effectively preserves the azimuthal characteristics of P-wave attributes; the technique of bend ray adaptive aperture pre-stack time migration insures to obtain the best image in each azimuth. Application of these processing methods guaranteed these attributes’ accuracy, and then improved the accuracy of fracture detection. After comparing and analyzing a variety of attributes, relative wave impedance (relative amplitude) attribute is selected to inverse the orientation of fracture medium; attenuation gradient and corresponding frequency of 85% energy are selected to inverse the intensity of fracture medium; then obtained the fracture distribution characteristics of lower Paleozoic and Precambrian in ZX ancient buried mountains. The results are good accord with the characteristics of faults system and well information in this area.

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.