任务分配问题的建模与求解

建立了极大极小任务分配问题的混合整数线性规划模型,提出一种矩阵作业解答,并与穷举解及混合整数线性规划解的计算复杂度进行了比较.理论分析和数值试验表明矩阵作业法对两类任务分配问题,极大极小和总体极小任务分配问题,有效地提供最优解.

智能机器人多传感器实验平台的图像伺服系统

结合智能机器人多传感器实验平台的研制工作 ,系统介绍了其图像伺服子系统的工作原理、系统构成以及实现方案 ,涉及图像处理、机器视觉和机器人跟踪控制等相关内容 ,并讨论了物体自动识别和抓取过程中需要注意的一些问题 .

基于模型的立体视觉的优化方法

针对 3D视觉检测中最典型的误差源图像中像素的量化误差 ,在已知模型的条件下 ,对双目立体视觉提出一种模型。利用该优化模型及改进的约束最小二乘法对未知模型的立体视觉解法做一些改进 ,可以明显减小系统的误差。仿真结果表明 ,与未知模型的立体视觉方法相比 ,在同样的量化误差条件下 ,该方法可将系统误差减小 50 %以上。

空值环境下关系数据库的更新 Ⅰ：扩展关系模型及基本运算

本文从空值语义及更新操作的关系出发，提出了一种新的扩展关系模型，用以组织更新操作下的含有空值的关系数据库中的信息．同时，定义了这种模型下的基本关系代数运算．为实现空值环境下关系数据库的数据更新奠定了基础。

CIMS环境下CAPP系统的设计

本文结合国家CIMS重点应用工厂——沈阳鼓风机厂实施 CIMS 的实际,对 CIMS 环境下的CAPP 系统的结构、采用的设计方法、涉及的关键技术,以及所需的基础工作进行了讨论.

Bézier曲线间最近距离的计算方法

针对Bzier曲线间最近距离计算问题,提出一种简捷、可靠的计算方法.该方法以Bernstein多项式算术运算为工具,建立Bzier曲线间最近距离的计算模型;然后充分利用Bzier曲面的凸包性质和de Casteljau分割算法进行求解.该方法几何意义明确,能有效地避免迭代初始值的选择和非线性方程组的求解,并可进一步推广应用于计算Bzier曲线/曲面间的最近距离.实验结果表明,该方法简捷、可靠且容易实现,与Newton-Raphson方法的融合可进一步提高该方法的运行速度.

改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

水热条件下钾长石碱性分解反应机理的实验研究

It is well known that our country is short of water-soluble potassium, but rich in insoluble potassium ores. Based on the work of the formers, using the orthogonal and monofactor experiments, the author optimized the production technology of micro-porous potassium silicon calcium mineral fertilizer by non-stirring hydrothermal chemical reaction when the alkaline earth booster CaO was available. The influences of temperature、time、reactant ratio and water-solid ratio on the dissolution rate of production’s elements were studied by orthogonal experiments, and the production technology was further optimized by monofactor experiments. By XRD、SEM、EDS and dissolving experiments, it was systematically studied that the effects of the reactant ratio、reaction time and reaction temperature on the properties of the production obtained by the hydrothermal reaction between KAlSi3O8 and CaO. The results showed that：when changing of the reaction condition, the reaction productions included tobermorite、 hibschite、α-C2SH and K2Ca(CO3)2; among which, K2Ca(CO3)2 was not the first production containing potassium, but K2Ca(CO3)2 was synthesized by the reaction among KOH、Ca(OH)2 and CO2. Whether the phase was synthesized was related to not only the reaction condition, but also their physicochemical properties; when the reaction condition was changed, the changes of different phases were different. The results of XRD and dissolution rate experiments explained the dissolution characteristic of every element of hydrothermal productions very well, and the relation between the dissolution rate of element and the phase of productions poured a good illumination on the production technology. The results of SEM and EDS showed that: hydrogarnet looked like spherical, and its surface was covered by some productions including K phase and Ca、Si phase; but the morphology of tobermorite was platy or lamellar or needlelike, and parts of Si in the structure of tobermorite were substituted by Al，and some K+ cations were inserted into the Ca interlayer of tobermorite at the same time. It was the first time that the interface between KAlSi3O8 and Ca(OH)2 was observed directly by SEM and EDS after the hydrothermal reaction, and the mechanism of hydrothermal reaction of KAlSi3O8 and Ca(OH)2 was further discussed. These results indicated that: the Ca-KAlSi3O8 intermediate compound was formed at first, and some K was released into the solution and KOH was produced at the same time; the C-S-H phase appeared before hydrogarnet, and then hydrogarnet was synthesized when the chemical reaction was carried on; if the reaction was carried on furthermore, α-C2SH、tobermorite and other C-S-H phases of different atom ratio appeared. The author found that the structure of KAlSi3O8 would be more drastically destroyed if there were some reactants, such as Ca(OH)2 which reacted with KAlSi3O8 and new phases were formed after the hydrothermal reaction between KAlSi3O8 and alkaline solution of equal ionic strength was finished. With the combination of calcination and hydrothermal reaction methods, the dissolution rate of products were greatly improved when the hydrothermal reaction was carried out after KAlSi3O8 and CaCO3 were calcined. Furthermore, the author has tentatively explored how to evaluate the effects of the differences of the activity of lime on the dissolution properties of hydrothermal products.

软土亚塑性模型与一维非线性固结理论研究

With the great development of Tianjing New Coastal District economy, people need more land to build and live. Land subsidence, which is caused by its special engineering geological conditions, has restricted the further development in the district. Soft soil consolidation is main factor of land subsidence ;thus , on the basis of consolidation theory, the paper make further study on soft soils one-dimension nonlinear consolidation which contains two parts:(1) the nonlinear consolidation of permeability coefficient and compressibility coefficient changing with time and depth, which means real one-dimension nonlinear consolidation;(2) the non-homogeneous consolidation of permeability coefficient and compressibility coefficient only changing with depth. Firstly, nonlinear characteristics of soft soils are elaborated. Hypoplastic theory is introduced to establish a modified soft soils nonlinear constitutive model; the nonlinear governing equation of compressibility coefficient is built, and the nonlinear characteristics of compressibility coefficient are analyzed. Secondly, Considering Load Fluctuation and soil thickness changing ,the consolidation characteristics of single layer is discussed in the paper; meanwhile, on the basis of the Davis and Raymond’s hypothesis and single layer nonlinear consolidation equation, the doubled-layer one-dimension nonlinear consolidation equation is also derived. The solution of the equation is obtained by analytical method, and the consolidation characteristics of doubled-layer soft soil nonlinear theory is also analyzed. Finally, based on assumption that permeability coefficient and compressibility coefficient is varying along depth, single layer soil one-dimension non-homogeneous consolidation differential equation is derived; and the approximate solution is obtained. Furthermore, the single layer non-homogeneous consolidation is extended to double layer non-homogeneous consolidation theory. By using parabolic differential scheme, the matrix equation is established; and the solution of the matrix equation is obtained by chase method. Consolidation characteristics of soil soft single (double) layer non-homogeneous consolidation theory and Terzaghi’s theory are also discussed.

水-岩土反应的力学与环境效应研究

The research on mechanical effects of water-rock and soil interaction on deformation and failure of rocks and soils involves three aspects of mechanics, physics and chemistry. It is the cross between geochemistry and rock mechanics and soil mechanics. To sum up, the mechanical effects of water-rock and soil interaction is related to many complex processes. Research in this respect has been being an important forward field and has broad prospects. In connection with the mechanism of the effects of the chemical action of water-rock on deformation and failure of rocks and soils, the research significance, the present state, the developments in this research domain are summarized. Author prospects the future of this research. The research of the subject should be possessed of important position in studying engineering geology and will lead directly to a new understand on geological hazard and control research. In order to investigation the macroscopic mechanics effects of chemical kinetics of water-rock interaction on the deformation and failure, calcic rock, red sandstone and grey granite reacting chemically with different aqueous solution at atmospheric temperature and atmospheric pressure are uniaxially compressed. The quantitative results concerning the changes of uniaxially compressive strength and elastic modulus under different conditions are obtained. It is found that the mechanical effects of water on rock is closely related to the chemical action of water-rock or the chemical damage in rock, and the intensity of chemical damage is direct ratio to the intensity of chemical action in water-rock system. It is also found that the hydrochemical action on rock is time-dependent through the test. The mechanism of permeation and hydrochemical action resulting in failure of loaded rock mass or propagation of fractures in rocks would be a key question in rock fracture mechanics. In this paper, the fracture mechanical effects of chemical action of water-rock and their time- and chemical environment-dependent behavior in grey granite, green granite, grey sandstone and red sandstone are analyzed by testing K_(IC) and COD of rock under different conditions. It is found that: ①the fracture mechanical effect of chemical action of water-rock is outstanding and time-dependent, and high differences exist in the influence of different aqueous solution, different rocks, different immersion ways and different velocity of cycle flow on the fracture mechanical effects in rock. ②the mechanical effects of water-rock interaction on propagation of fractures is consistent with the mechanical effects on the peak strength of rock. ③the intensity of the mechanical fracture effects increases as the intensity of chemical action of water-rock increases. ④iron and calcium ion bearing mineral or cement in rock are some key ion or chemical composition, and especially iron ion-bearing mineral resulting in chemical action of water-rock to be provided with both positive and negative mechanical effects on rock. Through the above two tests, we suggest that primary factors influencing chemical damage in rock consist of the chemical property of rock and aqueous solution, the structure or homogeneity of rocks, the flow velocity of aqueous solution passing through rock, and cause of formation or evolution of rock. The paper explores the mechanism on the mechanical effects of water-rock interaction on rock by using the theory of chemistry and rock fracture mechanics with chemical damage proposed by author, the modeling method and the energy point of view. In this paper, the concept of absorbed suction between soil grains caused by capillary response is given and expounded, and the relation and basic distinction among this absorbed suction, surface tension and capillary pressure of the soil are analyzed and established. The law of absorbed suction change and the primary factors affecting it are approached. We hold that the structure suction are changeable along with the change of the saturation state in unsaturated soils. In view of this, the concept of intrinsic structure suction and variable structure suction are given and expounded, and this paper points out: What we should study is variable structure suction when studying the effective stress. By IIIy κHH's theory of structure strength of soils, the computer method for variable structure suction is analyzed, the measure method for variable structure suction is discussed, and it reach the conclusions: ①Besides saturation state, variable structure suction is affected by grain composition and packing patter of grains. ②The internal relations are present between structure parameter N in computing structure suction and structure parameter D in computing absorbed suction. We think that some problems exit in available principle of effective stress and shear strength theory for unsaturated soil. Based on the variable structure suction and absorbed suction, the classification of saturation in soil and a principle of narrow sense effective stress are proposed for unsaturated soils. Based on generalized suction, the generalized effective stress formula and a principle of generalized effective stress are proposed for unsaturated soils. The experience parameter χ in Bishop's effective stress formula is defined, and the principal factors influencing effective stress or χ. The primary factor affecting the effective stress in unsaturated soils, and the principle classifying unsaturated soils and its mechanics methods analyzing unsaturated soils are discussed, and this paper points out: The theory on studying unsaturated soil mechanics should adopt the micromechanics method, then raise it to macromechanics and to applying. Researching the mechanical effects of chemical action of water-soil on soil is of great importance to geoenvironmental hazard control. The texture of soil and the fabric of soil mass are set forth. The tests on physical and mechanical property are performed to investigate the mechanism of the positive and negative mechanical effects of different chemical property of aqueous solution. The test results make clear that the plastic limit, liquid limit and plasticity index are changed, and there exists both positive and negative effects on specimens in this test. Based on analyzing the mechanism of the mechanical effects of water-soil interaction on soil, author thinks that hydrochemical actions being provided with mechanical effects on soil comprise three kinds of dissolution, sedimentation or crystallization. The significance of these tests lie in which it is recognized for us that we may improve, adjust and control the quality of soils, and may achieve the goal geological hazard control and prevention.The present and the significance of the research on environmental effects of water-rock and soil interaction. Various living example on geoenvironmental hazard in this field are enumerated. Following above thinking, we have approached such ideals that: ①changing the intensity and distribution of source and sink in groundwater flow system can be used to control the water-rock and soil interaction. ②the chemical action of water-rock and soil can be used to ameliorate the physical and mechanical property of rocks and soils. Lastly, the research thinking and the research methods on mechanical effects and environmental effects of water-rock and soil interaction are put forward and detailed.

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

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.

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

济阳坳陷及其南缘奥陶系碳酸盐岩潜山储集性及控制因素

The main research area of this thesis is Jiyang Depression in the Bohaiwan Basin and its southern margin. The object formation is Ordovician carbonate. The research is based on the outcrop observation and measurement of Ordovician carbonate and the drilling data of the oilfield. The internal reservoir characteristics of carbonate buried hill and its distribution were studied by comprehensive methods of sedimentology, reservoir geology and structural geology and technics of cathodoluminescence(CL)3electron microprobe,casting and C O isotope analysis etc. The influence depth of paleokarst facies formed during the Paleozoic is discriminated as 36-84m. The sollution porosity is well developed in paleokarst facies of Ordovician carbonate and is an important type of internal reservoir of buried hill. It may be infered that the fractures may be formed mainly during the Mesozoic and Cenozoic, they were not developed during the early Paleozoic when only micro-fractures might be created. The carbon and oxigen isotope analysis shows that the calcite cements in the fractures of Ordovician carbonate and secondary solution pores were related with meteoric water and three stages of fractures were divided. The reservoir space of Ordovician carbonate are mainly secondary porosity, cavern and fracture. The development of structural fracture was controlled by the lithology and tectonic background. More fractures exist in dolomite than that in limestone. There are also more fractures near the fault and the axis of fold. The development of porous reservoir is mainly controlled by the lithology and diagenesis, especially dolomitization and dissolution. It also results in the heterogeneity vertically. So the lithology is the basic factor for the forming of internal reservoir of buried hill and the tectogenesis and diagenesis are key factors to improve it. The porosity in carbonate might experienced solution-cementation-resolution or recementation. The porosity evolution history was a kind of historical dynamic equilibrium. The internal reservoir of Ordovician carbonate is the comprehensive result of constructive and/or destructive diagenesis. The worm's eye maps of the early Paleozoic and middle-upper Proterozoic were plotted. It was inferred that the paleostress field evoluted from NNW to NW during the Mesozoic and Cenozoic. Three types of buried hills can be divided: C-P/Pzi, Mz/ Pzi and E/ Pzi. The unconformity of the buried hill of E/ Pzi type, comparatively, was formed and reconstructed latestly, t he p orous r eservoir c ould b e w ell p reseved. T his c ondition w as v ery favorable t o t he migration and accumulation of oil and gas and could form upstanding association of source-reservoir-cap rocks. The buried hills of Mz/ Pzi and C-P/Pz] type were took second place.

A Modern Differential Geometric Approach to Shape from Shading

How the visual system extracts shape information from a single grey-level image can be approached by examining how the information about shape is contained in the image. This technical report considers the characteristic equations derived by Horn as a dynamical system. Certain image critical points generate dynamical system critical points. The stable and unstable manifolds of these critical points correspond to convex and concave solution surfaces, giving more general existence and uniqueness results. A new kind of highly parallel, robust shape from shading algorithm is suggested on neighborhoods of these critical points. The information at bounding contours in the image is also analyzed.

Interface of chip-based capillary electrophoresis-inductively coupled plasma-atomic emission spectrometry (CE-ICP-AES)

An interface of chip-based capillary electrophoresis (CE)-inductively coupled plasma-atomic emission spectrometry (ICP-AES) that is based on cross-flow nebulization has been developed. A polydimethylsiloxane (PDMS) CE-chip with conventional cross channel layout was used. A stainless steel tube was placed orthogonal to the exit of the CE separation channel for cross flow nebulization. A supplementary flow of buffer solution at the channel exit was used to improve nebulization efficiency. Two capillaries were inserted into the CE chip near the inlet of the separation channel for sample and buffer solution injection. Syringe pumps were used to manipulate the flow rate and flow direction of the sample, buffer, and supplementary buffer solution. Peak broadening due to the shape (bulb and tube-shaped) and size of the spray chambers was studied. The smaller tube-shaped spray chamber was used because of smaller peak broadening effect due to aerosol transport. The nebulization and transport efficiency of the CE-ICP interface was approximately 10%. Ba2+ and Mg2+ ions were eluted from the CE-chip within 30 s. Resolution of the Ba2+ and Mg2+ peaks was 0.7 using the chip-based CE-ICP-AES system.