Changes in Petrophysical Properties of the Stone Surface due to Past Conservation Treatments in Archaeological Sites of Merida (Spain)

Petrophysical properties, such as porosity, permeability, density or anisotropy de-termine the alterability of stone surfaces from archaeological sites, and therefore, the future preservation of the material. Others, like superficial roughness or color, may point out changes due to alteration processes, natural or man-induced, for ex-ample, by conservation treatments. The application of conservation treatments may vary some of these properties forcing the stone surface to a re-adaptation to the new conditions, which could generate new processes of deterioration. In this study changes resulting from the application of consolidating and hydrophobic treatments on stone materials from the Roman Theatre (marble and granite) and the Mitreo’s House (mural painting and mosaics), both archaeological sites from Merida (Spain), are analyzed. The use of portable field devices allows us to perform analyses both on site and in la-boratory, comparing treated and untreated samples. Treatments consisted of syn-thetic resins, consolidating (such as tetraethoxysilane TEOS) and hydrophobic products. Results confirm that undesirable changes may occur, with consequences ranging from purely aesthetic variations to physical, chemical and mechanical damages. This also permits us to check limitations in the use of these techniques for the evaluation of conservation treatments.

Fracture patterns and petrophysical properties of carbonates undergoing regional folding : A case study from Kurdistan, N Iraq

Acknowledgements The authors thank the Ministry of Natural Resources in Iraqi Kurdistan Region for permission to publish this paper. Gulf Keystone Petroleum Ltd. and HKN Energy Ltd. are acknowledged for providing the subsurface datasets. Great thanks to Colin Taylor at the University of Aberdeen for his assistance in the laboratory work. Thoughtful reviews by two anonymous referees improved the clarity of the paper. Graham Banks is thanked for his helpful and constructive review on a late version of the manuscript, which has significantly improved this paper.

Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low-shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional-wave velocity is strongly correlated with the shear-wave velocity in this suite of rocks. The ratio V-P/V-S varies significantly with change of both pore-fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot-flow and squirt-flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt-flow of fluid between the pores of various geometries. The contribution of the Biot-flow loss mechanism to the total loss is negligible. The compressional-wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low-shale, high-porosity sandstone reservoirs.

Durability of crystalline monumental stones in terms of their petrophysical characteristics

The durability of stone building materials is an issue of utmost importance in the field of monument conservation. In order to be able to preserve our built cultural heritage, the thorough knowledge of its constituent materials and the understanding of the processes that affect them are indispensable. The main objective of this research was to evaluate the durability of a special stone type, the crystalline stones, in correlation with their intrinsic characteristics, the petrophysical properties. The crystalline stones are differentiated from the cemented stones on the basis of textural features. Their most important specific property is the usually low, fissure-like porosity. Stone types of significant monumental importance, like the marble or granite belong to this group. The selected materials for this investigation, indeed, are a marble (Macael marble, Spain) and a granite (Silvestre Vilachán granite, Spain). In addition, an andesite (Szob andesite, Hungary) also of significant monumental importance was selected. This way a wide range of crystalline rocks is covered in terms of petrogenesis: stones of metamorphic, magmatic and volcanic origin, which can be of importance in terms of mineralogical, petrological or physical characteristics. After the detailed characterization of the petrophysical properties of the selected stones, their durability was assessed by means of artificial ageing. The applied ageing tests were: the salt crystallization, the frost resistance in pure water and in the presence of soluble salts, the salt mist and the action of SO2 in the presence of humidity. The research aimed at the understanding of the mechanisms of each weathering process and at finding the petrophysical properties most decisive in the degradation of these materials. Among the several weathering mechanisms, the most important ones were found to be the physical stress due to crystallization pressure of both salt and ice, the thermal fatigue due to cyclic temperature changes and the chemical reactions (mostly the acidic attack) between the mineral phases and the external fluids. The properties that fundamentally control the degradation processes, and thus the durability of stones were found to be: the mineralogical and chemical composition; the hydraulic properties especially the water uptake, the permeability and the drying; the void space structure, especially the void size and aperture size distribution and the connectivity of the porous space; and the thermal and mechanical properties. Because of the complexity of the processes and the high number of determining properties, no mechanisms or characteristics could be identified as typical for crystalline stones. The durability or alterability of each stone type must be assessed according to its properties and not according to the textural or petrophysical classification they belong to. Finally, a critical review of standardized methods is presented, based on which an attempt was made for recommendations of the most adequate methodology for the characterization and durability assessment of crystalline stones.

Petrophysical investigation on sediment core CRP-3 from the Ross Sea, Antarctica

A suite of petrophysical properties - velocity, resistivity, bulk density, porosity, and matrix density - was measured on 88 core plugs from the CRP-3 drillhole. Core-plug bulk densities were used to recalibrate both whole-core and downhole bulk density logs. Core-plug measurements of matrix density permit conversion of the whole-core and downhole bulk density logs to porosity. Both velocity and formation factor (a normalized measure of resistivity) are strongly correlated with porosity. The velocity/porosity pattern is similar to that for the lower part of CRP-2A and is consistent with the empirical relationship for sandstones. Core-plug and whole-core measurements of P-wave velocity at atmospheric pressure exhibit excellent agreement. Measurements of velocity as a function of pressure indicate a significantly higher velocity sensitivity to pressure than has been observed at CRP-1 and CRP-2A; rebound or presence of microcracks at CRP-3 may be responsible. The percentage difference between velocities at in situ pressures and atmospheric pressures increases downhole from 0% at the seafloor to 9% at the bottom. This pattern can be used to correct whole-core velocity data, measured at atmospheric pressure, to in situ velocities for depth-to-time conversion and associated comparison to the seismic profile across the drillsite

(Table 2) Petrophysical, petrological, and geochemical characteristics of ODP Hole 163-990A sediments

Ocean Drilling Program Hole 990A penetrated 131 m of subaerially emplaced Paleocene flood basalts on the Southeast Greenland margin with a recovery of 74%. Shipboard P-wave velocity (Vp), density, and magnetic susceptibility were measured with 2- to 15-cm intervals on the core. Individual flow units were divided into four zones based on the observed petrophysical characteristics. From the top, these are Zone I (<7 m thick with a Vp of ~2.5 km/s), Zone II (3-5 m thick with a strongly increasing Vp from 2.5 to 5.5 km/s), Zone III (up to 20 m thick with a Vp of ~5.5-6.0 km/s), and Zone IV (<2 m thick with a strongly decreasing Vp from 6.0 to 2.5 km/s). Eighteen samples were selected from three of the fourteen penetrated basalt units for geochemical, petrological, and petrophysical studies focusing on the altered, low-velocity upper lava Zones I and II. Zone I is strongly altered to >50% clay minerals (smectite) and iron hydroxides, and the petrophysical properties are primarily determined by the clay properties. Zone II is intermediately altered with 5%-20% clay minerals, where the petrophysical properties are a function of both the degree of alteration and porosity variations. Shipboard and shore-based measurements of the same samples show that storage permanently lowers the elastic moduli of basalt from Zones I to III. This is related to the presence of even small quantities of swelling clays. The data show that alteration processes are important in determining the overall seismic properties of flood basalt constructions. The degree and depth of alteration is dependent on the primary lava flow emplacement structures and environment. Thus, the interplay of primary emplacement and secondary alteration structures determine the elastic properties of basalt piles. Rock property theories for sand-clay systems are further used to model the physical property variations in these altered crystalline rocks.

Caracterização do reservatório Marlim por inversão acústica

Exploracionistas tem grande interesse em sistemas turbidíticos, pois em geral estes compõem prolíficos plays exploratórios. No entanto, estes potenciais reservatórios se encontram muitas vezes perto ou abaixo de resolução sísmica. Dessa forma, no processo de inversão, é importante integrar os dados sísmicos com qualidade de resolução obtida a partir do seu pré-condicionamento, para que o resultado final possua características detalhadas das camadas. O pré-condicionamento possibilita melhora na resolução dos dados sísmicos, através da atenuação dos ruídos aleatórios. Como objetivo final, foi realizada a inversão acústica em dados sísmicos post stack, migrados em tempo, a um sistema de turbiditos na Bacia de Campos. O principal objetivo da inversão sísmica é transformar o dado de reflexão em propriedades petrofísicas quantitativas. A inversão para impedância acústica é comumente utilizada para predição de porosidade. O fluxo de trabalho proposto foi dividido em cinco estágios principais: pré-condicionamento sísmico do dado 3D, correlação poço-sísmica, construção do modelo de baixa frequência, inversão do dado, e estimativa da porosidade. Os resultados mostraram que o cubo de impedância acústica invertido possui resolução muito superior quando comparado com o dado em amplitude sísmica, possibilitando melhor visualização das feições geológicas do Campo de Marlim. Além de suas limitações, como desconsiderando os efeitos das variações de fluido e variações litológicas complexas sobre a relação porosidade/impedância, o método fornece uma ferramenta confiável para exploração sísmica. Detalhes mais precisos das propriedades petrofísicas podem ser obtidos através de métodos de inversão mais sofisticados, a partir de dados pre stack.

Caracterização de reservatórios não convencionais Shale Gas na Formação Barreirinha Bacia do Amazonas

Os folhelhos pretos devonianos da Formação Barreirinha caracterizamse pela alta radioatividade na porção basal, grande extensão areal, espessura e profundidade de soterramento variável que vão de exposição na superfície até mais de 3000 m. Eles são as principais rochas geradoras do sistema petrolífero convencional da Bacia do Amazonas, e recentemente foram consideradas como promissores plays de gás não convencional. Folhelhos são geralmente caracterizados por uma matriz fechada, que faz com que sejam relativamente impermeáveis em relação ao fluxo de gás, a menos que ocorram fraturas, e dependendo das suas características geológicas e geoquímicas podem funcionar com um Sistema Petrolífero autossuficiente, atuando tanto como rocha fonte, quanto como reservatório de gás (reservatório Shale Gas). Assim, o gás natural termogênico ou biogênico gerado pode ser armazenado em folhelhos ricos em matéria orgânica na forma livre, adsorvida, ou em estado dissolvido. Em contraste com os sistemas petrolíferos convencionais, reservatórios Shale Gas, possuem mecanismos de aprisionamento e armazenamento únicos, sendo necessária a utilização de técnicas de avaliações específicas. No entanto, folhelhos prolíficos geralmente podem ser reconhecidos a partir de alguns parâmetros básicos: arquitetura geológica e sedimentar, propriedades geoquímicas e petrofísicas e composição mineralógica. Tendo em vista a carência de pesquisas de caráter descritivo, com cunho exploratório dos folhelhos geradores da Formação Barreirinha, esta dissertação tem como objetivo introduzir uma metodologia de identificação de intervalos de folhelho gerador com potencial para reservatório Shale Gas. Começando com uma investigação regional sobre o contexto geológico e sedimentar, seguido de uma avaliação abrangente enfocando as características geoquímicas, petrofísicas e litofácies dos folhelhos a partir da integração de parâmetros obtidos de perfis geofísicos de poço, análises geoquímicas e aplicação dos conceitos de Estratigrafia de Sequencia.

黄骅坳陷港西开发区河流相储层构型与剩余油分布模式研究

In the production tail of oilfield, water-cut is very high in thick channel sand oil reservoir, but recovery efficiency is relative low, and recoverable remaining oil reserves is more abundant, so these reserves is potential target of additional development. The remaining oil generally distributed with accumulation in certain areas, controlled by the reservoir architecture that mainly is the lateral accretion shale beddings in the point bar, so the study of reservoir architecture and the remaining oil distribution patterns controlled by architecture are very significant. In this paper, taking the Minghuazhen formation of Gangxi oilfield as a case, using the method of hierarchy analysis, pattern fitting and multidimensional interaction, the architecture of the meandering river reservoir is precisely anatomized, and the remaining oil distribution patterns controlled by the different hierarchy architecture are summarized, which will help to guide the additional development of oil fields. Not only is the study significant to the remaining oil forecasting, but also it is important for the theory development of reservoir geology. With the knowledge of sequence correlation and fluvial correlation model, taking many factors into account, such as combination of well and seismic data, hierarchical controlling, sedimentary facies restraint, performance verification and 3-D closure, an accurate sequence frame of the study area was established. On the basis of high-resolution stratigraphic correlation, single layer and oil sand body are correlated within this frame, and four architecture hierarchies, composite channel, single channels, point bars and lateral accretion sandbody are identified, The result indicates that Minghuazhen Formation of Gangxi oilfield are dominated by meandering river deposition, including two types of channel sandbodies, narrow band and wide band channel sandbody, and each of them has different characteristics of facies variation laterally. Based on the identification of composite channel, according to the spatial combination patterns and identified signs of single channel, combined with channel sandbody distribution and tracer material data, single channel sandbodies are identified. According to empirical formula, point-bar scales of the study area are predicted, and three identification signs are summarized, that is, positive rhythm in depositional sequence, the maximum thick sand and near close to the abandoned channel, and point bars are identified. On the basis of point bar recognition, quantitative architecture models inner point bar are ascertained, taking the lateral accretion sand body and lateral accretion shale beddings in single well as foundation, and quantitative architecture models inner point bar as guidance, and result of tracer material data as controlling, the the lateral accretion sand body and lateral accretion shale beddings are forecasted interwell, so inner architecture of point bar is anatomied. 3-D structural model, 3-D facies model and 3-D petrophysical properties models are set up, spatial distribution characteristics of sedimentary facies and petrophysical properties is reappeared. On the basis of reservoir architecture analysis and performance production data, remaining oil distribution patterns controlled by different hierarchy architecture units, stacked channel, single channel and inner architecture of point bar, are summarized, which will help to guide the additional development of oil fields.

东海陆架盆地复杂天然气藏地震勘探判别技术研究

Exploration study proves that East sea shelf basin embeds abundant hydrocarbon resources. However, the exploration knowledge of this area is very low. Many problems in exploration are encountered here. One of them is that the gas reservoir of this area, with rapid lateral variation, is deeply buried. Correlation of Impendence between sandstone, gas sand and shale is very poor. Another problem is that the S/N ratio of the seismic data is very low and multiples are relatively productive which seriously affect reservoir identification. Resolution of the seismic data reflected from 2500-3000 meter is rather low, which seriously affects the application of hydrocarbon direct identification (HDI) technology. This research established a fine geological & geophysical model based on drilling、well logging、geology&seismic data of East sea Lishui area. A Q value extraction method from seismic data is proposed. With this method, Q value inversion from VSP data and seismic data is performed to determine the subsurface absorption of this area. Then wave propagation and absorption rule are in control. Field acquisition design can be directed. And at the same time, with the optimization of source system, the performance of high resolution seismic acquisition layout system is enhanced. So the firm foundation is ensured for east sea gas reservoir exploration. For solving the multiple and amplitude preserving problems during the seismic data processing, wave equation pre-stack amplitude preservation migration and wave equation feedback iteratively multiple attenuation technologies are developed. Amplitude preservation migration technology can preserve the amplitude of imaging condition and wave-field extrapolation. Multiple removing technology is independent of seismic source wavelet and velocity model, which avoiding the weakness of Delft method. Aiming at the complicated formation condition of the gas reservoir in this area, with dissecting typical hydrocarbon reservoir, a series of pertinent advanced gas reservoir seismic identification technologies such as petrophysical properties analyzing and seismic modeling technology、pre-stack/post-stack joint elastic inversion, attribute extraction technology based on seismic non-stationary signal theory and formation absorption characteristic and so on are studied and developed. Integrated analysis of pre-stack/post-stack seismic data, reservoir information, rock physics and attribute information is performed. And finally, a suit of gas reservoir identification technology is built, according to the geological and geophysical characteristics of this area. With developed innovative technologies, practical application and intergrated interpretation appraisal researches are carried out in Lishui 36-1.The validity of these technologies is tested and verified. Also the hydrocarbon charging possibility and position of those three east sea gas exploration targets are clearly pointed out.

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

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.

高分辨率非线性岩性物性联合反演方法和应用

In the prediction of complex reservoir with high heterogeneities in lithologic and petrophysical properties, because of inexact data (e.g., information-overlapping, information-incomplete, and noise-contaminated) and ambiguous physical relationship, inversion results suffer from non-uniqueness, instability and uncertainty. Thus, the reservoir prediction technologies based on the linear assumptions are unsuited for these complex areas. Based on the limitations of conventional technologies, the thesis conducts a series of researches on various kernel problems such as inversions from band-limited seismic data, inversion resolution, inversion stability, and ambiguous physical relationship. The thesis combines deterministic, statistical and nonlinear theories of geophysics, and integrates geological information, rock physics, well data and seismic data to predict lithologic and petrophysical parameters. The joint inversion technology is suited for the areas with complex depositional environment and complex rock-physical relationship. Combining nonlinear multistage Robinson seismic convolution model with unconventional Caianiello neural network, the thesis implements the unification of the deterministic and statistical inversion. Through Robinson seismic convolution model and nonlinear self-affine transform, the deterministic inversion is implemented by establishing a deterministic relationship between seismic impedance and seismic responses. So, this can ensure inversion reliability. Furthermore, through multistage seismic wavelet (MSW)/seismic inverse wavelet (MSIW) and Caianiello neural network, the statistical inversion is implemented by establishing a statistical relationship between seismic impedance and seismic responses. Thus, this can ensure the anti-noise ability. In this thesis, direct and indirect inversion modes are alternately used to estimate and revise the impedance value. Direct inversion result is used as the initial value of indirect inversion and finally high-resolution impedance profile is achieved by indirect inversion. This largely enhances inversion precision. In the thesis, a nonlinear rock physics convolution model is adopted to establish a relationship between impedance and porosity/clay-content. Through multistage decomposition and bidirectional edge wavelet detection, it can depict more complex rock physical relationship. Moreover, it uses the Caianiello neural network to implement the combination of deterministic inversion, statistical inversion and nonlinear theory. Last, by combined applications of direct inversion based on vertical edge detection wavelet and indirect inversion based on lateral edge detection wavelet, it implements the integrative application of geological information, well data and seismic impedance for estimation of high-resolution petrophysical parameters (porosity/clay-content). These inversion results can be used to reservoir prediction and characterization. Multi-well constrains and separate-frequency inversion modes are adopted in the thesis. The analyses of these sections of lithologic and petrophysical properties show that the low-frequency sections reflect the macro structure of the strata, while the middle/high-frequency sections reflect the detailed structure of the strata. Therefore, the high-resolution sections can be used to recognize the boundary of sand body and to predict the hydrocarbon zones.

鄂尔多斯盆地子洲气田二叠系山西组储层特征与天然气富集规律研究

In this thesis, detailed studies on the sedimentology and petrophysical properties of reservoir rocks in the Shan#2 Member of Shanxi Formation, Zizhou gas field of Ordos Basin, are carried out, based on outcrop description, core description, wireline log interpretation and analysis of petrophysical properties. In the context of stratigraphic division scheme of the Upper Paleozoic in Ordos Basin, the Shan#2 Member is further divided into three subintervals: the Shan#23, Shan#22 and Shan#21, based on the marker beds，depositional cycles, wireline log patterns. Subaqueous deltaic-front distributary channels, distributed from the south to north, is identified，which is the main reservoir sand bodies for gases of Shan#2 Member at Zizhou gas field. Quartzose and lithic-quartzose sandstones, commonly with a high volume of cement, but a low volume of matrix, are the major reservoir rocks in the studied area. All sandstones have been evolved into the late diagenetic stage (referred to as diagenetic stage B) during the burial, experiencing compaction, cementation, replacement and dissolution, in which the compaction and cementation could have reduced the porosity, while dissolution could have improved the petrophysical properties. The pore types in the reservoirs are dominated by intergranular-solutional, intergranular-intercrystal and intercrystal-solutional porosity. According to the parameters and capillary pressure curves of test samples, five types of pore texture (I-V) are differentiated, in which types II and III pore textures displayed by low threshold pressure-wide pore throat and moderate threshold pressure-moderately wide pore throat, exist widespread. Sandstone reservoirs in the studied area are characterized by exceptionally low porosity and permeability, in which the petrophysical properties of those in Shan#23 horizon are relatively better. The petrophysical property of reservoirs was influenced both by the sedimentation and diagenesis. In general, the coarse quartzose sandstones deposited in subaqueous distributary channels show the best petrophysical property, which tends to be worse as the grain size decreases and lithic amount increases. Three types of gas reservoirs in Shan#23 horizon are classified according to petrophysical properties (porosity and permeability), which could have been influenced by the initial depositional facies, diagenesis and tectonics. On the basis of the study on the geological conditions of reservoirs in the area, it is concluded that sedimetary facies, diagenesis and tectonic actions can provide an important foundation for gas pool formation, which can also control the accumulation and distribution of gas reservoirs.

油藏动态模型和剩余油分布模式

This paper belong to national "973" technological project undertaken by Shengli Oilfield. Work area is composed of turbidite reservoir of S212 and delta reservoir of S283 of Sheng2 unit in Shengtuo Oilfield that has a 36 years water injection development history. Change of the macroscopic, microscopic and filterational parameters and its mechanism have been studied in the 4 water-cut stages i.e. the primary, moderate, high and supper-high stage by using multi-disciplinary theories and approaches, computer techniques and physical simulation comprehensively. Dynamic reservoir models to different water-cut stages have been established. The study of forming mechanism and distribution of residual oil revealed the main types and spatial distribution of residual oil in different water-cut stages and the distribution mode has also been built up. Macroscopic, microscopic and filterational parameters selecting principle, optimizing and selecting standard, matching standard and laws and related database of various dynamic parameters in different water-cut stages have been established, which laid good basis for revealing reservoir macroscopic, microscopic and filterational parameters' dynamic change and residual oil distribution. The study indicated that in general, the macroscopic, microscopic and filterational parameters will slowly increase and become better in both shallow turbidite and delta reservoirs with the increasing of water cut, but different reservoirs have their own characteristics and change laws. Parameters of I~2 unit, whose petrophysical properties are better, increase more quickly than 8~3, whose petrophysical properties are more unfavorable. The changes was relatively quickly in high water-cut stage, while relatively slowly from primary to moderate and from high to supper-high water-cut stage. This paper firstly put forward that reservoir macroscopic, microscopic and filterational parameters are controlled by dynamic geological function of reservoir fluid, which is considered the major reason of reservoir parameters' dynamic changes and residual oil formation and distribution during reservoir development. Physical simulation of filterational parameters verified that forming mechanism and distribution of residual oil in different water-cut stages are also controlled by dynamic geological function of reservoir fluid. The idea of fluid geological function during reservoir development developed the theory of development geology, and has important practical values. This paper firstly constructed dynamic geological and mathematical models and five modes of residual oil distribution in Shengtuo Oilfield, and achieved four-dimensional forecast of residual oil distribution in different watercut stages. Dynamic changes and mechanism of macroscopic, microscopic and fliterational parameters of reservoir and their change process have been revealed. Forecast of residual oil distribution has been achieved by computers. This paper established the related theories, approaches and techniques for residual oil study, characterization and in different water-cut stages, and realized dynamic forecast of residual oil. It gained remarkable economic benefit and social effect in guiding field development. These theories and techniques had important meaningfulness for residual oil prediction in the terrestrial faulted basins not only in Shengli Oilfield but also in the east of China. Furthermore, this study has developed the theory of development geology.