A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

九寨沟湖泊湿地植物群落与植物地理研究

九寨沟湖泊湿地在维持九寨沟的生态平衡中起着重要的作用，在旅游产业的发展下，湿地生态系统及生物多样性面临着较大的威胁。尽管九寨沟湿地具有重要的生态价值，但目前对其研究尚比较薄弱。湿地植物群落和植物地理研究可以为湿地资源的可持续利用和监测保护提供科学依据。作者从2004年8月到2007年11月对九寨沟湿地的植物物种组成、地理分布、优势植物群落的结构、生长动态、湿地土壤种子库进行了调查研究。主要结果如下： 1. 九寨沟湿地物种组成、地理分布特点及湿地植物群落特点 九寨沟湿地共有苔藓植物8科13属16种，维管植物为48科107属199种。九寨沟湿地植物的地理成份较为丰富，维管植物在科级水平上有7种地理分布型（变型），在属级水平上有13种地理分布型（变型）， 在种级水平上共有29种地理分布型（变型）。九寨沟湿地植物以温带成份和我国特有成份为主，同时兼有热带、亚热带成份和环极—高山成份。九寨沟湿地植物的分布表现出明显的垂直地带性和水平地带性。湿地植物群落可划分21个群落类型，不同植物群落类型的物种多样性及物种组成存在较大的差异。九寨沟湿地植物的物种多样性和群落多样性以及较高的生产力特征，是维持其湿地生态景观多样性和稳定性的基础。 2. 土壤、水环境、海拔等对湿地植物的分布及生物多样性的影响 九寨沟湿地土壤、水等环境因子存在较大的差异。帕米尔苔草和宽叶香蒲等群落的凋落物较多，土壤有机碳、土壤总磷较高，可能是九寨沟湿地的重要土壤碳库。 九寨沟湿地植物沿水环境梯度的分布规律表现为：沉水植物（轮藻—篦齿眼子菜，水苦荬，杉叶藻）——挺水植物（水木贼，芦苇，宽叶香蒲）——湿生草本（苔草、节节草、披散木贼）——湿生灌木（柳灌丛，小檗灌丛）等。海拔也影响湿地植物的物种组成。 水深对物种多样性有影响，水深与物种丰富度负相关。随着水深的增加，水木贼、芦苇、杉叶藻、宽叶香蒲等群落的物种多样性下降；在长期淹水和季节性淹水的地方，水木贼群落物种多样性存在显著差异。土壤总氮与水木贼群落物种丰富度正相关。 3. 土壤营养元素、水环境对植物生长的影响 水深影响湿地植物生物量的分配。芦苇无性系分株在47 cm水深的环境中单株平均生物量最大；在干滩地中（地面水深0 cm），叶生物量百分比最大，而茎生物量百分比最小，茎的生物量百分比和生长速率随水深的增加而增加；在较干的滩地生境中，开花率、花序的生物量百分比明显大于水较深的生境。 水深与水木贼地上生物量负相关，但水木贼地上生物量在长期淹水和季节性淹水的地方没有显著的差异。在水浅的地方，杉叶藻、水木贼、芦苇等植物群落中，其他伴生物种的生物量占样方总生物量的百分比较大。 土壤有机碳、土壤总氮、土壤总磷等对湿地植物生物量的影响比较大：宽叶香蒲地上生物量与土壤总磷正相关；水木贼地上生物量与土壤总氮正相关；杉叶藻地上生物量与土壤有机碳正相关。 水深、土壤营养成分对湿地植物高度、密度等有影响。水木贼的平均高度在季节性淹水的地方比长期淹水的地方低，平均密度在长期淹水的地方比季节性淹水的地方低；除了5月份，其他观察月份水木贼的密度都与水深负相关，同时与土壤有机碳正相关。另外，芦苇密度与土壤有机碳含量正相关，宽叶香蒲密度与水深负相关，帕米尔苔草高度与土壤有机碳负相关。 4. 优势植物群落的动态变化 在优势植物群落中，优势种的高度、密度、盖度、生物量等在群落中占绝对优势。除五花海，水木贼群落的物种组成、高度、生物量在两年间没有显著的变化。芦苇群落的物种丰富度在近两年有所增加。 湿地植物生长表现为明显的季节动态，生长的峰值大多在7月－8月。优势植物群落的物候与水文周期有关。湿地植物群落的物种组成和密度，可以作为对湿地监测和保护的生物指示。 5. 九寨沟湿地土壤种子库特征及其在湿地生物多样性恢复中的作用 水深和现存植被物种丰富度可以解释湿地土壤种子库的变化。水深可以解释表层物种丰富度45%的变化。现存植被物种丰富度可以分别解释10 cm土层、2－5 cm土层及5－10 cm土层土壤种子库45%、48%和25%的变化。 湿地土壤种子库的密度为0－15945粒m-2, 种子库中共发现23个物种。现存植被优势物种和种子库优势物种不同。各层土壤种子库密度和物种丰富度并不存在显著的差异，但第二层土壤种子库密度最大。海拔、现存植被优势种盖度、土壤总磷、土壤总氮、土壤有机碳对湿地土壤种子库的密度和垂直结构没有影响。土壤种子库物种丰富度小于地上植被物种丰富度。湿地土壤种子库与地上植被的相关性不大。在浅水区域，湿地土壤种子库在湿地植被恢复中有一定作用。但在深水区域，保护现存植被更重要。 The lakeshore wetlands are valuable ecological units of the Jiuzhaigou lakes. Pressure for travel industry development pose a continuing and severe threat to the biodiversity-support function of the wetland system. Despite the ecological importance of wetlands in Jiuzhaigou, they are so far poorly studied. Both general plant communties and biogeographical studies are needed in order to attain basis for sustainable use the wetland resources and adequate protection of these areas. The present study was undertaken to examine aquatic plants distribution and the species compositon, structure and growth dynamics of their communities with variations of environmental factors along altitudes, water depth and soil properities gradients in Jiuzhaigou. Analysis of field survey data collected during August 2004 and November 2007 in lakeshore wetlands in Jiuzhaigou National Nature Reserve, Sichuan, China. The results were as following: (i) Species composition and biogeography in wetland vegetation 8 families, 13 genus, 16 species of moss and 48 families, 107 genus and 199 species of vascular plants in Jiuzhaigou wetlands were found. The floristic compositions were abundunt. Ten geographical distribution types at family level, 13 geographical distributions types at generic level and 29 geographical distribution types at specific level in vascular plants were found. Most species in Jiuzhaigou wetlands are temperate elements and Chinese endemic elements, with a few of tropical and subtropical and some circumarctic elements. And the plant distributions show clear vertical and horizontal patterns. There were 21 major wetland plant community types. Species composition and species richness in different plant communities are different. The species diversity and plant community diversity and their high biomass are the basis for the diversity and stability of wetland landscapes in Jiuzhaigou. (ii) Water depth, soil nutrients and altitudes influence on the species diversity and plant distribution. Total phosphorous and organic cabon in soil were higher in C. pamiernensis and T. latifolia communities, where are important cabon reservoirs in Jiuzhaigou wetlands. Along gradients of water depth, among populations of the dominant plant species present: submerged macrophytes (Chara vulgaris, Potagemonton pectinatus, Veronica anagalis-aquatica，Hippuris vulgaris), emergent macrophytes (Equisetum fluviatile, Phragamites australis, Typha latifolia), helophytes (Carex pamirensis )and shrubs (Salix sp., Berberis sp. ). Altitudes influence on the assemblage of plant communities. Water depth negatively correlated with species richness. Specie richness showed differences between permanently flooded sites and seasonally flooded sites in E. fluvatile communities. And total nitrogen in soil was negatively correlated with species richness in E. fluviatile communities. Altitudes show no significant influence on species richness, but in fact, through our analyses, they do have influence on the assemblage of wetland plants. (iii) Water depth, soil nutrients influence on the plant growth Water depth influences the biomass allocation in Phragmities australis. The average aboveground biomass of a single ramet (4.2 g) was the largest in the habitat with water level 47 cm above the soil surface. At the habitat with water level under soil surface 15 cm (－15 cm), the leaf biomass percentage (of the total ramet biomass) was the largest (46.1%), and the height and percentage of ramose ramets ( with branches on stem )(of the total ramets in a plot) were found obviously different. The deeper in water, the larger the biomass percentage and growth rate of stems were. The flowering rate and biomass of panicles were greater in shallow water than those in deep water. Water depth negatively correlated with aboveground biomass of E. fluviatile. However, above-ground biomass of E. fluviatile showed no significant difference between permanently flooded sites and seasonally flooded sites. But in shallow water, more biomasses of accompanying species were found in dominant plant communities such as H. vulgaris communities, E. fluviatile communities and P. australis communities. Water depth, soil nutrients influence on shoot density and shoot length of wetland plants. The shoot density of E. fluviatile was correlated to water depth in all growth months. Annual average density was significantly lower at permanently flooded sites than at seasonally flooded sites. But the annual average shoot length was significantly lower at seasonally flooded sites than at permanently flooded sites. (iv) Growth dynamics of dominant communities in Jiuzhaigou wetland The shoot length and shoot density, coverage and biomass of domiant species were dominated in plant communities. The species composition increased in P. australis communities in recent two years. The species richness in E. fluviatile communities showed no difference between 2005 and 2007. The above-ground biomass and shoot density in Five－flower Lake from July 2005 to July 2007 were significantly different, while in other sites, the differences were not significant. Shoot height, shoot density and above-ground biomass showed significant seasonal changes in all sites. Growth dynamics correlated with the cycle of water levels in lakes. Most plants growth parameters peaked at July or August. The biomass of T. latifolia peaked in August. But the shoot length of T. latifolia in deeper water peaked in July. The shoot length of E. fluviatile increased significantly from May to August except in seasonally flooded sites in Arrow-bamboo Lake. The species composition of communities and shoot density can be used as bioindicators in Jiuzhaigou wetland. (v) Soil seed bank in Jiuzhaigou wetland and its role in vegetation restoration Seed density in all soil layer samples was negatively correlated to water depth. Water depth can explain 45% variance of species richness in surface layer in sediment. Species richness in extant vegetation can explain 45%, 48%, 25% variance of species richness in total 10 cm and in 2－5 cm and 5－10 cm layer sediment respectively. Mean seed densities in wetlands ranged from 0 to 15945 m–2. A total of 23 species germinated in seed bank. The dominant species in seed bank and extant vegetation showed great difference. The total number of species and seedlings that germinated in different layers was not significantly different. But the second layer had the greatest seed density. In shallow water, seed bank can contribute to vegetation restoration, while in deeper water, protection of extant vegetation may be a better strategy.

Heavy-ion conformal irradiation in the shallow-seated tumor therapy terminal at HIRFL

Basic research related to heavy-ion cancer therapy has been done at the Institute of Modern Physics (IMP), Chinese Academy of Sciences since 1995. Now a plan of clinical trial with heavy ions has been launched at IMP. First, superficially placed tumor treatment with heavy ions is expected in the therapy terminal at the Heavy Ion Research Facility in Lanzhou (HIRFL), where carbon ion beams with energy up to 100 MeV/u can be supplied. The shallow-seated tumor therapy terminal at HIRFL is equipped with a passive beam delivery system including two orthogonal dipole magnets, which continuously scan pencil beams laterally and generate a broad and uniform irradiation field, a motor-driven energy degrader and a multi-leaf collimator. Two different types of range modulator, ripple filter and ridge filter with which Guassian-shaped physical dose and uniform biological effective dose Bragg peaks can be shaped for therapeutic ion beams respectively, have been designed and manufactured. Therefore, two-dimensional and three-dimensional conformal irradiations to tumors can be performed with the passive beam delivery system at the earlier therapy terminal. Both the conformal irradiation methods have been verified experimentally and carbon-ion conformal irradiations to patients with superficially placed tumors have been carried out at HIRFL since November 2006.

Physical property measurement of therapeutic carbon ion beam in the shallow-seated tumor therapy terminal at HIRFL

For the first time the physical properties of therapeutic carbon-ion beam supplied by, the shallow-seated tumor therapy terminal at the Heavy Ion Research Facility in Lanzhou (HIRFL) are measured. For a 80.55MeV/u C-12 ion beam delivered to the therapy terminal, the homogeneity of irradiation fields is 73.48%, when the beam intensity varied in the range of 0.001-0.1nA (i.e. 1 X 10(6) - 1 X 10(8) particles per second). The stability of the beam intensity within a few minutes is estimated to be 80.87%. The depth-dose distribution of the beam at the isocenter of the therapy facility is measured, and the position of the high-dose Bragg peak is found to be located at the water-equivalent depth of 13.866mm. Based on the relationship between beam energy and Bragg peak position, the corresponding beam energy at the isocenter of the therapy terminal is evaluated to be 71.71MeV/u for the original 80.55MeV/u C-12 ion beam, which consisted basically with calculation. The readout of the previously-used air-free ionization chamber regarding absorbed dose is calibrated as well in this experiment. The results indicate that the performance of the therapy facility should be optimized further to meet the requirements of clinical trial.

Shallow gas in the muddy sediments of Eckernförde Bay, Germany

The seafloor of central Eckernförde Bay is characterised by soft muddy sediments that contain free methane gas. Bubbles of free gas cause acoustic turbidity which is observed with acoustic remote sensing systems. Repeated surveys with subbottom profiler and side scan sonar revealed an annual period both of depth of the acoustic turbidity and backscatter strength. The effects are delayed by 3–4 months relative to the atmospheric temperature cycle. In addition, prominent pockmarks, partly related to gas seepage, were detected with the acoustic systems. In a direct approach gas concentrations were measured from cores using the gas chromatography technique. From different tests it is concluded that subsampling of a core should start at its base and should be completed as soon as possible, at least within 35 min after core recovery. Comparison of methane concentrations of summer and winter cores revealed no significant seasonal variation. Thus, it is concluded that the temperature and pressure influences upon solubility control the depth variability of acoustic turbidity which is observed with acoustic remote sensing systems. The delay relative to the atmospheric temperature cycle is caused by slow heat transfer through the water column. The atmospheric temperature cycle as ‘exiting function’ for variable gas solubility offers an opportunity for modelling and predicting the depth of the acoustic turbidity. In practice, however, small-scale variations of, e.g., salinity, or gas concentration profile in the sediment impose limits to predictions. In addition, oceanographic influences as mixing in the water column, variable water inflow, etc. are further complications that reduce the reliability of predictions.

Petroleum System in Deepwater Basins of the Northern South China Sea

The northern South China Sea margin has experienced a rifting stage and a post-rifting stage during the Cenozoic. In the rifting stage, the margin received lacustrine and shallow marine facies sediments. In the post-rifting thermal subsidence, the margin accumulated shallow marine facies and hemipelagic deposits, and the deepwater basins formed. Petroleum systems of deepwater setting have been imaged from seismic data and drill wells. Two kinds of source rocks including Paleogene lacustrine black shale and Oligocene-Early Miocene mudstone were developed in the deepwater basin of the South China Sea. The deepwater reservoirs are characterized by the deep sea channel rill, mass flow complexes and drowned reef carbonate platform. Profitable capping rocks on the top are mudstones with huge thickness in the post-rifting stage. Meanwhile, the faults developed during the rifting stage provide a migration path favourable for the formation of reservoirs. The analysis of seismic and drilling data suggests that the joint structural and stratigraphic traps could form giant hydrocarbon fields and hydrocarbon reservoirs including syn-rifting graben subaqueous delta, deepwater submarine fan sandstone and reef carbonate reservoirs.

Characteristics of shallow gas hydrate in Okhotsk Sea

Multidisciplinary field investigations were carried out in Okhotsk Sea by R/V Akademik M.A. Lavrentyev (LV) of the Russian Academy of Sciences (RAS) in May 2006, supported by funding agencies from Korea, Russia, Japan and China. Geophysical data including echo-sounder, bottom profile, side-scan-sonar, and gravity core sample were obtained aimed to understand the characteristics and formation mechanism of shallow gas hydrates. Based on the geophysical data, we found that the methane flare detected by echo-sounder was the evidence of free gas in the sediment, while the dome structure detected by side-scan sonar and bottom profile was the root of gas venting. Gas hydrate retrieved from core on top of the dome structure which was interbedded as thin lamination or lenses with thickness varying from a few millimeters to 3 cm. Gas hydrate content in hydrate-bearing intervals visually amounted to 5%-30% of the sediment volume. This paper argued that gases in the sediment core were not all from gas hydrate decomposition during the gravity core lifting process, free gases must existed in the gas hydrate stability zone, and tectonic structure like dome structure in this paper was free gas central, gas hydrate formed only when gases over-saturated in this gas central, away from these structures, gas hydrate could not form due to low gas concentration.

Polygonal faults and their implications for hydrocarbon reservoirs in the southern Qiongdongnan Basin, South China Sea

The Polygonal faults were identified in Qiongdongnan Basin, South China Sea, by using the technique of time coherent slice and horizon flattening of high-resolution 3D seismic data. These polygonal faults occur in three tiers of the upper Meishan Formation and the Huangliu Formation. The faults have lengths of 150-1500 m, spacings of 50-3000 m, throws of 10-40 m and dips of 50-90 degrees. Tectonic evolution in the Qiongdongnan Basin can be divided into a rifted stage and a post-rifted stage. Tectonic faults are widely distributed in the rifted sequences, but are not well developed in the post-rifted stage. Few faults in the post-rifted sequences might suggest the absence of a migration pathway for hydrocarbon or other fluids. However, the existence of polygonal faults in the post-rifted sequences can serve as the pathway and promote the hydrocarbon migration and accumulation in the Qiongdongnan Basin during the post-rifted stage. (C) 2010 Elsevier Ltd. All rights reserved.

Geophysical evaluation methods for buried hill reservoirs in the Jiyang superdepression of the Bohai Bay basin, eastern China

The Jiyang superdepression is one of the richest hydrocarbon accumulations in the Bohai Bay basin, eastern China. Comprehensive seismic methods have been used in buried hill exploration in Jiyang to describe these fractured reservoirs better. Accurate seismic stratigraphic demarcation and variable-velocity mapping were applied to reveal the inner structure of the buried hills and determine the nature of the structural traps more precisely. Based on the analysis of rock properties and the characteristics of well-developed buried hill reservoirs, we have successfully linked the geology and seismic response by applying seismic forward technology. Log-constrained inversion, absorption coefficient analysis and tectonic forward-inversion with FMI loggings were applied to analyse and evaluate the buried hill reservoirs and gave satisfying results. The reservoir prediction was successful, which confirmed that the comprehensive utilization of these methods can be helpful in the exploration of buried hill reservoirs.

Shallow water tidal constituents in the Bohai Sea and the Yellow Sea from a numerical adjoint model with TOPEX/POSEIDON altimeter data

A numerical adjoint model with TOPEX/POSEIDON (T/P) altimeter data was set up to investigate the shallow water tidal constituents in the Bohai Sea and the Yellow Sea. Shallow water tidal constituents W-4, MS4 and M-6) in the Bohai Sea and the Yellow Sea were first extracted from nearly 10 years of T/P data and then assimilated into a nonlinear barotropic tidal model by using adjoint method in order to fully describe the tides in this area. The general patterns of M-4 and MS4 solutions were in good agreement with those of Kang et al. (Cont. Shelf. Res. IS (1998) 739.) and Lefevre et al., (J. Geophys. Res. 105 (2000) 8707.). The RMS values for the principal constituents and coastal constituents were obviously less than those calculated by Kang et al. (1998) and Lefevre et al. (2000). It was shown that the calculated tidal constituents charts obtained in the present study were more accurate than those in other models. In the future the model will be applied to other coastal areas and some semi-enclosed seas. (C) 2004 Elsevier Ltd. All rights reserved.

The nonlinear evolution equation on the shallow water wave

Based on the variation principle, the nonlinear evolution model for the shallow water waves is established. The research shows the Duffing equation can be introduced to the evolution model of water wave with time.

A SINGULAR EXACT SOLUTION FOR NONLINEAR SHALLOW-WATER ON A ROTATING EARTH

In this paper, we present an exact solution for nonlinear shallow water on a rotating planet. It is a kind of solitary waves with always negative wave height and a celerity smaller than linear shallow water propagation speed square-root gh. In fact, it propagates with a speed equal to (1 + a/h) square-root gh(1 + a/h) where a is the negative wave height. The lowest point of the water surface is a singular point where the first order derivative has a discontinuity of the first kind. The horizontal scale of the wave has actually no connection with the water depth.

Intra-habitat heterogeneity of microbial food web structure under the regime of eutrophication and sediment resuspension in the large subtropical shallow Lake Taihu, China

Planktonic microbial community structure and classical food web were investigated in the large shallow eutrophic Lake Taihu (2338 km(2), mean depth 1.9 m) located in subtropical Southeast China. The water column of the lake was sampled biweekly at two sites located 22 km apart over a period of twelve month. Site 1 is under the regime of heavy eutrophication while Site 2 is governed by wind-driven sediment resuspension. Within-lake comparison indicates that phosphorus enrichment resulted in increased abundance of microbial components. However, the coupling between total phosphorus and abundance of microbial components was different between the two sites. Much stronger coupling was observed at Site 1 than at Site 2. The weak coupling at Site 2 was mainly caused by strong sediment resuspension, which limited growth of phytoplankton and, consequently, growth of bacterioplankton and other microbial components. High percentages of attached bacteria, which were strongly correlated with the biomass of phytoplankton, especially Microcystis spp., were found at Site 1 during summer and early autumn, but no such correlation was observed at Site 2. This potentially leads to differences in carbon flow through microbial food web at different locations. Overall, significant heterogeneity of microbial food web structure between the two sites was observed. Site-specific differences in nutrient enrichment (i.e. nitrogen and phosphorus) and sediment resuspension were identified as driving forces of the observed intra-habitat differences in food web structure.

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

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