砂泥岩储层测井产能预测方法研究

Productivity prediction is very important in the exploration and development of oilfields. Using well log data to predict productivity is a front-line technology， which is key issue in petroleum exploration phase. The essential factors of productivity prediction is building practical models and correcting various causes to improve precision of prediction parameters. Any errors of parameters selections can affect the calculation of productivity prediction; therefore， how to improve research means and calculation accuracy is an important task of productivity prediction. Theory and case examples are deeply and comprehensively studied in the paper. Based on the theory of mud-filtrate invasion and experimental results， the damage of drilling， cementing， perforating，acidizing and fracturing were investigated. The damage depth was quantitatively evaluated by log data， based on this， the processing results of reservoir sensitivity were used to analysis quantitatively the damage of reservoir. The productivity prediction and reservoir damage were initiatively incorporated according to well logging， and the precision of productivity prediction was effectively improved. The method of NMR was explored to calculate the fluid viscosity on the basis of reservoir physical method， and the differences between the two methods were compared in the paper. From the theory fluid flow in porous media， various of theoretical models of production prediction were explored and several practical models were consided， such as productivity index method， improved productivity index method， improved Bearder method， SVM and so on. The characteristic and the application scope of these methods were studied. The inflow productivity and outflow productivity were incorporated and nodal analysis method was used to forecast wellhead yield， thus achieved scientifically production. On the applied background of conventional logging suite， the applying of special items or new logging method which is practical in the research area were studied， the logging suite was further optimized， and the precision of forecast was improved. On the basis of the modeling and the calculation of parameters， these methods were verified and analyzed， and the reconstruct principle was also built for block reservoir. The research block was processed by these methods and compared with testing data. Based on above the research， a technological system which is practical for shaly sand profiles in Shengli Oilfield was built. The system can reach commercialized degree，and satisfied the need of exploration and development of the oilfield.

油田高含水期剩余油分布规律研究

After systemic investigation of the techniques，route lines and mechanisms about the remaining oil，the dynamic migration and congregation behavior of the remaining oil are discussed on base of interaction between flowing and enriching of water and oil.After the micro-scope modeling of the fluid flow in porous media and the changes in petrol-physical properties of the flowing system, the characters of fluid fields and the dynamic distribution of oil are discussed, among which the preference-flowing is focused on. Based on the preference-flowing in porous media, the concept of the preference-flowing channels is developed. According to above, heterogeneous distribution of water and oil in the field and dynamic mechanism of remaining oil are all obvious. media can be divided into three kinds, directional, stochastic, arbitrary porous media. The main research results are as following: 1. Treating the characteristic parameters such as permeability, porosity and wettability as regional parameter, the fluid field with high water-cut has been established by geostatistical method, among which the difference of flowing pores and the changes of its petrol-physical properties during flooding are studied. 2. The flow process of water and oil are recurrent in physical simulation experiments, in which the mechanisms and phenomena are caught and analyzed. Fluid flow mechanics in porous media with preference-flowing channels have been studied. 3. The mutual coupling between water and oil is induced and the mathematical evolution equations including this interaction were built. . 4. Through coupling effect between flowing water and oil, the dynamic migration and congregation behavior of remaining oil depend upon this coupling. 5. Coupling between water and oil act as driving force and trapping force for the remaining oil. The coupling model of thesis has been verified by simplified the numerical model and compared results with Ng35 oil reservoir in Gudao oil field, it has important theoretical and application values for improving precision of remaining oil and production performance prediction, and is a new method for studying the mechanics of remaining oil in channeled porous media has been established. Key words：flow field，high water-cut，coupling，dominant flow in porous media，remaining oil

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

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

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

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.

地震勘探三维可视化方法研究与实现

This thesis mainly studies the technologies of 3-D seismic visualization and Graphic User Interface of seismic processing software. By studying Computer Graphics and 3-D geological modeling, the author designs and implements the visualization module of seismic data processing software using OpenGL and Motif. Setting seismic visualization flow as the subject, NURBS surface approximation and Delaunay Triangulation as the two different methods, the thesis discusses the key algorithms and technologies of seismic visualization and attempts to apply Octree Space Partitioning and Mip Mapping to enhance system performance. According to the research mentioned above, in view of portability and scalability, the author adopts Object-oriented Analysis and Object-oriented Design, uses standard C++ as programming language, OpenGL as 3-D graphics library and Motif as GUI developing tool to implement the seismic visualization framework on SGI Irix platform. This thesis also studies the solution of fluid equations in porous media. 2-D alternating direction implicit procedure has been turned into 3-D successive over relaxation iteration, which possesses such virtues as faster computing speed, faster convergence rate, better adaptability to heterogeneous media and less memory demanding.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.

孔隙介质中油气二次运移动力学机制研究

Study of dynamical mechanism of hydrocarbon secondary migration is the key research project of China Petroleum and Chemical Corporation in the ninth "Five-Year Plan", and this research is the hot and difficult spot and frontline in the domain of reservoir forming study in recent years. It is a common recognition that the dynamical mechanism of hydrocarbon's secondary migration is the important factor to control the reservoir type, distribution and oil/gas abundance. Therefore, to study this mechanism and establish the modes of hydrocarbon's migration and accumulation in different reservoirs under different conditions are of great theoretical meaningfulness and practical value on both developing the theory and method of hydrocarbon migration/accumulation dynamics in terrestrial rift-subsidence lacustrine basins and guiding the exploration and production. A laboratory for physical simulation of hydrocarbon's secondary migration/accumulation mechanism has been build up. 12 types of physical simulation tests to determine the volume of oil/gas migration and accumulation within these 3 series of plentiful hydrocarbon sources, different hydrocarbon abundance and pore level have been carried out under the guide of multidisciplinary theories, applying various methods and techniques, and 24 migration/accumulation modes have been established. The innovative results and recognition are as follows: 1, The oil/gas migration and accumulation modes for sandstones of moderate, fine grain size and silt in these six paleo depositional environments of shallow lake, fluvial, lacustrine, fluvial-deltaic, turbidite-delta, and salty-semi salty lake have been established. A new view has been put forward that the oil/gas volumetric increment during their migration and accumulation in different porous media of different rocks has similar features and evolution history. 2. During oil/gas migration and accumulation in different grain-sized porous media or different reservoirs, all the volumetric increment had experienced three period of rapid increasing, balanced and slower increasing and limited increasing. The dynamical process of oil/gas secondary migration and accumulation has been expounded. 3 The two new concepts of "source supply abundance" and "source supply intensity" have been proposed for the first time, and the physical simulation for hydrocarbon's migration, accumulation and forming a reservoir has been realized. 4, Source supply abundance is the important factor to control the accumulated volume of oil phase in the porous media. It is impossible to accumulate large amount of hydrocarbon volume in an open boundary system when the source supply abundance is low, i.e. impossible to form reservoirs of high productivity. 5 The above 12 types of physical simulation tests indicated that enough energy (pressure) of the oil sources is the decisive factor to ensure hydrocarbon's entering, flowing and accumulating through porous media, and both oil and gas phase will accumulate into the favorable places nearest to the oil sources. 6 The theory, method and related techniques for physical simulation of hydrocarbon's secondary migration/accumulation mechanism have been formed and applied to the E&P of Shengtuo rollover anticline and Niuzhuang turbidite lithological reservoirs. 7 This study developed the theory and method of hydrocarbon migration/accumulation dynamics in terrestrial rift-subsidence lacustrine basins, and the benefits and social effect are remarkable.

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

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

花岗岩风化土微观结构及其工程地质特性研究

Saprolite is the residual soil resulted from completely weathered or highly weathered granite and with corestones of parent rock. It is widely distributed in Hong Kong. Slope instability usually happens in this layer of residual soil and thus it is very important to study the engineering geological properties of Saprolite. Due to the relic granitic texture, the deformation and strength characteristics of Saprolite are very different from normal residual soils. In order to investigate the effects of the special microstructure on soil deformation and strength, a series of physical, chemical and mechanical tests were conducted on Saprolite at Kowloon, Hong Kong. The tests include chemical analysis, particle size analysis, mineral composition analysis, mercury injection, consolidation test, direct shear test, triaxial shear test, optical analysis, SEM & TEM analysis, and triaxial shear tests under real-time CT monitoring.Based on the testing results, intensity and degree of weathering were classified, factors affecting and controlling the deformation and strength of Saprolite were identified, and the interaction between those factors were analyzed.The major parameters describing soil microstructure were introduced mainly based on optical thin section analysis results. These parameters are of importance and physical meaning to describe particle shape, particle size distribution (PSD), and for numerical modeling of soil microstructure. A few parameters to depict particle geometry were proposed or improved. These parameters can be used to regenerate the particle shape and its distribution. Fractal dimension of particle shape was proposed to describe irregularity of particle shapes and capacity of space filling quantitatively. And the effect of fractal dimension of particle shape on soil strength was analyzed. At the same time, structural coefficient - a combined parameter which can quantify the overall microstructure of rock or soil was introduced to study Saprolite and the results are very positive. The study emphasized on the fractal characteristics of PSD and pore structure by applying fractal theory and method. With the results from thin section analysis and mercury injection, it was shown that at least two fractal dimensions Dfl(DB) and Df2 (Dw), exist for both PSD and pore structure. The reasons and physical meanings behind multi-fractal dimensions were analyzed. The fractal dimensions were used to calculate the formation depth and weathering rate of granite at Kowloon. As practical applications, correlations and mathematical models for fractal dimensions and engineering properties of soil were established. The correlation between fractal dimensions and mechanical properties of soil shows that the internal friction angle is mainly governed by Dfl 9 corresponding to coarse grain components, while the cohesion depends on Df2 , corresponding to fine grain components. The correlations between the fractal dimension, friction angle and cohesion are positive linear.Fractal models of PSD and pore size distribution were derived theoretically. Fragmentation mechanism of grains was also analyzed from the viewpoint of fractal. A simple function was derived to define the theoretical relationship between the water characteristic curve (WCC) and fractal dimension, based on a number of classical WCC models. This relationship provides a new analytical tool and research method for hydraulic properties in porous media and solute transportation. It also endues fractal dimensions with new physical meanings and facilitates applications of fractal dimensions in water retention characteristics, ground water movement, and environmental engineering.Based on the conclusions from the fractal characteristics of Saprolite, size effect on strength was expressed by fractal dimension. This function is in complete agreement with classical Weibull model and a simple function was derived to represent the relationship between them.In this thesis, the phenomenon of multi-fractal dimensions was theoretically analyzed and verified with WCC and saprolite PSD results, it was then concluded that multi-fractal can describe the characteristics of one object more accurately, compared to single fractal dimension. The multi-fractal of saprolite reflects its structural heterogeneity and changeable stress environment during the evolution history.

遵义市岩溶水地下水流和污染物运移数值模拟研究

地下水数值模拟技术已成为评估人类活动对地下水质和量的影响、评价地下水资源、预测地下水污染发展趋势等的最主要的方法和手段。喀斯特含水层由于含水介质和地下水流场的非均质性和各向异性，对其进行地下水流的数值模拟一直是水文地质学界的难题。 遵义市地处我国西南喀斯特发育区，为贵州省第二大工业城市，属重度缺水地区，地下水资源的开发利用极大缓解了区内的缺水危机。但长期以来，由于对地下水资源的开发利用缺乏合理的统筹规划和强有力的管理，引发了一些环境地质问题，如地下水降落漏斗、岩溶塌陷、地下水质恶化等。因此选择遵义市进行地下水流和污染物运移数值模拟研究具有理论和实际意义。 通过ArcGIS平台建立了研究区的水文地质信息数据库，对研究区地下水的水位动态以及水化学特征做了简单分析。概要总结和阐述了高桥－河溪坝块段的自然地理、地质概况和水文地质条件，建立了水文地质概念模型；在水文地质概念模型的基础上，利用Groundwater Vistas软件建立了枯水期和丰水期的二维非均质各向异性稳定流模型，三维有限差分地下水流模拟程序MODFLOW用于模拟地下水水流，三维溶质运移模块MT3DMS用于模拟污染物在对流弥散情况下的迁移。根据分析和模拟结果可以得出如下几点结论： 1、基岩裂隙水水位峰值滞后大气降水峰值2～3个月，属渐峰型动态；岩溶水水位、地下河出口和泉流量变化步调与降水强度一致，对降水响应敏感。 2、对NH4+、NO3-、NO2-、SO42-、Mn五种组分含量进行了时空分析，结果表明地下水污染物的含量可能受人为活动输入物质的不均匀性和降雨等各方面因素控制，各组分每年的污染面积不一致，没有明显的规律性；受污染的一般是岩溶水，尤其是在石灰岩溶洞、地下河强烈发育而三废排放量大的居民集中地区面积较大。 3、为了有效地进行地下水资源管理，论文对高桥－河溪坝岩溶含水系统进行了一定的概化，将岩溶含水介质近似作为等价多孔介质（Equivalent Porous Media, EPM）模型来进行研究，采用MODFLOW的六个子程序模拟含水层系统的源汇项：降水子程序包RCH模拟降水入渗量、井流子程序包WEL模拟抽水量、通用水头子程序包GHB模拟侧向补给/排泄量、排水沟渠子程序包DRN模拟地下河出口流量、河流子程序包RIV模拟河流与地下水的交换量和已知水头边界子程序包CHD。从水位观测点和地下水位等势面两者结合来校正模型，结果表明能够达到相应国家标准规定的要求。因此EPM模型是可以适用于我国西南喀斯特地区的地下水流模拟的。 4、通过稳定流模型识别了枯水期和丰水期的渗透系数。在高桥和茅草铺附近渗透系数较高，枯水期介于100～400 m/d，而丰水期在高桥最高可达到3220m/d；其余单元渗透系数低于100 m/d，大多数小于10m/d。总体来说，由于丰水期含水层的饱水度大，渗透系数要高于枯水期。 5、通过地下水均衡计算，确定了各补给项和排泄项的水量。枯水期最重要的补给来源是研究区东北角的侧向补给量，占总补给量的70％，人工开采是最大的排泄项；丰水期最重要的补给源是西部的已知水头边界，占总补给量的49％，东北角的侧向补给量是第二补给源，占39％，地下河出口是最主要的排泄方式，达到排泄总量的74％。 6、对水文地质参数和源汇项敏感度分析的结果表明，不管是枯水期还是丰水期，对研究区水位影响最大的是渗透系数，外部源汇项中则是抽水量对地下水流形态的影响最大。 7、研究区岩溶地下水流速很大，污染物的运移是一个对流占绝对优势的问题，弥散的作用则相对很小。通过在茅草铺地区假设污染源，用MT3DMS程序模拟了地下水污染物在时间和空间上的迁移特征。结果发现：污染羽的形状和扩散方位主要受地下水流场的控制，而污染物的浓度与水量多少相关。

Thermocapillary Convection In An Inhomogeneous Porous Layer

A new model consisting of an inhomogeneous porous medium saturated by incompressible fluid is investigated. We focus on the effects of inhomogeneity for the streamline patterns and instabilities of the system. Influences of the 'mean porosity' and gradient of distributions of porosity are also emphasized. The results cannot be obtained by studying the media with constant porosity as carried out by other researchers, and have not been discussed before.