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

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.

石油二次运移机理物理模拟研究

Oil and gas migration is very important for theoretical hydrocarbon geology study and exploration practice, but related research is weak. Physical simulation is a main method to study oil migration. Systematic experiments were done to quantitatively describe the migration patterns, path characters and oil saturation by adjusting the possible dynamic factors respectively. The following conclusions were drawn. 1. Darcy velocity and pore throat diameter were calculated according to seepage cross-sectional area and glass beads arrangement. With such normalized Darcy velocity and pore throat diameter, the date from one and two dimensional experiments can be reasonably drawn in two phase diagrams. It is found that the migration pattern can be identified using only one dimensionless number L which is defined as the ration of capillary number and Bond number. 2. Oil saturated in the pores between glass beads was used as calibration and oil saturation in the path was measured by magnetic resonance imaging. The results show that oil saturation in the center of migration path can reach 100%, is higher than oil saturation in the edge of migration path. 3. Percolation backbone during secondary oil migration was identified experimentally using Hele-Shaw cell. The backbone formed mainly because of the spatial variation of the cluster conductivity caused by oil saturation heterogeneity, main resistant force change, and path shrinkage and snap-off. Percolation backbone improves hydrocarbon migration efficiency and is a favorable factor for reservoir forming. 4. In the three dimensional filling models, the thickness of the secondary migration path is mall. It is only 2.5cm even for the piston pattern. Inclination of the model is the main influencing factor of the secondary path width.

松辽盆地南部岩性油藏成藏模式与分布规律

According to the basic geologic conditions, the paper is directed by the modem oil-gas accumulation theory and petroleum system in which typical oil pools are analyzed and the shape of lithologic trap and geologic factors are pointed out. The process during which oil and gas migrate from source rock to lithologic trap is rebuilt, and the accumulation model of oil pool is set up. With the comprehensive application of seismic geologic and log data and paying attention to the method and technology which is used to distinguish lithologic accumulation. Promising structural-lithofacies zones are got and the distribution rule of various lithologic accumulation is concluded. With making use of the biologic mark compound, different reservoirs are compared. As a result, the oil and gas in HeiDimiao come from Nenjiang Group's source rocks; in SaErTu from QingShenkou Group's and Nenjiang Group's, and in PuTaohua. GaoTaizi and FuYang from QingShankou Group's. According to the development and distribution of effective source rock, oil distribution and the comparison in the south of SongLiao basin, the characteristic of basin structure and reservoir distribution is considered, and then the middle-upper reservoir of SongLiao basin south are divided into two petroleum system and a complex petroleum system. Because of the characteristic of migration and accumulation, two petroleum systems can furtherly be divided into 6-7 sub-petroleum systems,20 sub-petroleum systems in all. As a result of the difference of the migration characteristic, accumulation conditions and the place in the petroleum system, the accumulation degree and accumulation model are different. So three accumulation mechanism and six basic accumulation model of lithologic trap are concluded. The distribution of lithologic pools is highly regular oil and gas around the generation sag distribute on favorable structural-lithofacies zones, the type of lithological pool vary regularly from the core of sandstone block to the upper zone. On the basic of regional structure and sedimentary evolution, main factors which control the form of trap are discovered, and it is the critical factor method which is used to discern the lithologic trap. After lots of exploration, 700km~2 potential trap is distinguished and 18391.86 * 10~4 tons geologic reserves is calculated. Oil-water distribution rule of pinch-out oil pool is put up on plane which is the reservoirs can be divided into four sections. This paper presented the law of distribution of oil and water in updip pinch-out reservoir, that is, hydrocarbon-bearing formation in plane can be divided into four zones: bottom edge water zone, underside oil and water zone, middle pure oil zone and above residual water zone. The site of the first well should be assigned to be middle or above pure oil zone, thus the exploration value of this type of reservoir can be recognized correctly. In accordance with the characteristics of seism and geology of low permeability thin sandstone and mudstone alternation layer, the paper applied a set of reservoir prediction technology, that is: (1)seism multi-parameter model identification; (2) using stratum's absorbing and depleting information to predict reservoir's abnormal hydrocarbon-bearing range. With the analysis of the residual resource potential and the research of two petroleum system and the accumulation model, promising objective zones are predicted scientifically. And main exploration aim is the DaRngZi bore in the west of ChangLin basin, and YingTai-SiFangZi middle-upper assembly in Honggang terrace.

柴达木盆地北缘地区油气成藏过程模拟与含油气系统评价

As an important part of petroleum exploration areas in the west of China, the north part of Qaidam basin is very promising in making great progress for petroleum discovery. But there are still many obstacles to overcome in understanding the process of petroleum formation and evaluation of oil & gas potential because of the complexity of geological evolution in the study area. Based upon the petroleum system theory, the process of petroleum formation is analyzed and the potential of oil & gas is evaluated in different petroleum systems by means of the modeling approach. The geological background for the formation of petroleum systems and the consisting elements of petroleum systems are described in detail. The thickness of strata eroded is estimated by means of vitrinite reflectance modeling, compaction parameter calculating and thickness extrapolating. The buried histories are reconstructed using the transient compaction model, which combines of forward and reverse modeling. The geo-history evolution consists of four stages - sedimentation in different rates with different areas and slow subsidence during Jurassic, uplifting and erosion during Cretaceous, fast subsidence during the early and middle periods of Tertiary, subsidence and uplifting in alternation during the late period of Tertiary and Quaternary. The thermal gradients in the study area are from 2.0 ℃/100m to 2.6 ℃/100m, and the average of heat flow is 50.6 mW/m~2. From the vitrinite reflectance and apatite fission track data, a new approach based up Adaptive Genetic Algorithms for thermal history reconstruction is presented and used to estimate the plaeo-heat flow. The results of modeling show that the heat flow decreased and the basin got cooler from Jurassic to now. Oil generation from kerogens, gas generation from kerogens and gas cracked from oil are modeled by kinetic models. The kinetic parameters are calculated from the data obtained from laboratory experiments. The evolution of source rock maturation is modeled by means of Easy %Ro method. With the reconstruction of geo-histories and thermal histories and hydrocarbon generation, the oil and gas generation intensities for lower and middle Jurassic source rocks in different time are calculated. The results suggest that the source rocks got into maturation during the time of Xiaganchaigou sedimentation. The oil & gas generation centers for lower Jurassic source rocks locate in Yikeyawuru sag, Kunteyi sag and Eboliang area. The centers of generation for middle Jurassic source rocks locate in Saishenteng faulted sag and Yuka faulted sag. With the evidence of bio-markers and isotopes of carbonates, the oil or gas in Lenghusihao, Lenghuwuhao, Nanbaxian and Mahai oilfields is from lower Jurassic source rocks, and the oil or gas in Yuka is from middle Jurassic source rocks. Based up the results of the modeling, the distribution of source rocks and occurrence of oil and gas, there should be two petroleum systems in the study area. The key moments for these two petroleum, J_1-R(!) and J_2-J_3, are at the stages of Xiaganchaigou-Shangyoushashan sedimentation and Xiayoushashan-Shizigou sedimentation. With the kinetic midels for oil generated from kerogen, gas generated from kerogen and oil cracked to gas, the amount of oil and gas generated at different time in the two petroleum systems is calculated. The cumulative amount of oil generated from kerogen, gas generated from kerogen and gas cracked from oil is 409.78 * 10~8t, 360518.40 * 10~8m~3, and 186.50 * 10~8t in J_1-R(!). The amount of oil and gas generated for accumulation is 223.28 * 10~8t and 606692.99 * 10~8m~3 in J_1-R(!). The cumulative amount of oil generated from kerogen, gas generated from kerogen and gas cracked from oil is 29.05 * 10~8t, 23025.29 * 10~8m~3 and 14.42 * 10~8t in J_2-J_3 (!). The amount of oil and gas generated for accumulation is 14.63 * 10~8t and 42055.44 * 10~8m~3 in J_2-J_3 (!). The total oil and gas potential is 9.52 * 10~8t and 1946.25 * 10~8m~3.

内蒙古乌拉特中旗北部宝音图群研究

The Baoyintu Group, lies in Wulate-zhongqi, Inner Mongolia, is a set of medium-grade metamorphic rock series which undergoes complex deformations. It consists of pelite schist, greenschist, plagioclase amphibolite quartzite and marble. The pelite schist is the main rock type and contains the classic medium pressure metamorphic minerals. The author divided Baoyintu group into five assemblages, investigated the rock association and plotted geological section of each assemblage in this area. Based on the systemically study of structural geology, petrology, geochemistry and mineralogy, the author reconstructs the protolith, sedimentary environment and tectonic evolution, discusses the mesoscopic and microscopic structure, metamorphism, geochemistry characters and the correlation between porphyroblast growth and deformation-metamorphism. There are three phase deformations in the research area: the earliest one occurred as the Baoyintu group deformed and metamorphosed and the main structure pattern is tight fold within layers during the Dl, large scale reversed fold and two phase faults (Fl fault and F2 fault) during the D2, and superimposed fold and F3 fault during D3. The F3 trancate the Wenduermian group of Silurian. The second and third phase deformation are relate to the orogenic event of late period of early Paleozoic. According to the rock association ,characteristics of the rocks and research of geochemistry, we get some information of the sedimentary environment and tectonic evolution of Baoyintu group. The source rocks are a set of terrigenous deposits-volcanic formation which reflect the history of the tectonic setting: stable- active-restable. And there are two sedimentary cycle from first assemblage to fifth assemblage: from first assemblage to fourth assemblage is a course of progression and the fifth assemblage is a start of regression. We also get the information of the P-T-t path by studying petrographies and calculating temperature and pressure. The path is not similiat to any classic type. And the interpretation is different from the traditional opinion. The P-T-t path reflects the dynamic course of convergence and uplift, magma underplating, back-arc extension and convergence of continental margin. Applying the theory of deformation partitioning to this area, the author discuss the relationship between deformation and porphyroblast growth, and get the conclusion of the sequence of deformation and metamorphism. At the first time we measure the distribution of chemical composition within the porphyroblast by XRF, confirm the theory of deformation partitioning quantitative and get new understanding about growth phase of porphyroblast and growth mode of porphyroblast: porphyroblast grow in the manner as "rose flower", the growth is controlled by the deformation. The elements distribution in porphyroblast reflects the growth manner and indicate history of metamorphism and deformation. So, we can deduce the metmorphism and deformation from the elements distribution in porphyroblast.

论大型拆离滑覆构造的特征、形成机制及其与油气成藏关系—以渤海湾盆地冀中坳陷饶阳凹陷为例

The studies of this paper is an important part of the "ninth five" science&technology-tackling project of CNPC -The oil and gas distribution regulation and the aims of explortion in jizhong depression. Basing on the former research results, with the materials of regional structural setting, major tectonic movements, bi-and tri-dimension seismic sections, oil well sections and reservoir sections, this paper involves studies of tectonic evolution, sedimentarv evolution, magma movement and reservoir prediction. The existence of huge stripping and gliding nappe is proved in the RaoYang Sag for the frist time. The properties, development, evolution and the relationship with reservoir of the stripping and gliding nappe are discussed in details in this paper. It is also talked about the affects of stripping and nappes to oil and gas exploration theoretically and practically in the paper. The marking attributes of the stripping and gliding nappe includes stripping and gliding plane, two deformation systems, stratigraphic repeat and hiatus close to the stripping and gliding plane, and the deformation attributes in the front and back of stripping and gliding nappe. The RaoYang stripping and gliding nappes can be divided into different belts in north-south direction and different zones in east-west direction. RaoYang Stripping and gliging nappes took place in the late Paleogene period and before the sedimentation of Neogene period. The sliding direction is NWW. The sliding distance is about 6km. The geothermal gradient in the separating slump area is low and stable. The formation of the stripping and gliding nappes is due to the regional structural setting, the sediments of Paleogene system, the soft roof and the uneven rising movement of structure units. The evolution of the stripping and gliding nappes can be divided into the following stages: regional differential elevation and subsidence, unstable gravity and gravitational sliding, the frist wholly stripping faults and sliding stage, and the following second and third stripping faults and sliding stages. The identification of RaoNan stripping and gliding nappes has an important role on the research of regional structure and oil and gas exploration. Basing on the properties of stripping and gliding nappes, we can identtify the gliding fractures, ductile compressional folds, the front and back structures of gliding nappes and gliding plane covered structures. Combination with different reservoir forming conditions, these structures can lead to different categories of reservoirs.

孤北洼陷沙三-沙四段隐蔽油藏形成机制及其预测

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in Gubei depression, this paper builds high resolution sequence stratigraphic structure, sedimentary system, sandbody distribution, the effect of tectonic in sequence and sedimentary system evolution and model of tectonic-lithofacies. The pool formation mechanism of subtle trap is developed. There are some conclusions and views as follows. 1.With the synthetic sequence analysis of drilling, seismic, and well log, the highly resolution sequence structure is build in Gubei depression. They are divided two secondary sequences and seven three-order sequences in Shahejie formation. They are include 4 kinds of system traces and 7 kinds of sedimentary systems which are alluvial fan, under water fan, alluvial fan and fan-delta, fan-delta, lacustrine-fan, fluvial-delta-turbidite, lakeshore beach and bar, and deep lake system. Sandbody distribution is show base on third order sequence. 2.Based on a lot of experiment and well log, it is point out that there are many types of pore in reservoir with the styles of corrosion pore, weak cementing, matrix cementing, impure filling, and 7 kinds of diagenetic facies. These reservoirs are evaluated by lateral and profile characteristics of diagenetic facies and reservoir properties. 3.The effect of simultaneous faulting on sediment process is analyzed from abrupt slope, gentle slope, and hollow zone. The 4 kinds of tectonic lithofacies models are developed in several periods in Gubei depression; the regional distribution of subtle trap is predicted by hydro accumulation characteristics of different tectonic lithofacies. 4.There are 4 types of compacting process, which are normal compaction, abnormal high pressure, abnormal low pressure and complex abnormal pressure. The domain type is normal compaction that locates any area of depression, but normal high pressure is located only deep hollow zone (depth more than 3000m), abnormal low pressures are located gentle slope and faulted abrupt slope (depth between 1200~2500m). 5.Two types dynamic systems of pool formation (enclosed and partly enclosed system) are recognized. They are composed by which source rocks are from Es3 and Es4, cap rocks are deep lacustrine shale of Esl and Es3, and sandstone reservoirs are 7 kinds of sedimentary system in Es3 and Es4. According to theory of petroleum system, two petroleum systems are divided in Es3 and Es4 of Gubei depression, which are high or normal pressure self-source system and normal or low pressure external-source system. 6.There are 3 kinds of combination model of pool formation, the first is litholgical pool of inner depression (high or normal pressure self-source type), the second is fault block or fault nose pool in marginal of depression (normal type), the third is fault block-lithological pool of central low lifted block (high or normal pressure type). The lithological pool is located central of depression, other pool are located gentle or abrupt slope that are controlled by lithological, faulting, unconfirmed. 7.This paper raise a new technique and process of exploration subtle trap which include geological modeling, coring description and logging recognition, and well log constrained inversion. These are composed to method and theory of predicting subtle trap. Application these methods and techniques, 6 hydro objects are predicted in three zone of depression.

New method for analysis of Chinese herbal complex prescription and its application

Traditional Chinese medicine (TCM) is a great treasure of China, the analysis of which is an arduous task. The viewpoint that all chemical constituents of Chinese herbal complex prescription should be analyzed as a black box is elucidated for the first time. Intelligent multi-mode multi-column chromatographic system (IMMCC) with its hybrids is the basic method and HPLC Unified Method is the breakthrough for the black box analysis. Dang-Gui-Bu-Xue-Tang was selected as a typical TCM and a systematic separation method from non-aqueous mobile phase to pure water mobile phase was put forward in order to convert unknown sample to known sample. The a, c values and UV spectra of 66 components of Astragalus, 78 components of Angelica and 71 components of Dang-Gui-Bu-Xue-Tang were obtained. Intelligent optimization and peak identification method and software for complex samples were developed and the optimum multi-step multi-binary gradient curve of mobile phase for Astragalus was ascertained. The maximum error and minimum error of predicted retention time for all components of Astragalus are 8.62% and 0.05% respectively. All components of Astragalus were compared with those of Angelica and it is found that many components of Astragalus are the same as those of Angelica, while the contents of these components are different. Many components of Dang-Gui-Bu-Xue-Tang are also the same as those of Astragalus and Angelica with different contents.