Synthesis and ionic conductivity studies of solid polymer electrolytes based on modified alternating maleic anhydride copolymer with oligo(oxyethylene) side chains

Comb-like polymers (CP) based on modified alternating methyl vinyl ether/maleic anhydride copolymer with oligo-oxyethylene side chains of the type-O(CH2CH2O)(n)CH3 have been synthesized and characterized, and complexed with LiNO3 to form an amorphous polymer electrolyte. CP/salt complexes showed conductivity up to 10(-5) S/cm at room temperature. The temperature dependence of ionic conductivity suggests that the ion transport is controlled by segmental motion of the polymer, shown by linear curves obtained in Vogel-Tammann-Fulcher plots. The ionic conductivity maximum moves to a higher salt concentration as the temperature increases. IR results also indicate that the ester in CP might decompose at 140 degrees C and reproduce the maleic anhydride ring.

Lithium ion conducting polymer electrolytes based on alternating maleic anhydride copolymer with oligo-oxyethylene side chains

A comb polymer with oligo-oxyethylene side chains of the type -(CH2CB2O)(12)CH3 was prepared from methyl vinyl ether/maleic anhydride copolymer and poly (ethylene glycol) methyl ether. The polymer can dissolve LiClO4 salt to form homogeneous amorphous polymer electrolyte. The ac ion conduction was measured using the complex impedance method, and conductivities were investigated as functions of temperatures and salt concentration. The complexes were first found to have two classes of glass transition which increase with increasing salt content, The optimum conductivity attained at 25 degrees C is in the order of 5.50 x 10(-6)Scm(-1). IR spectroscopy was used to study the cation-polymer interaction.

A NEW SCHEME FOR ASSIGNMENT OF A CANONICAL CONNECTION TABLE

A new algorithm for deriving canonical numbering of atoms in a molecular graph has been developed. Some graph invariants, such as node properties, degree (connectivity), topological path, the smallest node ring index, etc., are encoded together to partit

Ship-board measurements and estimations of air-sea fluxes in the western tropical Pacific during TOGA COARE

Direct air-sea flux measurements were made on RN Kexue #1 at 40 degrees S, 156 degrees E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-alpha hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme, There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results, Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heal flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.

Wind-vector estimation for RADARSAT-1 SAR images: Validation of wind-direction estimates based upon geometry diversity

In this letter, a new wind-vector algorithm is presented that uses radar backscatter sigma(0) measurements at two adjacent subscenes of RADARSAT-1 synthetic aperture radar (SAR) images, with each subscene having slightly different geometry. Resultant wind vectors are validated using in situ buoy measurements and compared with wind vectors determined from a hybrid wind-retrieval model using wind directions determined by spectral analysis of wind-induced image streaks and observed by colocated QuikSCAT measurements. The hybrid wind-retrieval model consists of CMOD-IFR2 [applicable to C-band vertical-vertical (W) polarization] and a C-band copolarization ratio according to Kirchhoff scattering. The new algorithm displays improved skill in wind-vector estimation for RADARSAT-1 SAR data when compared to conventional wind-retrieval methodology. In addition, unlike conventional methods, the present method is applicable to RADARSAT-1 images both with and without visible streaks. However, this method requires ancillary data such as buoy measurements to resolve the ambiguity in retrieved wind direction.

大型激光叠焊系统研究与开发

介绍一种新型激光叠焊系统,利用该系统可将尺寸范围为12.0m×1.5m厚度不等的两块不锈钢板重叠焊接在一起。针对大尺寸不锈钢板尺寸大、易变形等特点提出了主—辅式压紧机构,有效地控制了焊接过程中大尺寸板材的变形并消除了板材间隙;设计了ufxyz型四轴联动方案:u轴与x轴同向,u轴驱动钢板进行大范围运动,x轴驱动激光头小范围运动。u轴与x轴的复合运动提高系统精度的同时提高了系统焊接动态响应速度。对实际焊接板材进行试验验证,结果表明使用该系统进行焊接的钢板可承受压强最大能达到6.1MPa。

地质导向钻井关键技术及方法

On the basis of analyzing the principle and realization of geo-steering drilling system, the key technologies and methods in it are systematically studied in this paper. In order to recognize lithology, distinguish stratum and track reservoirs, the techniques of MWD and data process about natural gamma, resistivity, inductive density and porosity are researched. The methods for pre-processing and standardizing MWD data and for converting geological data in directional and horizontal drilling are discussed, consequently the methods of data conversion between MD and TVD and those of formation description and adjacent well contrast are proposed. Researching the method of identifying sub-layer yields the techniques of single well explanation, multi-well evaluation and oil reservoir description. Using the extremum and variance clustering analysis realizes logging phase analysis and stratum subdivision and explanation, which provides a theoretical method and lays a technical basis for tracing oil reservoirs and achieving geo-steering drilling. Researching the technique for exploring the reservoir top with a holdup section provides a planning method of wellpath control scheme to trace oil and gas reservoir dynamically, which solves the problem of how to control well trajectory on condition that the layer’s TVD is uncertain. The control scheme and planning method of well path for meeting the demands of target hitting, soft landing and continuous steering respectively provide the technological guarantee to land safely and drill successfully for horizontal, extended-reach and multi-target wells. The integrative design and control technologies are researched based on geology, reservoir and drilling considering reservoir disclosing ratio as a primary index, and the methods for planning and control optimum wellpath under multi-target restriction, thus which lets the target wellpath lie the favorite position in oil reservoir during the process of geo-steering drilling. The BHA (bottomhole assembly) mechanical model is discussed using the finite element method, and the BHA design methods are given on the basis of mechanical analyses according to the shape of well trajectory and the characteristics of BHA’s structure and deformation. The methods for predicting the deflection rate of bent housing motors and designing their assemblies are proposed based on the principle of minimum potential energy, which can clearly show the relation between the BHA’s structure parameters and deflection rate, especially the key factors’ effect to the deflection rate. Moreover, the interaction model between bit and formation is discussed through the process of equivalent formation and equivalent bit considering the formation anisotropy and bit anisotropy on the basis of analyzing the influence factors of well trajectory. Accordingly, the inherence relationship among well trajectory, formation, bit and drilling direction is revealed, which lays the theory basis and technique for predicting and controlling well trajectory.

VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.

MB倒转在中国黄土中的记录及海陆对比

Chinese loess preserved in northwest and north China are famous for its fine grain size, high accumulate rate and high community and can be good archives for paleoclimate and paleomagnetic variation over the later Cenozoic, with which can be correlated well between marine sediments. Major geomagnetic chrons and long term paleoclimate changes in Quaternary are successfully extracted from Chinese loess-paleosols, as well as short-term geomagnetic excursions and climate instability of high resolution. Magneticstratigraphy based on paleogeomagnetic polarity reversal recorded in Chinese loess is a basic project in loess research since decades ago. True geomagnetic records and exact location of geomagnetic reversal boundary in section is the foundation of magneticstratigraphy. Matuyama-Brunhes (MB) reversal as the youngest one still remains divarication about exact location of its boundary (MBB). L8 and S8 of Luochuan and Xifeng located in the interior of Chinese Loess Plateau(CLP) and Mangshan in southeast part of CLP are chosen to make clear some problems which include magnetic mineral, process of MB reversal, location of MBB, downward displacement scale of magnetic reversal boundary, time lag of paleoclimate record by marine and loess, new correlation between Chinese loess and marine sediments. Rock-magnetic investigations carried on L8 and S8 show that the main mineral are ferrimagnetic assemblage consists of magnetite and maghemite in Luochuan and Xifeng, and magnetite in Mangshan, which all contains little hematite belongs to antiferromagnetic phases. The main carrier of nature remanet magnetism (NRM) is detrital magnetite with pseudo-single domain. Detailed paleomagnetic investigations display that there are several rapid reversals in direction during the process of MB reversal which started at the upper part of S8 and finished at the lower part of L8, and lasted about decades of centimeters to more than 100cm correspond to about 104 years. On the assumption that MBB is located in the middle part of the layer which recorded the very reversal, 11cm is considered as the scale of downward displacement for the MBB in Chinese loess after estimation through correlation between Luochuan and Mangshan records. So this study denies the theory of large scale displacement of MBB and large scale Lock-in depth of NRM acquired from Chinese loess. Time lag of paleoclimate records in terrestrial sediments and marine sediments is considered after reassessment of correlation between low field susceptibility of Chinese loess with marine oxygen isotope in benthic foraminifera. On the basis of traditional correlation between Chinese loess and marine oxygen isotope, this study document a new scheme which correlates L8 and S8 to MIS18 and MIS19, respectively.

松辽盆地南部腰英台区块低渗透储层整体压裂改造技术研究

The Yaoyingtai Block is located within the northeastern Changling Depression of southern Songliao Basin, where the reservoir sandstones are petrophysically characterized by very low permeability, which results in the low success probability of artificial fracturing, and the low oil yield by water injection in the course of oil production. In order to improve the situations as stated above, this research aims to work out an integral fracturing technology and strategy applicable to the low permeable reservoirs in Yaoyingtai Block. Under the guidance of geological theory, reservoir engineering and technology, the subsurface occurrences of natural and hydraulic fractures in the reservoirs are expected to be delineated, and appropriate fracturing fluids and proppants are to be optimized, based on the data of drilling, well logging, laboratory and field experiments, and geological data. These approaches lay the basis of the integral fracturing technology suitable for the low permeable reservoir in the study area. Based on core sample test, in-situ stress analysis of well logging, and forward and inversion stress field modeling, as well as fluid dynamic analysis, the maximum in-situ stress field is unraveled to be extended nearly along the E-W direction (clustering along N85-135°E) as is demonstrated by the E-W trending tensional fractures. Hydraulic fractures are distributed approximately along the E-W direction as well. Faulting activities could have exerted obvious influences on the distribution of fractures, which were preferentially developed along fault zones. Based on reservoir sensitivity analysis, integrated with studies on rock mechanics, in-situ stress, natural fracture distribution and production in injection-production pilot area, the influences of primary fractures on fracturing operation are analyzed, and a diagnostic technology for primary fractures during depressurization is accordingly developed. An appropriate fracturing fluid (hydroxypropyl guar gum) and a proppant (Yixing ceramsite, with a moderate-density, 0.45-0.9mm in size) applicable to Qingshankou Formation reservoir are worked out through extensive optimization analysis. The fracturing fluid can decrease the damage to the oil reservoir, and the friction in fracturing operation, improving the effect of fracturing operation. Some problems, such as sand-out at early stage and low success rate of fracturing operations, have been effectively solved, through pre-fracturing formation evaluation, “suspension plug” fracturing, real-time monitoring and limited-flow fracturing. Through analysis of fracture-bearing tight reservoir with variable densities and dynamic analysis of influences of well patterns on fracturing by using numerical simulation, a fracturing operation scheme for the Qingshankou Formation reservoir is proposed here as being better to compress the short factures, rather than to compress the long fractures during hydraulic fracturing. It is suggested to adopt the 450m×150m inverted 9-spot well pattern in a diamond shape with wells placed parallel to fractures and a half fracture length of 60-75m.

白垩纪超静磁带地球磁场古强度测定

The long-term variations in the strength of the geomagnetic field provide important constraints on the chemical-physical processes of the Earth’s interior. Especially, the intensity of the geomagnetic field during the Cretaceous normal superchron (CNS) is crucial to understand the geodynamo. But a paucity of paleointensity further limits to obtain essential knowledge interior process in the deep earth. In order to improve the experimental efficiency, this study tried to apply two new rock magnetic methods (FORC diagram and low-temperature demagnetization technique) to determine the paleointensity. First, some problems in the theory and technique in paleointensity experiments were introduced. A combined palaeomagnetic and geochronologic study was further conducted on a basaltic lava sequence at Jianchang in Liaoning Province, northeastern China. Radiometric 40Ar/39Ar dating indicates that the volcanism occurred at about 119 Ma within the marine anomaly C34n in Cretaceous normal superchron (CNS). Rock magnetic investigations show that pseudo-single domain (PSD) titanium-poor titanomagnetite is dominant in the studied lava flows. Both stepwise thermal and alternating field demagnetizations isolate the well-defined normal characteristic remanent magnetization (ChRM) in three independent lava flows with a mean direction of D/I = 6.0/51.9 degree(a95 = 12.3degree). Palaeointensity was determined using the modified Thellier method with systematic partial thermoremanent magnetization (pTRM) checks on total 72 samples, but only 10 samples exhibit ideal linear behavior on the Arai-plot in the temperature interval of 300-560 C and yield an average paleointensity of (25.8+/-1.4)uT. In addition, slopes of the line defined by the initial and the final points on the Arai-plot for the other 18 samples with characteristic PSD features give an average paleointensity estimation value of (24.8+/-1.9)uT. The consistency of these two approaches confidently demonstrates the fidelity of our results. The overall mean field strength determined using both approaches are thus estimated to be (25.2+/-0.7 )uT. This value corresponds to the virtual dipole magnetic moments (VDM) of (4.5+/-0.1)E22 Am^2, which is about half of the value of present field. This finding suggests that palaeointensity just at the onset of the CNS is characterized by a weak magnetic field.

鄂尔多斯盆地志丹地区延长组储层综合评价

Sedimentary provenance direction，sedimentary facies，reservoir geological characteristic，pore structure； physical property characteristic，reservoir classification and evaluation ，forthermore，favorable area area are forecasted of Yanchang formation in ZhiDan region， by mainly study on the data of field section observation ，core observation, well logging explaination and routine microscope slice identification，scanning Electron Microscope and reservoir analysis of lithology and physical property , Under the guidance of such advanced theories and methods as sedimentology，reservoir sedimentology，lithological oil pool and so on，in the thesis. The following fruits are mainly achieved in this paper： Yanchang formation stratum is divided and correlated in this entire region， and the characteristic of oil layer unit is detailed discussed , respectively. According to main marker bed and supplementary ones.and research result shows that the source of provenance direction during Yanchang Formation mianly is north-east. Delta and lake are mainly developed in study area ，sub-facies and micro-facies are divided，distribution of sedimentary micro-facies in plane and palaeogeographic evolution are described，and gentle slope type- shallow water delta depositional model is established. Fine-grain arkose sandstone is the main reservoir，and which have experienced such different degree diagenesis as compaction, cementation, replacement and dissolution, and in which compaction and cementation are mainly factors to reduce sandstone physical property and dissolution effectively improved physical property during burial diagenesis procedure. All reservoirs of Yanchang Formation have entered A period of late diagenetic stage according to scheme of diagenesis period division . Intergranular porosity，dissolution porosity，fissure porosity are main pore types. And porosity structure are analyse by mercury penetration capillary pressure parameter，fine-shortness type and fine- length throat type are mainly developed. as a whole，the reservoir， with the characteristic of porosity and permeability altering apparently，strong inhomogeneity ， is a medium- porosity and medium permeability one. In plane，higher- porosity and higher-permeability are corresponded well with distributary channel area， physical property and inhomogeneity are affected by both deposition and diagenesis，and distributary channel and underwater distributary channel are favorable facies . According to such characteristic as lithology，physical property，pore structure ，diagenesis and sandstone distribution, the sandy reservoir can be classified 4 types, and the main sandy in every oil layer unit are evaluated according to the standard. The analysis result of petroleum concentration rule shows that Yanchang Formation are with not only favourable oil source rock，reservoir，covering combination ，but also good entrapment condition in study area. Lithology and structure-lithology oil pool are mainly developed ，based on condition of favorable reservoir developments，accounting for deliverability and sandstone superface elevation，zone of profitabilitis are forecasted.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

河流成因储层沉积微相与剩余油分布及形成机理研究

Based on the study of fluvial sandstone reservoir in upper of Guantao group in Gudao and Gudong oilfields, this paper first introduces A.D.Miall's(1996a) architectural-element analysis method that was summarized from ground outcrop scale into the reservoir formation research of the study area, more subtly divides sedimentary microfacies and establishes sedimentary model of research area.on this base, this paper summarizes the laws of residual oil distribution of fluvial formation and the control effect of sedimentary microfacies to residual oil distribution, and reveals residual oil formation mechanism. These results have been applied to residual oil production, and the economic effect is good. This paper will be useful for residual oil research and production and enhancement of oil recovery in similar reservoir. The major conclusions of this paper are as follows. 1. Using the architectural-element analysis method to the core data, a interfacial division scheme of the first to the dixth scale is established for the studied fluvial formation. 2.Seven architectural-elements are divided in upper of Guantao group of study area. The sandstone group 5～1+2 of Neogene upper Gutao group belongs to high sinuous fine grain meandering river, and the sandstone group 6 is sandy braided river. 3. Inter layer, the residual oil saturation of "non-main layer" is higher than "main layer", but the residual recoverable reserve of former is larger. Therefore, "main layer" is the main body of residual oil distribution. The upper and middle part of inner layer has lower permeability and strong seeping resistance. Addition to gravity effect in process of driving, its driving efficiency is low; residual oil saturation is high. Because of controlling of inside non-permeable interlayer or sedimentary construction, the residual oil saturation of non-driving or lower driving efficiency position also is high. On plane, the position of high residual oil saturation mostly is at element LV, CS, CH (FF), FF etc, Which has lower porosity and permeability, as well as lens sand-body and sand-body edge that is not controlled by well-net, non-perfect area of injection and production, lower press difference resort area of inter-well diffiuent-line and shelter from fault, local high position of small structure. 4.Microscopic residual oil mainly includes the non-moved oil in the structure of fine pore network, oil in fine pore and path, oil segment in pore and path vertical to flow direction, oil spot or oil film in big pore, residual oil in non-connective pore. 5.The most essential and internal controlling factor of fluvial formation residual oil distribution is sedimentary microfacies. Status of injection and production is the exterior controlling factor of residual oil distribution. 6. The controlling effect of formation sedimentary microfacies to residual oil distribution indicates inter-layer vertical sedimentary facies change in scale of injection and production layer-series, planar sedimentary face change and inner-layer vertical sedimentary rhythm and interbed in single layer to residual oil distribution. 7. It is difficult to clear up the inter-layer difference in scale of injection and production layer-series. The using status of minor layer is not good and its residual oil saturation is high relatively. It is obvious that inter-layer vertical sedimentary facies changes control inter-layer residual oil distribution at the same or similar conditions of injection and production. For fluvial formation, this vertical sedimentary facies change mainly is positive gyration. Namely, from down to top, channel sediment (element CHL, LA) changes into over-bank sediment (element LV, CR, CS). 8. In water-injection developing process of transverse connecting fluvial sandstone oil formation, injection water always comes into channel nearby, and breaks through along channel and orientation of high pressure gradient, does not expand into side of channel until pressure gradient of channel orientation changes into low. It brings about that water-driving status of over-bank sedimentary element formation (LV, CR, CS) is not good, residual oil saturation is high. In non-connective abandoned channel element (CH) formation with channel, because this reverse is difficult to control by injection and production well-series, its using status is not good, even terribly not good, residual oil is enrichment. 9. The rhythm and sedimentary structure, sedimentary facies change in single sand body brings about vertical changes of formation character, growth character of inner layer interbed. These are important factor of controlling and affecting vertical water spread volume and inner layer residual oil forming and distribution in single sand body. Positive rhythm, is the principal part of fluvial sandstone inner layer sedimentary rhythm. Namely, from down to upside, rock grain granularity changes from coarse to fine, seeping ability changes from strong to feebleness. It brings about that water-driving status of inner layer upside is not good, residual oil saturation is high. Inner layer interbed has different degree affecting and controlling effect to seeping of oil and water. Its affecting degree lies on interbed thickness, extending scale, position, and jeted segment of production or injection well. The effect of interbed at upside of oil formation to oil and water seeping is less; the effect of interbed at middle of oil formation to oil and water seeping is more. 10. Indoor experiment and research indicate that wettability, permeability step, vertical permeability, position of Kmax and ratio of oil viscousity and water viscousity all have great effect on the water-driving recovery ratio. 11. Microscopic residual oil distribution is affected and controlled by formation pore network structure, pressure field distribution, and oil characteristic. 12.The residual oil forming mechanism: the over-bank sedimentary element and upper part of a positive rhythm sandstone have fine pore and throat network, permeability is low, displacement pressure of pore and throat is high. The water-driving power usually falls short of displacement pressure that brings about injection water does not spread into these pore and throat network, thereby immovable oil area, namely residual oil, is formed. At underside of channel sedimentary element and positive rhythm sandstone, porosity and permeability is relatively high, connecting degree of pore and throat is high, displacement pressure of pore and throat is low. Thereby injection water is easy to enter into pore and throat, driving oil in them. Because the pore space is irregular, the surface of pore wall is coarse and non-flat. That the oil locate on concave hole of pore wall and the dead angle of pore, and the oil attaches on surface of pore wall by surface tension, are difficult to be peeled off, becoming water-driving residual oil (remaining oil). On the other hand, Because flowing section lessens, flowing resistance increase, action of capillary fore, or seeping speed decreases at process of transfer at pass narrow throat path in the course carried by driving water. The "oil drop", "oil bead", or "oil segment" peeled off by driving water is difficult to carry and to drive out by water at less pressure difference. Thereby they are enclosed in pore to form discontinuous residual oil. 13.This results described above have been applied in nine develop blocks of Gudao and Gudong oilfield. Its applying effect is marked through local injection production adjustment, deploying replacement well, repair hole, replacement envelop, block off water and profile control etc. Relative method and technology can be applied to other oil production area of Shengli oilfield, and obtain better economic and societal effect.

青藏高原壳幔形变与层圈耦合数值模拟

Based on multi-principle (such as structures, tectonics and kinematics) exploratory data and related results of continental dynamics in the Tibetan plateau, the author reconstructed the geological-geophysical model of lithospherical structure and tectonic deformation, and the kinetics boundary conditions for the model. Then, the author used the numerical scheme of Fast Lagrangian Analysis of Continua (FLAC), to stimulate the possible process of the stress field and deformational field in the Tibetan plateau and its adjacent area, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. With the above-mentioned results, the author discussed the relationship between crustal movement in shallow layer and the deformational process in interior layers, and its possible dynamic constraints in deep. At the end of the paper, an integrative model has been put forward to explain the outline images of crust-mantle deformation and coupling in the Tibetan Plateau. (1) The characteristics of crust-mantle structure of the Tibetan plateau have been shown to be very complex, and vertical and horizontal difference is significant. The general characteristics of crust-mantle of the Tibetan plateau may be that it's layering in depth direction, and shows blocking from south to north and belting from east to west, mainly according to the results of about 20 seismic sections, such as wide-angle seismic profiles, CMP, seismic tomography and so on. (2) The crust had shortened about 2200km, while the shortening is different for different block from south to north in the Tibetan plateau. It is about 11.5mm/a in Himalayan block, about 9.0mm/a in Lhas-Gangdese block, about 7.0mm/a in Qiangtang block and Songpan-Ganzi-Kekexili block, about 8.0mm/a in Kunlun-Qaidam, and about ll.Omm/a in Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. Which - in demonstrates the shortening rate decreases from south to north, but this rate increases near the north edge of the Tibetan plateau. The crust thickening rate is about 0.4mm/a in the whole Tibetan plateau; and this rate is about 0.5mm/a in Himalayan block, about 0.4mm/a in Lhas-Gangdese block, about 0.3mm/a in Qiangtang block, about 0.2mm/a in Songpan-Ganzi-Kekexili block and about O.lmm/a in Kunlun-Qaidam-Qilian block, since the convergence-collision between the Indian continent and Eurasia continent about 50Ma ago. This implies that the thickening rate decreases in the blocks of the Tibetan plateau. From south to north, the displacement of eastern boundary in the Tibetan plateau is about 37mm/a in Himalayan block, about 45mm/a in Lhas-Gangdese block, about 47mm/a in Qiangtang block, about 43mm/a in Songpan-Ganzi-Kekexili block, and about 35mm/a in Kunlun-Qaidam-Qilian block, since the collision-matching between the Indian continent and Eurasia continent had happened about 50Ma ago. This implies that the rate of eastward displacement is biggest in the middle of plateau, and decreases to both sides. The transition of S-N compression stress field in Tibetan Plateau, since about 28Ma+ ago, may be caused by two reasons: On one hand, the movement direction of Eurasia continent changed from northward to southward about 28Ma± ago in the northern plateau. On the other hand, the front belt that is located between India continent's and Eurasia continent's convergence-collision, had moved southward to high Himalayan from Indus-Brahmaputra suture almost at the same time in southern plateau. Affected by the stress field, the earlier tectonics rotated clockwise, NE and NW conjugate strike-slip faults developed, and the SN rift formed. This indicated that the EW movement started. The ratio between upper crust and lower crust of different blocks from south to north in the Tibetan plateau during the process of deformation are as following: about 3.5~5:1 in Himalayan block, about 1~5: 3-4 (which is about 1:3o--4 in south and about 4~5:3 in north) in Lhas-Gangdese block, about 1:3~447mm/a in these blocks: Which is located to the north of Banggong-nujiang suture.