基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

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

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.

南海北部大陆边缘盆地深部地热特征用地球动力学演化

This report is a conclusion of the major research outcome during my post-doctoral residence of research and work. Its content covers the researches of the deep thermal characteristic and dynamics evolution beneath the northern margin basin of South China Sea. In this report, the each other action and effect between lithosphere ad mantle convection were regarded by the combine of deep and shallow study, subdivision from whole to part, and pay equal attention to determine the nature and fixed quantity. The investigative method we used in this report is geothermal and gravity methods. By the help of geological model and geophysics modeling, we calculated lithosphere thermal structure, rheology structure and mantle convection. Firstly, the report introduces concisely the purpose and the previous achievement to this research. Then, it analyzed the characteristic of heat flow on South China Sea. The structure of deep temperature and thermal has been calculated in some models of heat generation and conduction. The rock rheology structure also was computed by the relationship between temperature and viscosity. All these calculations were finished under the guidelines of combine with geology and geophysics. Meanwhile, the fields both deep mantle convection and small scale upper mantle convection are computed. Beside, the density and temperature disorder resulted by mantle convection were also computed with the convection field. After these, the report bring the contribution of local field of mantle convection, thermal construct and effective viscosity beneath the northern margin basin of South China Sea. And, base on the tectonic background and evolution feature, this report discussion the evolution mechanism of south China Sea and its northern margin basin. The end of this report, the main conclusion of this research was summarized and brings out.

复杂断块油藏形成机制及剩余油预测

The central uplift in the Huimin depression is famous for its large amounts of faults and small-scale fault-block area, and it is the famed typical complicated fault-block group oil & gas field in the whole world. After many years of rolling exploration and exploitation, many complex oil &gas field have been discovered in the central uplift, and won the splendent fruit. With the gradual deepening and development of the rolling exploitation, the exploration faces more and more difficulties. Therefore, it is important to reveal the forming mechanism and distributing rule of the complex fault-block reservoir, and to realize the forecast of the complex fault-block reservoir, sequentially, expedite the exploration step. This article applies the new multi-subject theory, method and technique such as structure geometry, kinematics, dynamics, structural stress field, fluid potential field, well logging record and constrained inversion of seismic records, coherence analysis, the seal mold and seal history of oil-bounded fault etc, and try to reveal the forming mechanism and distributing law of the complex fault-block reservoir, in result, implements the forecast of the fault-block reservoir and the remaining oil distributing. In order to do so, this article synthetically carries out structural estimate, reservoir estimate, fault sealing history estimate, oil-bearing properties estimate and residual. This article also synthetically researches, describes and forecast the complex fault-block in Huimin depression by use of the techniques, e.g. seismetic data post-stack processing technique, multi-component demarcating technique, elaborate description technique for the fault-block structure, technique of layer forecasting, fault sealing analysis technique, comprehensive estimate technique of fault-block, comprehensive analysis and estimate technique of remaining oil etc. The activities of the faults varies dramatically in the Huimin depression, and most of the second-class and the third-class faults are contemporaneous faults, which control the macroscopical distribution of the reservoir in the Huimin depression. The fourth-class faults cause the complication between the oil & gas among the fault-blocks. The multi-period strong activities of the Linyi fracture resulted in the vertical migration of large amount of oil & gas along with the faults. This is the main reason for the long vertical distribution properties near the Linyi fracture in the Huimin depression. The sealing ability of the fault is controlled by the property,size and direction of the main stress, the contact relationship of the both sides of the fault, the shale polluting factor, and the configuration relationship between the fault move period and the migration period of oil & gas. The article suggest four fault-sealing modes in the research zone for the first time, which establishes the foundation for the further forecast of the complex fault-block reservoir. Numerical simulation of the structural stress field reveals the distribution law and the evolvement progress of the three-period stress field from the end of the Dongying period to the Guantao period to nowadays. This article puts forward that the Linyi and Shanghe regions are the low value of the maximum main stress data. This is combined with the fault sealing history estimate, then multi-forming-reservoir in the central uplift is put forward. In the Shanghe oilfield, the article establishes six reservoir geological modes and three remaining oil distributing modes(the plane, the inside layer and the interlayer), then puts forward six increase production measure to enhance the remaining oil recovery ratio. Inducting the exploitation of oilfield, it wins notable economic effects and social effects.

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

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

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Towards Intelligent Structures: Active Control of Buckling

The buckling of compressively-loaded members is one of the most important factors limiting the overall strength and stability of a structure. I have developed novel techniques for using active control to wiggle a structural element in such a way that buckling is prevented. I present the results of analysis, simulation, and experimentation to show that buckling can be prevented through computer-controlled adjustment of dynamical behavior.sI have constructed a small-scale railroad-style truss bridge that contains compressive members that actively resist buckling through the use of piezo-electric actuators. I have also constructed a prototype actively controlled column in which the control forces are applied by tendons, as well as a composite steel column that incorporates piezo-ceramic actuators that are used to counteract buckling. Active control of buckling allows this composite column to support 5.6 times more load than would otherwise be possible.sThese techniques promise to lead to intelligent physical structures that are both stronger and lighter than would otherwise be possible.

Rapid Losses of Surface Elevation following Tree Girdling and Cutting in Tropical Mangroves.

Discussion Conclusions Materials and Methods Acknowledgments Author Contributions References Reader Comments (0) Figures Abstract The importance of mangrove forests in carbon sequestration and coastal protection has been widely acknowledged. Large-scale damage of these forests, caused by hurricanes or clear felling, can enhance vulnerability to erosion, subsidence and rapid carbon losses. However, it is unclear how small-scale logging might impact on mangrove functions and services. We experimentally investigated the impact of small-scale tree removal on surface elevation and carbon dynamics in a mangrove forest at Gazi bay, Kenya. The trees in five plots of a Rhizophora mucronata (Lam.) forest were first girdled and then cut. Another set of five plots at the same site served as controls. Treatment induced significant, rapid subsidence (−32.1±8.4 mm yr−1 compared with surface elevation changes of +4.2±1.4 mm yr−1 in controls). Subsidence in treated plots was likely due to collapse and decomposition of dying roots and sediment compaction as evidenced from increased sediment bulk density. Sediment effluxes of CO2 and CH4 increased significantly, especially their heterotrophic component, suggesting enhanced organic matter decomposition. Estimates of total excess fluxes from treated compared with control plots were 25.3±7.4 tCO2 ha−1 yr−1 (using surface carbon efflux) and 35.6±76.9 tCO2 ha−1 yr−1 (using surface elevation losses and sediment properties). Whilst such losses might not be permanent (provided cut areas recover), observed rapid subsidence and enhanced decomposition of soil sediment organic matter caused by small-scale harvesting offers important lessons for mangrove management. In particular mangrove managers need to carefully consider the trade-offs between extracting mangrove wood and losing other mangrove services, particularly shoreline stabilization, coastal protection and carbon storage.

A hardware proof of concept of a sailing robot for ocean observation

M J Neal, A hardware proof of concept of a sailing robot for ocean observation, IEEE Journal of Oceanic Engineering, 2006 accepted for publication RAE2008

Evidence for the Flare Trigger Site and Three-Dimensional Reconnection in Multiwavelength Observations of a Solar Flare

Fletcher, L., Metcalf, T.R., Alexander, D., Brown, D.S. and Ryder, L.A., 2001, Evidence for the flare trigger site and 3D reconnection in multi-wavelength observations of a solar flare, Astrophysical Journal, 554, 451-463.

Understanding student information behavior in relation to electronic information services: lessons from longitudinal monitoring and evaluation Part 1

Rowley, J.& Urquhart, C. (2007). Understanding student information behavior in relation to electronic information services: lessons from longitudinal monitoring and evaluation Part 1. Journal of the American Society for Information Science and Technology, 58(8), 1162-1174. Sponsorship: JISC

Children drawing toy commercials: re-imagining television production features

Griffiths, M. (2005). Children drawing toy commercials: re-imagining television production features. Visual Communication. 4(1), pp.21-37. RAE2008

Extremely long baseline interplanetary scintillation measurements of solar wind velocity

Breen, Andrew; Fallows, R. A.; Thomasson, P.; Bisi, M. M., 'Extremely long baseline interplanetary scintillation measurements of solar wind velocity', Journal of Geophysical Research (2006) 111(A8) pp.A08104 RAE2008

Biomass?Diversity Responses and Spatial Dependencies in Disturbed Tallgrass Prairies

Javier G. P. Gamarra and Ricard V. Sole (2002). Biomass-diversity responses and spatial dependencies in disturbed tallgrass prairies. Journal of Theoretical Biology, 215 (4) pp.469-480 RAE2008