非均质储层叠前地震反演技术研究及应用

(1) I research on the relationship between elastic parameters, lithology and liquid. It is a physical base for pre-stack seismic inversion. I research all kinds of approximate expressions of Zoeppritz function. Then the relation of all kinds of approximate expressions can be confirmed. The geological model of water sand and gas sand in different depth was designed. Moreover I research on precision of all kinds of approximate expressions. (2) In process of seismic data which aim at amplitude recovery and apply in pre-stack seismic inversion, I advance to adopt double flow chart for different aim. Pre-stack noise elimination, real amplitude recovery and NMO correction of long offset are the key taches. (3) I made a systemic expatiate for the thinking and applicability about all kinds of expressions of elastic impedance. And mathematical model was applied to compare the precision with all kinds of expressions of elastic impedance. I propose a new pre-stack simultaneous inversion which is based on the Zoppritz function and simulated annealing algorithm. This method can ensure calculation precision of reflection coefficient from different incident angle and get a global optimum solution. Therefore this method improves the precision of pre-stack seismic inversion. (4) The object function of P-S wave pre-stack simultaneous inversion was established. I compared the precision and convergence between simultaneous inversion and P-wave inversion. And the results show that simultaneous inversion is superior to P-wave inversion. Through the study of AVO event of transformed wave, AVO characters of different kinds of gas sand were analyzed. (5) I carried out the study work of pre-stack seismic inversion for carbonate reservoir in middle of Tarim basin and sand shale reservoir in Sulige Area of Erdos Basin. The method and technology in this paper was applied to practical work. And I made a prediction for heterogeneous reservoir. Moreover it acquires a good application effect. Key Word: reflection coefficient， amplitude recovery， pre-stack seismic inversion， Heterogeneous reservoir，prediction.

残留盆地分布综合地球物理研究——以大港地区为例

Macro-distribution of residual basins is a basic question in residual basin research，the main object of macro-distribution study is to build strata framework, compute thickness of residual strata and analyze characteristics of residual basins. With the guidance of the theory of integrated geology and geophysical research, the paper assembled series of methods and established the technical chart based on gravity and magnetic data, with restriction of geology, seismic and drilling data. Based on potential field data processing and analysis, forward and inverse computation, region potential field analysis and potential field separation, etc. it computed depth of gravity/magnetic basement and got strata framework. It had got effective results in the research of macro-distribution of residual basin research in the Dagang area. It did the wavelet transform of gravity/magnetic data with multi-kind of wavelet basis using a trou algorithm. From comparison of processing result and their spectral of wavelet analysis, up continuation and filter method, the wavelet approximation is better to fit the regional potential field, and it is an effective method to separate gravity/magnetic effect caused by deep geology bodies. The experiment of matching pursuit shows that te transform domain methods have great advantage in potential data analysis. From the integrated geophysical study of rock property study, gravity/magnetic basement inversion and fault system analysis of the Dagang area, it gets the strata framework and the thickness of pre-Cenozoic residual strata. Comprehensive study with gravity and magnetotelluric profile inversion and interpretation, three prospect plays of macro-distribution of residual basins are fingered out. It has great residual strata thickness in the northern part of Chengning Uplift and there is thrust fault in the deep zone and good up-Paleozoic hydrocarbon source rocks in this area. With integrated analysis, this area will be the most prospective hydrocarbon location of pre-Cenozoic residual basins.

松辽盆地古龙地区葡萄花油层组沉积成岩相研究

The sediment and diagenesis process of reservoir are the key controlling factors for the formation and distribution of hydrocarbon reservoir. For quite a long time, most of the research on sediment-diagenesis facise is mainly focusing on qualitative analysis. With the further development on exploration of oil field, the qualitative analysis alone can’t meet the requirements of complicated requirements of oil and gas exploreation, so the quantitative analysis of sediment-diagenesis facise and related facies modling have become more and more important. On the basis of the research result from stratum and sediment on GuLong Area Putaohua Oil Layer Group, from the basic principles of sedimentology, and with the support from the research result from field core and mining research results, the thesis mainly makes the research on the sediment types, the space framework of sands and the evolution rules of diagenesis while mainly sticking to the research on sediment systement analysis and diagenetic deformation, and make further quantitative classification on sediment-diageneses facies qualitatively, discussed the new way to divide the sediment-diagenesis facies, and offer new basis for reservoir exploration by the research. Through using statistics theory including factor analysis, cluster analysis and discriminant analysis, the thesis devided sediment-diagenesis facies quantitatively. This research method is innovative on studying sediment-diagenesis facies. Firstly, the factor analysis could study the main mechanism of those correlative variables in geologic body, and then could draw a conclusion on the control factors of fluid and capability of reservoir in the layer of studying area. Secondly, with the selected main parameter for the cluster analysis, the classification of diagenesis is mainly based on the data analysis, thus the subjective judgement from the investigator could be eliminated, besides the results could be more quantitative, which is helpful to the correlative statistical analysis, so one could get further study on the quantitative relations of each sediment-diagenesis facies type. Finally, with the reliablities of discriminant analysis cluster results, and the adoption of discriminant probability to formulate the chart, the thesis could reflect chorisogram of sediment-diagenesis facies for planar analysis, which leads to a more dependable analytic results.According to the research, with the multi-statistics analysis methods combinations, we could get quantitative analysis on sediment-diagenesis facies of reservoir, and the final result could be more reliable and also have better operability.

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

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.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

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

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.

石油二次运移的物理与数值模拟实验及其应用

Facing the problems met in studies on predominant hydrocarbon migration pathways, experiments and numerical simulating were done in this thesis work to discuss the migration mechanisms. The aim is to analyze quantitatively the pathway pattern in basin scale and to estimate the hydrocarbon loss on the pathway that offer useful information for confirming the potential hydrocarbon accumulation. Based on our understandings on hydrocarbon migration and the fluid dynamic theory, a series of migration experiments were designed to observe the phenomena where kerosene is used as draining phase driven only by buoyancy force that expulses pore water. These experiments allow to study the formation of migration pathways, the distribution of non-wetting oil along these pathways, and the re-utilizing of previously existing pathways marked by residual traces etc. The types of pattern for migration pathways may be characterized by a phase diagram using two dimensionless numbers: the capillary number and the Bond number. The NMR technique is used to measure the average saturation of residual oil within the pathways. Based our experiment works and percolation concept, a numerical simulation model were proposed and realized. This model is therefore called as BP (Buoyancy Percolation) simulator, since buoyancy is taken as the main driving force in hydrocarbon migration. To make sure that BP model is applicable to simulate the process of oil secondary migration, the experimental phenomena are compared with those simulated with BP model by fractal method, and the result is positive. After then, we use BP simulator to simulate the process of migration of oil in the porous media saturated with water at different scale. And the results seem similar to those cited in literatures. In addition, our software is applied in Paris basin to predict the pathway of hydrocarbon migration happened in the Middle Jurassic reservoirs. It is found that the results obtained with our BP model are generally agree with Hindle (1997) and Bekeles'(1999), but our simulated migration pathway pattern and migration direction seem more reasonable than theirs.

Analysis of traumatic injuries presenting to a referral hospital emergency department in Moshi, Tanzania.

BACKGROUND: Injuries represent a significant and growing public health concern in the developing world, yet their impact on patients and the emergency health-care system in the countries of East Africa has received limited attention. This study evaluates the magnitude and scope of injury related disorders in the population presenting to a referral hospital emergency department in northern Tanzania. METHODS: A retrospective chart review of patients presenting to the emergency department at Kilimanjaro Christian Medical Centre was performed. A standardized data collection form was used for data abstraction from the emergency department logbook and the complete medical record for all injured patients. Patient demographics, mechanism of injury, location, type and outcomes were recorded. RESULTS: Ten thousand six hundred twenty-two patients presented to the emergency department for evaluation and treatment during the 7-month study period. One thousand two hundred twenty-four patients (11.5%) had injuries. Males and individuals aged 15 to 44 years were most frequently injured, representing 73.4% and 57.8%, respectively. Road traffic injuries were the most common mechanism of injury, representing 43.9% of injuries. Head injuries (36.5%) and extremity injuries (59.5%) were the most common location of injury. The majority of injured patients, 59.3%, were admitted from the emergency department to the hospital wards, and 5.6%, required admission to an intensive care unit. Death occurred in 5.4% of injured patients. CONCLUSIONS: These data give a detailed and more robust picture of the patient demographics, mechanisms of injury, types of injury and patient outcomes from similar resource-limited settings.

Mutations in NGLY1 cause an inherited disorder of the endoplasmic reticulum-associated degradation pathway.

PURPOSE: The endoplasmic reticulum-associated degradation pathway is responsible for the translocation of misfolded proteins across the endoplasmic reticulum membrane into the cytosol for subsequent degradation by the proteasome. To define the phenotype associated with a novel inherited disorder of cytosolic endoplasmic reticulum-associated degradation pathway dysfunction, we studied a series of eight patients with deficiency of N-glycanase 1. METHODS: Whole-genome, whole-exome, or standard Sanger sequencing techniques were employed. Retrospective chart reviews were performed in order to obtain clinical data. RESULTS: All patients had global developmental delay, a movement disorder, and hypotonia. Other common findings included hypolacrima or alacrima (7/8), elevated liver transaminases (6/7), microcephaly (6/8), diminished reflexes (6/8), hepatocyte cytoplasmic storage material or vacuolization (5/6), and seizures (4/8). The nonsense mutation c.1201A>T (p.R401X) was the most common deleterious allele. CONCLUSION: NGLY1 deficiency is a novel autosomal recessive disorder of the endoplasmic reticulum-associated degradation pathway associated with neurological dysfunction, abnormal tear production, and liver disease. The majority of patients detected to date carry a specific nonsense mutation that appears to be associated with severe disease. The phenotypic spectrum is likely to enlarge as cases with a broader range of mutations are detected.

Durability of antiretroviral therapy and predictors of virologic failure among perinatally HIV-infected children in Tanzania: a four-year follow-up.

BACKGROUND: In Tanzania, HIV-1 RNA testing is rarely available and not standard of care. Determining virologic failure is challenging and resistance mutations accumulate, thereby compromising second-line therapy. We evaluated durability of antiretroviral therapy (ART) and predictors of virologic failure among a pediatric cohort at four-year follow-up. METHODS: This was a prospective cross-sectional study with retrospective chart review evaluating a perinatally HIV-infected Tanzanian cohort enrolled in 2008-09 with repeat HIV-1 RNA in 2012-13. Demographic, clinical, and laboratory data were extracted from charts, resistance mutations from 2008-9 were analyzed, and prospective HIV RNA was obtained. RESULTS: 161 (78%) participants of the original cohort consented to repeat HIV RNA. The average age was 12.2 years (55% adolescents ≥12 years). Average time on ART was 6.4 years with 41% receiving second-line (protease inhibitor based) therapy. Among those originally suppressed on a first-line (non-nucleoside reverse transcriptase based regimen) 76% remained suppressed. Of those originally failing first-line, 88% were switched to second-line and 72% have suppressed virus. Increased level of viremia and duration of ART trended with an increased number of thymidine analogue mutations (TAMs). Increased TAMs increased the odds of virologic failure (p = 0.18), as did adolescent age (p < 0.01). CONCLUSIONS: After viral load testing in 2008-09 many participants switched to second-line therapy. The majority achieved virologic suppression despite multiple resistance mutations. Though virologic testing would likely hasten the switch to second-line among those failing, methods to improve adherence is critical to maximize durability of ART and improve virologic outcomes among youth in resource-limited settings.

Preoperative predictors of extended hospital length of stay following total knee arthroplasty.

The purpose of this study was to identify the preoperative predictors of hospital length of stay after primary total knee arthroplasty in a patient population reflecting current trends toward shorter hospitalization and using readily obtainable factors that do not require scoring systems. A single-center, multi-surgeon retrospective chart review of two hundred and sixty consecutive patients who underwent primary total knee arthroplasty was performed. The mean length of stay was 3.0 days. Among the different variables studied, increasing comorbidities, lack of adequate assistance at home, and bilateral surgery were the only multivariable significant predictors of longer length of stay. The study was adequately powered for statistical analyses and the concordance index of the multivariable logistic regression model was 0.815.

Orientation, Flow, and Clogging in a Two-Dimensional Hopper: Ellipses vs. Disks

Two-dimensional (2D) hopper flow of disks has been extensively studied. Here, we investigate hopper flow of ellipses with aspect ratio $\alpha = 2$, and we contrast that behavior to the flow of disks. We use a quasi-2D hopper containing photoelastic particles to obtain stress/force information. We simultaneously measure the particle motion and stress. We determine several properties, including discharge rates, jamming probabilities, and the number of particles in clogging arches. For both particle types, the size of the opening, $D$, relative to the size of particles, $\ell$ is an important dimensionless measure. The orientation of the ellipses plays an important role in flow rheology and clogging. The alignment of contacting ellipses enhances the probability of forming stable arches. This study offers insight for applications involving the flow of granular materials consisting of ellipsoidal shapes, and possibly other non-spherical shapes.

Clinical utility of a Web-enabled risk-assessment and clinical decision support program.

PURPOSE: Risk-stratified guidelines can improve quality of care and cost-effectiveness, but their uptake in primary care has been limited. MeTree, a Web-based, patient-facing risk-assessment and clinical decision support tool, is designed to facilitate uptake of risk-stratified guidelines. METHODS: A hybrid implementation-effectiveness trial of three clinics (two intervention, one control). PARTICIPANTS: consentable nonadopted adults with upcoming appointments. PRIMARY OUTCOME: agreement between patient risk level and risk management for those meeting evidence-based criteria for increased-risk risk-management strategies (increased risk) and those who do not (average risk) before MeTree and after. MEASURES: chart abstraction was used to identify risk management related to colon, breast, and ovarian cancer, hereditary cancer, and thrombosis. RESULTS: Participants = 488, female = 284 (58.2%), white = 411 (85.7%), mean age = 58.7 (SD = 12.3). Agreement between risk management and risk level for all conditions for each participant, except for colon cancer, which was limited to those <50 years of age, was (i) 1.1% (N = 2/174) for the increased-risk group before MeTree and 16.1% (N = 28/174) after and (ii) 99.2% (N = 2,125/2,142) for the average-risk group before MeTree and 99.5% (N = 2,131/2,142) after. Of those receiving increased-risk risk-management strategies at baseline, 10.5% (N = 2/19) met criteria for increased risk. After MeTree, 80.7% (N = 46/57) met criteria. CONCLUSION: MeTree integration into primary care can improve uptake of risk-stratified guidelines and potentially reduce "overuse" and "underuse" of increased-risk services.Genet Med 18 10, 1020-1028.

Nueva carta agroclimática de horas de frío en la Argentina

p.219-225