Modeling the dynamic response of wrap-faced reinforced soil retaining walls

This paper describes the development of a numerical model for simulating the shaking table tests on wrap-faced reinforced soil retaining walls. Some of the physical model tests carried out on reinforced soil retaining walls subjected to dynamic excitation through uniaxial shaking tests are briefly discussed. Models of retaining walls are constructed in a perspex box with geotextile reinforcement using the wraparound technique with dry sand backfill and instrumented with displacement sensors, accelerometers, and soil pressure sensors. Results showed that the displacements decrease with the increase in number of reinforcement layers, whereas acceleration amplifications were not affected significantly. Numerical modeling of these shaking table tests is carried out using the Fast Lagrangian Analysis of Continua program. The numerical model is validated by comparing the results with experiments on physical models. Responses of wrap-faced walls with varying numbers of reinforcement layers are compared. Sensitivity analysis performed on the numerical models showed that the friction and dilation angle of backfill material and stiffness properties of the geotextile-soil interface are the most affecting parameters for the model response.

Numerical Simulation of Pulsed Laser Transformation Hardening

A novel pulsed laser surface processing technology is introduced, which can make use of the spatial and temporal profile of laser pulse to obtain ideal hardening parameters. The intensity distribution of laser pulse is spatially and temporally controlled by using laser shape transformation technology. A 3D numerical model including multi-phase transformations is established to explore material microstructure evolution induced by temperature field evolution. The influences of laser spatial-temporal profiles on hardening parameters are investigated. Different from the continuous laser processing technology, results indicate that spatial and temporal profiles are important factors in determining processing quality during pulsed laser processing method.

Análise de critérios de dimensionamento de carga em dutos de concreto instalados em vala na situação de recobrimento mínimo.

Perante diversas situações da engenharia são utilizadas formulações empíricas de dimensionamento baseadas em dados de campo e experiência profissional que definem muito o caráter subjetivo da metodologia padrão de projeto. O presente trabalho de pesquisa aborda os diversos métodos de obtenção dos esforços gerados em dutos enterrados submetidos a cargas dinâmicas e estáticas e sua posterior reavaliação através de modelagem numérica com o programa Plaxis 3D. Os métodos analíticos não convencionais foram comparados com o método padrão de cálculo sendo que o mesmo demonstrou ter uma boa precisão mesmo sem considerar outros fatores importantes como a parcela de resistência devida à coesão do solo e sua deformabilidade. A modelagem numérica demonstrou o conservadorismo do método de Marston e o subdmensionamento do espraiamento em prisma devido aos efeitos locais ocasionados pela adoção do recobrimento mínimo e sobrecarga dinâmica elevada. Também se observou, através da modelagem 3D, que a utilização dos dois métodos clássicos favorecem a obtenção de resultados dentro da razoabilidade.Verificou-se também, como resultado desta pesquisa, que a proposta de um método clássico modificado permite uma melhor aproximação da carga que atinge o duto.

Simulação numérica de dutos enterrados, submetidos à perda de apoio e elevação localizada.

Os recentes desastres ocorridos no país, como o rompimento da adutora em Campo Grande e os desastres relacionados às enchentes urbanas, mostram a necessidade de desenvolvimento de pesquisas científicas que auxiliem na compreensão e no dimensionamento das estruturas projetadas para atender a demanda da população. Os métodos analíticos e experimentais mais utilizados possuem algumas limitações de ordem teórica ou prática. Por outro lado, os métodos numéricos, capazes de simular etapas construtivas e envolver materiais com diferentes modelos constitutivos numa mesma análise, buscam atender às necessidades práticas dos projetos de geotecnia e, ao mesmo tempo, complementam os modelos analíticos e experimentais. Nesse trabalho foram realizadas comparações entre resultados obtidos em ensaios experimentais e resultados extraídos do modelo computacional, buscando aumentar a compreensão sobre a interação solo-estrutura em relação à distribuição de tensões mobilizadas e aos deslocamentos e deformações provocados. A simulação numérica foi feita com a utilização do PLAXIS/3D, software de análise geotécnica baseado no método dos elementos finitos. Os ensaios foram confeccionados na Escola de Engenharia de São Carlos/USP por Costa (2005) e envolveram dutos enterrados submetidos à perda de apoio ou elevação localizada. O estudo experimental foi realizado através de modelos físicos compostos por um maciço de areia pura, contendo um tubo repousando sobre um alçapão no centro do vão. Os modelos físicos foram equipados com instrumental capaz de medir as deflexões e as deformações específicas ao longo do duto, além das tensões totais no maciço de solo circundante e na base do equipamento.

Numerical modelling of dynamic response of loss modulated monolithic multisection mode-locked semiconductor lasers

Monolithic multisection mode-locked semiconductor lasers with an integrated distributed Bragg reflector (DBR) have recently been demonstrated to generate stable picosecond pulses at high repetition rates suitable for optical communication systems. However, there has been very little theoretical work on understanding the physical mechanisms of the device and on optimisation of the absorber modulator design. This article presents numerical modeling of the loss modulated mode-locking process in these lasers. The model predicts most aspects experimentally observed within this type of device, and the results show the output waveform, optical spectrum, instantaneous frequency chirp, and stable operating range.

On a classical problem of acoustic wave scattering by a free vortex: Numerical modelling

The scattering of sound from a point source by a Rankine vortex is investigated numerically by solving the Euler equations with the novel high-resolution CABARET method. For several Mach numbers of the vortex, the time-average amplitudes of the scattered field obtained from the numerical modeling are compared with the theoretical scaling laws' predictions. Copyright © 2009 by Sergey Karabasov.

Modeling nonlinearities in MEMS oscillators

We present a mathematical model of a microelectromechanical system (MEMS) oscillator that integrates the nonlinearities of the MEMS resonator and the oscillator circuitry in a single numerical modeling environment. This is achieved by transforming the conventional nonlinear mechanical model into the electrical domain while simultaneously considering the prominent nonlinearities of the resonator. The proposed nonlinear electrical model is validated by comparing the simulated amplitude¿frequency response with measurements on an open-loop electrically addressed flexural silicon MEMS resonator driven to large motional amplitudes. Next, the essential nonlinearities in the oscillator circuit are investigated and a mathematical model of a MEMS oscillator is proposed that integrates the nonlinearities of the resonator. The concept is illustrated for MEMS transimpedance-amplifier-based square-wave and sine-wave oscillators. Closed-form expressions of steady-state output power and output frequency are derived for both oscillator models and compared with experimental and simulation results, with a good match in the predicted trends in all three cases. © 1986-2012 IEEE.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

基于广义正交多项式迭积微分算子的地震波场数值模拟研究

Seismic Numerical Modeling is one of bases of the Exploratory Seismology and Academic Seismology, also is a research field in great demand. Essence of seismic numerical modeling is to assume that structure and parameters of the underground media model are known, simulate the wave-field and calculate the numerical seismic record that should be observed. Seismic numerical modeling is not only a means to know the seismic wave-field in complex inhomogeneous media, but also a test to the application effect by all kinds of methods. There are many seismic numerical modeling methods, each method has its own merits and drawbacks. During the forward modeling, the computation precision and the efficiency are two pivotal questions to evaluate the validity and superiority of the method. The target of my dissertation is to find a new method to possibly improve the computation precision and efficiency, and apply the new forward method to modeling the wave-field in the complex inhomogeneous media. Convolutional Forsyte polynomial differentiator (CFPD) approach developed in this dissertation is robust and efficient, it shares some of the advantages of the high precision of generalized orthogonal polynomial and the high speed of the short operator finite-difference. By adjusting the operator length and optimizing the operator coefficient, the method can involve whole and local information of the wave-field. One of main tasks of the dissertation is to develop a creative, generalized and high precision method. The author introduce convolutional Forsyte polynomial differentiator to calculate the spatial derivative of seismic wave equation, and apply the time staggered grid finite-difference which can better meet the high precision of the convolutional differentiator to substitute the conventional finite-difference to calculate the time derivative of seismic wave equation, then creating a new forward method to modeling the wave-field in complex inhomogeneous media. Comparing with Fourier pseudo-spectral method, Chebyshev pseudo-spectral method, staggered- grid finite difference method and finite element method, convolutional Forsyte polynomial differentiator (CFPD) method has many advantages: 1. Comparing with Fourier pseudo-spectral method. Fourier pseudo-spectral method (FPS) is a local operator, its results have Gibbs effects when the media parameters change, then arose great errors. Therefore, Fourier pseudo-spectral method can not deal with special complex and random heterogeneous media. But convolutional Forsyte polynomial differentiator method can cover global and local information. So for complex inhomogeneous media, CFPD is more efficient. 2. Comparing with staggered-grid high-order finite-difference method, CFPD takes less dots than FD at single wave length, and the number does not increase with the widening of the studying area. 3. Comparing with Chebyshev pseudo-spectral method (CPS). The calculation region of Chebyshev pseudo-spectral method is fixed in , under the condition of unchangeable precision, the augmentation of calculation is unacceptable. Thus Chebyshev pseudo-spectral method is inapplicable to large area. CFPD method is more applicable to large area. 4. Comparing with finite element method (FE), CFPD can use lager grids. The other task of this dissertation is to study 2.5 dimension (2.5D) seismic wave-field. The author reviews the development and present situation of 2.5D problem, expatiates the essentiality of studying the 2.5D problem, apply CFPD method to simulate the seismic wave-field in 2.5D inhomogeneous media. The results indicate that 2.5D numerical modeling is efficient to simulate one of the sections of 3D media, 2.5D calculation is much less time-consuming than 3D calculation, and the wave dispersion of 2.5D modeling is obviously less than that of 3D modeling. Question on applying time staggered-grid convolutional differentiator based on CFPD to modeling 2.5D complex inhomogeneous media was not studied by any geophysicists before, it is a fire-new creation absolutely. The theory and practices prove that the new method can efficiently model the seismic wave-field in complex media. Proposing and developing this new method can provide more choices to study the seismic wave-field modeling, seismic wave migration, seismic inversion, and seismic wave imaging.

冀东油田柳南地区新近系油藏剩余油形成及分布规律研究

This thesis is based on the research project of Study on the Geological Characteristics and Remaining Oil Distribution Law of Neogene Reservoirs in Liunan Area, which is one of the key research projects set by PetroChina Jidong Oilfield Company in 2006. The determination of remaining oil distribution and its saturation changes are the most important research contents for the development and production modification of oilfields in high water-cut phases. Liunan oilfield, located in Tangshan of Hebei Province geographically and in Gaoliu structural belt of Nanpu sag in Bohai Bay Basin structurally, is one of the earliest fields put into production of Jidong oilfield. Focusing on the development problems encountered during the production of the field, this thesis establishes the fine geological reservoir model through the study of reservoir properties such as fine beds correlation, sedimentary facies, micro structures, micro reservoir architecture, flow units and fluid properties. Using routine method of reservoir engineering and technology of reservoir numerical modeling, remaining oil distribution in the target beds of Liunan area is predicted successfully, while the controling factors of remaining oil distribution are illustrated, and the model of remaining oil distribution for fault-block structure reservoirs is established. Using staged-subdivision reservoir correlation and FZI study, the Strata in Liunan Area is subdivided step by step; oil sand body data-list is recompiled; diagram databases are established; plane and section configuration of monolayer sandstone body, and combination pattern of sandstone bodys are summarized. The study of multi-level staged subdivision for sedimentary micro-facies shows that the Lower member of Minghuazhen formation and the whole Guantao formation in Liunan Area belong to meandering river and braided river sedimentary facies respectively, including 8 micro facies such as after point bar, channel bar, channel, natural levee, crevasse splay, abandoned channel, flood plain and flood basin. Fine 3D geological modeling is performed through the application of advanced software and integration of geological, seismic logging and reservoir engineering data. High resolution numerical simulation is performed with a reserve fitting error less than 3%, an average pressure fitting fluctuation range lower than 2Mpa and an accumulate water cut fitting error less than 5%. In this way, the distribution law of the target reservoir in the study area is basically recognized. Eight major remaining oil distribution models are established after analysis of production status and production features in different blocks and different layers. In addition, fuzzy mathematics method is used to the integreted evaluation and prediction of abundant remaining oil accumulation area in major production beds and key sedimentary time units of the shallow strata in Liunan Area and corresponding modification comments are put forward. In summary, the establishment of fine reservoir geological model, reservoir numerical simulation and distribution prediction of remaining oil make a sound foundation for further stimulation of oilfield development performance.

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

塔里木盆地轮南地区三叠系油气地质模型及油气预测研究

As an important measure to understand oil and gas accumulation during petroleum exploration and development, Petroleum geological model is an integrated system of theories and methods, which includes sedimentology, reservoir geology, structural geology, petroleum geology and other geological theories, and is used to describe or predict the distribution of oil and gas. Progressive exploration and development for oil and gas is commonly used in terrestrial sedimentary basin in China for the oil and gas generation, accumulation and exploitation are very intricate. It is necessary to establish petroleum geological model, adaptive to different periods of progressive exploration and development practice. Meanwhile there is lack of an integrated system of theories and methods of petroleum geological model suitable for different exploration and development stages for oil and gas, because the current different models are intercrossed, which emphasize their different aspects. According to the characteristics of exploration and development for the Triassic oil and gas pool in Lunnan area, Tarim Basin, the Lunnan horst belt was selected as the major study object of this paper. On the basis of the study of petroleum geological model system, the petroleum geological models for different exploration and development stages are established, which could be applied to predict the distribution of oil and gas distribution. The main results are as follows. (1) The generation-accumulation and exploration-development of hydrocarbon are taken as an integrated system during the course of time, so petroleum exploration and development are closely combined. Under the guidance of some philosophical views that the whole world could be understood, the present writer realizes that any one kind of petroleum geological models can be used to predict and guide petroleum exploration and development practice. The writer do not recognize that any one kind of petroleum geological models can be viewed as sole model for guiding the petroleum exploration and development in the world. Based on the differences of extents and details of research work during various stage of exploration and development for oil and gas, the system of classification for petroleum geological models is established, which can be regarded as theoretical basis for progressive petroleum exploration and development. (2) A petroleum geological model was established based on detailed researches on the Triassic stratigraphy, structure, sedimentology and reservoir rocks in the Lunnan area, northern Tarim Basin. Some sub-belt of hydrocarbon accumulation in the Lunnan area are divided and the predominate controlling factors for oil and gas distribution in the Lunnan area are given out. (3) Geological models for Lunnan and Jiefangqudong oil fields were rebuilt by the combinations of seismology and geology, exploration and development, dynamic and static behavior, thus finding out the distribution of potential zones for oil and gas accumulations. Meanwhile Oil and gas accumulations were considered as the important unit in progressive exploration and development, and the classification was made for Lunnan Triassic pools. Petroleum geological model was created through 3D seismic fine interpretation and detailed description of characteristics of reservoir rocks and the distribution of oil and gas, especially for LN3 and LN26 well zones. The possible distribution of Triassic oil traps and their efficiency in the Lunnan area has been forecasted, and quantitative analysis for original oil(water) saturation in oil pools was performed. (4) The concept of oil cell is proposed by the writer for the first time. It represents the relatively oil-rich zones in oil pool, which were formed by the differences of fluid flows during the middle stage of reservoir development. The classification of oil cells is also given out in this paper. After the studies of physical and numerical modeling, the dominant controlling factors for the formation of various oil cells are analyzed. Oil cells are considered as the most important hydrocarbon potential zones after first recovery, which are main object of progressive development adjustment and improvement oil recovery. An example as main target of analysis was made for various oil cells of Triassic reservoir in the LN2 well area. (5) It is important and necessary that the classification of flow unit and the establishment of geological model of flow unit based on analysis of forecast for inter-well reservoir parameters connected with the statistical analysis of reservoir character of horizontal wells. With the help of self-adaptive interpolation and stochastic simulation, the geological model of flow units was built on the basis of division and correlation of flow units, with which the residual oil distribution in TIII reservoir in the LN2 well area after water flooding can be established.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

孤南-三合村洼陷含油气系统研究

The development petroleum geology has made people from studying and studying and predicting in statically and respectively the pool-forming conditions of an area such as oil source bed, reservoir, overlying formation, migration, trap and preservation, etc. to regarding these conditions as well as roles of generation, reservation and accumulation as an integrated dynamic evolution development system to do study .Meanwhile apply various simulating means to try to predict from quantitative angle. Undoubtedly, the solution of these questions will accumulate exploration process, cut down exploration cost and obtain remarkable economic and social benefits. This paper which take sedimentology ,structural geology and petroleum geology as guides and take petroleum system theory as nucleus and carry out study thinking of beginning with static factor and integration of point and face as well as regarding dynamic state factor as factor and apply study methods of integration of geology, Lab research and numerical modeling proceed integrated dissect and systematic analysis to GuNan-SanHeCun depression. Also apply methods of integration of sequence stratigraphy, biostratigraphy, petrostratigraphy and seismic data to found the time-contour stratigraphic framework and reveal time-space distribution of depositional system and meantime clarify oil-source bed, reservoir and overlying distribution regular patterns. Also use basin analysis means to study precisely the depositional history, packed sequences and evolution. Meanwhile analyze systematically and totally the fracture sequence and fault quality and fault feature, study the structural form, activity JiCi and time-space juxtaposion as well as roles of fault in migration and accumulation of oil and gas of different rank and different quality fault. Simultaneously, utilize seismic, log, analysis testing data and reservoir geology theory to do systematic study and prediction to GuNan-SanHeCun reservoir, study the reservoir types macroscopic distribution and major controlling factors, reservoir rock, filler and porosity structural features as well as distribution of reservoir physical property in 3D space and do comprehensive study and prediction to major controlling and influential factors of reservoir. Furthermore, develop deepingly organic geochemistry comprehensive study, emphasis on two overlaps of oil source rock (ESI, ES3) organic geochemistry features, including types, maturity and spatial variations of organic matter to predict their source potential .Also apply biological marks to proceed oil-to-source correlation ,thereby establish bases for distribution of petroleum system. This study recover the oil generation history of oil source rocks, evaluate source and hydrocarbon discharge potential ,infer pool-forming stages and point out the accumulation direction as well as discover the forming relations of mature oil-source rock and oil reservoir and develop research to study dynamic features of petroleum system. Meanwhile use systematic view, integrate every feature and role of pool forming and the evolution history and pool-forming history, thereby lead people from static conditions such as oil source bed, reservoir, overlying formation, migration, trap and preservation to dynamically analyzing pool-forming process. Also divide GuNan-SanHeCun depression into two second petroleum system, firstly propose to divide second petroleum system according to fluid tress, structural axis and larger faults of cutting depression, and divide lower part of petroleum system into five secondary systems. Meanwhile establish layer analysis and quantitative prediction model of petroleum model, and do quantitative prediction to secondary petroleum system.

高含水油藏复合驱剩油分布研究-以羊三木油田馆陶组油藏复合驱为例

Based on the study of the combined flooding test block of Guantao formation in Third faulted block of Yangsanmu oil field, this paper carries out the integration of reservoir precise characterization for very high water cut reservoir, establishes precise 3D geologic model for high water cut development period and states the changing law of the reservoir architecture dtiring development by combined flooding. Then, by subdivided the thick oil reservoir, the study of remaining oil saturation monitoring in fiber glass cased well and tracer monitoring is developed. According the study of multiple constrained combined flooding reservoir numerical simulation, remaining oil distribution are predicted, the methods architecture of predicting remaining oil distribution are established for fluvial facies reservoir at late development stage, develops plan is designed and adjustment associating technologies for enhancing oil recovery. On these base, related measures for tapping the potential are given, it is verified and optimized through the field former test and the good economic effect is achieved . The major achievements of this paper are as follows. The changing law of the reservoir architecture and it's property parameters is revealed, The result indicates that the temperature-pressure of the injecting material and the interaction effect of the injecting material and reservoir petrography are the main factors of the dynamic changes of the reservoir architecture. The quantitative reservoir geologic model, which is tallied with dynamic reservoir parameters of the study area, is established. Subdivided the thick oil reservoir is very important for the study of the remaining oil distribution within the thick oil reservoir. Subdivided the thick oil reservoir technology, which consists of six technologies as follow: micro-cyclic divided, flow unit method, architectural element method, high resolution log technology, high resolution-process technology for normal logging data and using the production data is presented. 3. It is established dynamic monitoring system of remaining oil saturation quantitative research which are inner and interlayer remaining oil saturation from time-lapse logging in fiber glass cased well, inter-well remaining oil saturation from the technology of isotopic tracer monitoring technology, and 4d remaining oil saturation distribution from combined flooding numerical modeling integrated by production datao The forming mechanism of remaining oil for polymer flooding and alkali/polymer combined flooding is clarified, and the plane and vertical distribution law of remaining oil after combined flooding is revealed. Predicting methods and technologies for the combined flooding reservoir of fluvial facies is developed. Combined flooding has been achieved good displacement result in the pilot of Third fault block in Yangsanmu oil field, and accumulated types of important parameters and optimum plans, this technology of combined flooding is expected to increase recovery ratio by 4.77%.