基于Multi-Agent的生产过程建模与仿真技术研究

某些流程行业由于采用按配方进行分组加工的模式组织生产,在排产时存在多条路径调度优化的问题,应用一般的优化算法对于现场在线调度难以给出满意结果,而基于Agent的过程仿真在解决离散、非线性系统模拟方面有显著的优势,本文采用Agent的方法对生产过程建模,然后对方案组内的备选方案进行仿真,通过对比各方案的仿真结果找到最优的方案作为执行方案,为现场的优化排产提供决策支持。

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

陕北地区石家湾-李家岔探区三叠系延长组石油成藏条件分析

Shijiawan –Lijiacha area, lying on the northeastern part of the Shanbei Slope of Ordos Basin, was selected as studying area. The previous explorations proved that the 2nd segment and 6th segment of the Yanchang Formation are the most important oil-bearing formations. It is indicated that the sedimentary facies and reservoir characteristics restricted the hydrocarbon accumulation regularity. Therefore, with petrology methodologies, such as outcrop observation, core description, geophysical logging interpretation, thin section determination, scanning electron microscope, as well as rock property analysis, the reservoirs was were systematically studied and characterized. The sedimentary micro-facies, seals, reservoir-seal combines, migration pathways and entrapping modes were taken into account. The author tempted to establish a base for further studies on reservoirs and on petroleum geology, and to provide some reliably geological evidences for later prospect activities. It was found that the sediments in the 2nd and 3rd segments of the Yanchang Formation in Shijiawan –Lijiacha area were deposited in braided rivers, and most sandy-bodies were identified as channel sandbars. The 4+5th and 6th segments were principally deposited in deltaic-plain environment, consisting of corresponding sub-facies such as distributary channels, natural levee, crevasse-splay and marsh. The skeleton sandy-bodies were identified as sandy sediments of distributary channels. The sand grains in reservoir in studied area possess generally low mineralogical maturity and moderate structural maturity, and the form of pores may be classified into intergranular types and dissolved types. Most reservoirs of Yanchang Formation in Shijiawan –Lijiacha area belong to extreme low-porosity low-permeability ones (type III), and the 2nd sediments belongs to low permeability one (type II) and the 6th segment belong to super low-permeability one(type Ⅳ). The reservoirs in the 2nd segment behave more heterogeneous than those in the 6th segment. The statistic analysis results show that, for 6th and 4+5th segments, the high quality reservoir-seal combines may be found everywhere in the studied area except in the northwest and the southwest parts; and for 1st and 2nd segments, in the northeast, central and southwest parts Petroleum migration happened in the duration of the Early Cretaceous period in both lateral and vertical directions. The migration paths were mainly constructed by permeable sandy-bodies. The superimposed channel sandy-bodies consist of the principal part of the system of carriers. the vertical fractures, that may travel through the seals between reservoirs, offered the vertical paths for migrating oil. It may be synthesized that oil coming from south kitchens migrated first laterally in carriers in the 6th segment. When arrived at the studied area, oil will migration laterally or/and vertical within both the sandy-bodies and fractures, in a climbing-stair way. The results demonstrate that the oil was entrapped in traps structure-lithology and/or lithology traps. In some cases, the hydrodynamic force may help to trap oil. Accumulation of oil in the area was mainly controlled by sedimentary facies, seals, structure, and heterogeneity of reservoir in the 2nd, 4+5th and 6th segments. Especially, the oil distributions in both the 2nd and 6th segments were obviously influenced by seals in the 4+5th segment. The existence of seals in 1st segment seems important for accumulation in the 2nd segment.

井间地震速度层析成像方法与软件系统

Crosshole Seismic tomography has been broadly studied and applied in the fields of resource exploration and engineering exploration because of its special observing manner and better resolution than normal seismic exploration. This thesis will state the theory and method of Crosshole Seismic tomography. Basing on the previous studies，the thesis studied the initial velocity model，ray-tracing method, and developed the three-dimension tomography software. All the cells that a ray passes through are of the same velocities if the paths from transmitters to receivers are straight. The cells that the each ray passes through are recorded, and rays that pass through each cell are calculated. The ray average velocity which passes through a cell is set as the cell velocity. Analogously we can make a initial node velocity model because the velocity sum is calculated on the all cells which own to a certain node, and the cell number is summed about each nodes，the ratio of the velocity sum to the all cells number is set as the node velocity. The inversion result from the initial node velocity model is better than that of the average velocity model. Ray-bending and Shortest Path for Rays (SPR) have shortcomings and limitations respectively. Using crooked rays obtained from SPR rather than straight lines as the starting point can not only avoid ray bending converging to the local minimum travel time path, but also settle the no smooth ray problem obtained by SPR. The hybrid method costs much computation time, which is roughly equal to the time that SPR expends. The Delphi development tool based on the Object Pascal language standard has an advantage of object-oriented. TDTOM (Three Dimensions Tomography) was developed by using Delphi from the DOS version. Improvement on the part of inversion was made, which bring faster convergence velocity. TDTOM can be used to do velocity tomography from the first arrival travel time of the seismic wave, and it has the good qualities of friendly user interface and convenient operation. TDTOM is used to reconstruct the velocity image for a set of crosshole data from Karamay Oil Field. The geological explanation is then given by comparing the inversion effects of different ray-tracing methods. High velocity zones mean the cover of oil reservoir, and low velocity zones correspond to the reservoir or the steam flooding layer.

旋转轴对称弱各向异性介质模型的射线追踪方法

Ray tracing is a rapid and effective method for wave field calculation. Not only in the field of seismic-wave theory, but also in the field of seismic inversion and migration imaging,the seismic ray tracing method has become one of the most important methods. In anisotropic media, group velocity and phase velocity have different propagation directions. The seismic wave propagates along the direction of group velocity , it does not depend on the direction of phase velocity. Ray angle is a complex function with respect to phase angle, it is difficult to measure and calculate. But most rocks are weak anisotropic, so the expression of phase velocity can be simplified greatly. Based on the approximate expression of phase velocity this thesis for rotating axisymmetric weak anisotropic media deduces an expression of the partial derivative of phase velocity and an expression of group velocity with the method of linear approximation. This paper uses the fourth order Runge-Kutta method together with the two-dimensional interpolation and linear interpolation to obtain the parameters of the physical locations. At last the paths of seismic wave in rotating axisymmetric weak anisotropic media are computed. According to the analysis of the computational results, it indicates that the method developed in this paper has strong adaptability, high computational efficiency and high accuracy for rotating axisymmetric weak anisotropic media.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

东喜马拉雅构造结新生代碰撞变形与岩石圈演化

The Eastern Himalayan Syntaxis (EHS) is one of the strongest deformation area along the Himalayan belt resulted from the collision between Indian plate and the Eurasian Plate since the 50～60Ma, and has sensitivity tracked and preserved the whole collisional processes. It should depend on the detail geological investigations to establish the deformational accommodate mode, and the uplift history, to elucidate the deep structure and the crust-mantle interaction of the Tibet Plateau of the EHS. The deep-seated (Main Mantle Thrusts) structures were exhumed in the EHS. The MMT juxtapose the Gangdese metamorphic basement and some relic of Gangdese mantle on the high Himalayan crystalline series. The Namjagbawa group which is 1200～1500Ma dated by U/Pb age of zircon and the Namla group which is 550Ma dated by U/Pb age of zircon is belong to High Himalayan crystalline series and Gangdese basement respectively. There is some ophiolitic relic along the MMT, such as metamorphic ocean mantle peridotite and metamorphic tholeiite of the upper part of ocean-crust. The metamorphic ocean mantle peridotites (spinel-orthopyroxene peridotite) show U type REE patterns. The ~(87)Sr/~(86)Sr ratios were, 0.709314～0.720788, and the ~(143)Nd/~(144)Nd ratios were 0.512073～0.512395, plotting in the forth quadrant on the ~(87)Sr/~(86)Sr-~(143)Nd/~(144)Nd isotope diagram. Some metamorphic basalt (garnet amphibolite) enclosures have been found in the HP garnet-kynite granulite. The garnet amphibolites can be divided two groups, the first group is deplete of LREE, and the second group is flat or rich LREE, and their ~(87)Sr/~(86)Sr, ~(143)Nd/~(144)Nd ratios were 0.70563～0.705381 and 0.512468～0.51263 respectively. Trace element and isotopic characteristics of the garnet amphibolites display that they formed in the E-MORB environment. Some phlogolite amphibole harzburgites, which exhibit extensive replacement by Phl, Amp, Tc and Dol etc, were exhumed along the MMT. The Phl-Amp harzburgites are rich in LREE and LILE, such as Rb, K etc, and depletes Eu (Eu~* = 0.36 ～ 0.68) and HFSE, such as Nb, Ta, Zr, Hf, P, Ti etc. The trace element indicate that the Phl-Amp harzburgites have island arc signature. Their ~(87)Sr/~(86)Sr are varied from 0.708912 to 0.879839, ~(143)Nd/~(144)Nd from 0.511993 to 0.512164, ε Nd from- 9.2 to - 12.6. Rb/Sr isochrone age of the phlogolite amphibole harzburgite shows the metasomatism took place at 41Ma, and the Amp ~(40)Ar/~(39)Ar cooling age indcate the Phl-Amp harzburgite raising at 16Ma. There is an intense crust shortening resulted from the thrust faults and folds in the Cayu block which is shortened more 120km than that of the Lasha block in 35～90Ma. With the NE corner of the India plate squash into the Gangdese arc, the sinistral Pai shear fault and the dextral Aniqiao shear fault on the both sides of the Great bent of Yalun Zangbu river come into active in 21～26Ma. On the other hand, the right-lateral Gongrigabu strike-slip faults come into activity at the same period, a lower age bound for the Gongrigabu strike-slip fault is estimated to be 23～24Ma from zircon of ion-probe U/Pb thermochronology. The Gongrigabu strike-slip faults connect with the Lhari strike-slip fault in the northwestern direction and with the Saganing strike-slip at the southeastern direction. Another important structure in the EHS is the Gangdese detachment fault system (GDS) which occurs between the sedimental cover and the metamorphic basement. The lower age of the GDS is to be 16Ma from the preliminary 40Ar/39Ar thermochronology of white mica. The GDS is thought to be related to the reverse of the subducted Indian crust and the fast uplift of the EHS. Structural and thermochronology investigation of the EHS suggest that the eastern Tibet and the western Yunnan rotated clockwise around the EHS in the period of 35～60Ma. Later, the large-scale strike-slip faults (RRD, Gaoligong and Saganing fault) prolongate into the EHS, and connect with the Guyu fault and Gongrigabu fault, which suggest that the Indianchia block escape along these faults. Two kind of magmatic rocks in the EHS have been investigated, one is the mantle-derived amphibole gabbro, dioposide diorite and amphibole diorite, another is crust origin biotit-garnet adamellite, biotit-garnet granodiorite and garnet-amphibole-biotite granite. The amphibole gabbro dioposite diorite and amphibole diorite are rich in LREE, and LILE, such as Ba, Rb, Th, K, Sr etc, depleted in HFSE, such as Nb, Ta, Zr, Hf, Ti etc. The ratio of ~(87)Sr/~(86)Sr are from 0.7044 to 0.7048, ~(143)Nd/~(144)Nd are from 0.5126 to 0.5127. The age of the mantle origin magamatic rocks, which result from the partial melt of the raising and decompression anthenosphere, is 8Ma by ~(40)Ar/~(39)Ar dating of amphibole from the diorite. The later crust origin biotite-garnet adamellite, biotite-garnet granodiorite and garnet-amphibole-biotite granite are characterized by aboudance in LREE, and strong depletion of Eu. The ratios of ~(87)Sr-~(86)Sr are from 0.795035 to 0.812028, ~(143)Nd/~(144)Nd from 0.51187 to 0.511901. The ~(40)Ar/~(39)Ar plateau age of the amphibole from the garnet-amphibole-biotite granite is 17.5±0.3Ma, and the isochrone age is 16.8±0.6Ma. Their geochemical characteristics show that the crust-derived magmatic rocks formed from partial melting of the lower curst in the post-collisional environment. A group of high-pressure kaynite-garnet granulites and enclave of high-pressure garnet-clinopyroxene grnulites and calc-silicate grnulites are outcroped along the MMT. The peak metamorphic condition of the high-pressure granulites yields T=800～960 ℃, P=1.4～1.8Gpa, corresponding the condition of 60km depth. The retrograde assemblages of the high-pressure grnulites occur at the condition of T=772.3～803.3 ℃, P=0.63～0.64Gpa. The age of the peak metamorphic assemblages are 45 ～ 69Ma indicated by the zircon U/Pb ion-plobe thermochronology, and the retrograde assemblage ages are 13～26Ma by U/Pb, ~(40)Ar/~(39)Ar thermochronology. The ITD paths of the high-pressure granulites show that they were generated during the tectonic thickening and more rapid tectonic exhumation caused by the subducting of the Indian plate and subsequent break-off of the subducted slab. A great deal of apatite, zircon and sphene fission-track ages, isotopic thermochronology of the rocks in the EHS show that its rapid raising processes of the EHS can be divided into three main periods. There are 35～60Ma, 13～25Ma, 0～3Ma. 3Ma is a turn in the course of raising in the EHS which is characterized by abruptly acceleration of uplifting. The uplift ratios are lower than 1mm .a~(-1) before 3Ma, and higher than 1mm .a~(-1) with a maximum ratio of 30mm .a~(-1) since 3Ma. The bottom (knick point) of the partial anneal belt is 3.8km above sea level in the EHS, and correspond to age of 3Ma determined by fission-track age of apatite. The average uplift ratio is about 1.4 mm .a~(-1) below the knick point. The EHS has raised 4.3km from the surface of 2.36km above sea level since 3Ma estimated by the fossil partial anneal belt of the EHS. We propose a two-stage subduction model (B+A model) basing on Structural, thermochronological, magmatical, metamorphic and geophysical investigations of the EHS. The first stage is the subduction of the Indian continental margin following after the subduction of the Tethys Ocean crust and subsequent collision with the Gangdese arc, and the second stage is the Indian crust injecting into the lower crust and upper mantle of the Tibet plateau. Slab break-off seems to be occurred between these two stages.

油气运聚动力学的模拟分析——以准噶尔盆地为例

The petroleum migration, happening in the geologic past, is the very important and complex dynamic processes in the petroleum systems. It plays a linking role among all static factors in a system. The accumulation is in fact the result of the petroleum migration. For the petroleum geology, the dynamics research of the petroleum migration refers to the mechanism and process research, as well as the use of the quantitative methods. In this thesis, combining with the qualitative analysis and quantitative modeling, the author manages to discuss theoretically some key problems dealing with migration processes, which have not been solved yet, and to apply the studied results in petroleum system analysis in actual basins. The basin analysis offers the base of the numerical modeling for geological phenomena occurring in sedimentary basins, that consists of the sedimentary facies analysis, the section reconstructing technique, eroded thickness estimating, etc. The methods to construct the geologic model, which is needed in the research of oil and gas migration and accumulation, are discussed. The basin analysis offers also the possibility for the latter modeling works to get and select the parameters, such as stratum's thickness, age, stratigraphy etc. Modeling works were done by using two basin modeling softwares: Basin_Mod and TPC_Mod. The role of compaction during the secondary migration and the heterogeneity of migrating paths within the clastic carrier are modeled. And the conclusions were applied in the migration studies in the Jungaer Basin, lying on the Northwest part of the China. To construct a reliable migration model, the author studied the characteristics of the sedimentation, the pore fluid pressure evolution, as well as the distribution and the evolution of fluid potential, following the tectonic evolution of the Jungaer Basin. The geochemical prospecting results were used to evidence and to calibrate the migration processes: the oil-source correlation, the distribution of the properties of oil, gas and water. Finally, two important petroleum systems, Permian one and Jurassic one were studied and identified, according, principally, to the studies on the petroleum migration within the Jungaer Basin. Since the oil, as well as the gas, moves mainly in separate phase during the secondary migration, their migrating behaviors would be determined by the dynamics conditions of migration, including the driving forces and pathways. Based on such a consideration, the further understandings may be acquired: the roles played by permeable carriers and low-permeable source rock would be very different in compaction, overpressure generation, petroleum migration, and so on. With the numerical method, the effect of the compaction on the secondary migration was analyzed and the results show that the pressure gradient and the flux resulted from compaction are so small that could be neglected by comparing to the buoyancy of oil. The main secondary migration driving forces are therefore buoyancy and capillary within a hydrostatic system. Modeling with the commercial software-Basin_Mod, the migration pathways of petroleum in clastic carriers seem to be inhomogeneous, controlled by heterogeneity of the driving force, which in turn resulted from the topography of seals, the fabrics and the capillary pressure of the clastic carriers. Furthermore, the direct and indirect methods to study fault-sealing properties in the course of migration were systemically summarized. They may be characterized directly by lithological juxtaposition, clay smear and diagenesis, and indirectly the comparing the pressures and fluid properties in the walls at two apartments of a fault. In Jungaer Basin, the abnormal pressures are found in the formations beneath Badaowan or Baijantan Formation. The occurrence of the overpressure seems controlled by the stratigraphy. The rapid sedimentation, tectonic pressuring, clay sealing, chemical diagensis were considered as the principal pressuring mechanisms. The evolution of fluid pressure is influenced differently at different parts of the basin by the tectonic stresses. So the basin appears different pressure evolution cycles from each part to another during the geological history. By coupling the results of thermal evolution, pressure evolution and organic matter maturation, the area and the period of primary migration were acquired and used to determine the secondary migration time and range. The primary migration in Fengcheng Formation happened from latter Triassic to early Jurassic in the main depressions. The main period of lower-Wuerhe Formation was at latter Jurassic in Changji, Shawan and Pen-1-jing-xi Depression, and at the end of early Cretaceous in Mahu Depression. The primary migration in Badaowan and Sangonghe Formation is at the end of early-Cretaceous in Changji Depression. After then, the fluid potential of oil is calculated at the key time determined from area and time of the primary migration. Generally, fluid potential of oil is high in the depressions and low at the uplifts. Synthetically, it is recognized that the petroleum migration in the Jungaer Basin is very complex, that leads us to classify the evolution of petroleum systems in Northwestern China as a primary stage and a reformed one. The remigration of accumulated petroleum, caused by the reformation of the basin, results in the generation of multiple petroleum systems. The faults and unconformities are usually the linkers among the original petroleum systems. The Permian petroleum system in Jungaer Basin is such a multiple petroleum system. However, the Jurassic petroleum system stays still in its primary stage, thought the strong influences of the new tectonic activities.

吐合盆地鲁克沁构造带断块油藏形成机制及油藏预测

Lukeqin arc belt is a compound structure generated by multi-movements and composed of 6 sub-structural zones, which are connected by Huoyanshan Mountain. General characteristics of the arc belt are multi-patterns of structure, multi-phases for petroleum, multi-types of trap and multi-layers for reservoirs. As a part of the eastern Lukeqin arc belt located on the south of Taibei depression, Lukeqin structural zone behaves as a complex faulted-fold zone, in which the formation and distribution of hydrocarbons are controlled by structures. As the dominant source of dynamics for the second migration of hydrocarbon, structure stress field is closely related with the potentials of hydrodynamics. Results derived from the simulations of stress field by finite element method indicate that the northwest tending faults prefer seal to the northeast tending ones. The reason is that the northwest tending faults were squeezed more strongly than the northeast tending ones. Therefor, the northeast tending faults become always the paths for oil to migrate southeastward. Lukeqin structural zone is the main site for oil to concentration because it is surrounded by high stress. Situated on the front of the foreland basin of Turpan-Hami, Lukeqing arc belt is a dam to hold back the southward migrating oil from Shengbei depression. The axis line of Shenquan-Shengnan-Yanmuxi, Lukeqin and Yubei controls the migrating paths and concentrating process of oil and gas. Results derived from stress simulation and structure analyses indicate consistently that both Yubei and Lukeqin structural zones are the favorite areas for oil to migrate. The generally southward paths for oil to migrate out of Taibei depression can be two ways. One of them is from Taibei depression to Yubei structural zone and the other is from Taibei depression to Lukeqin structural zone. By the both ways, oil migrated upward along the faults and southeastward along the structural axis to concentrate in either Permian or Triassic system. The newly ascertained path for oil migration, which is accurately southeastward instead of coarsely southward, indicates the directions for further explorations on the compound Lukeqin block zone. Five kinds of seal models of fault are all found in Lukeqin block zone by studying the seal features of faults occurred in the zone. Having studied the fault seal and their controlling factors by fuzzy set method, the paper deems that the northwest tended faults are better than the northeast tended ones for oil to concentrate. The most important factors to decide the seal extent of faults in this zone are the characteristics of main stress and fluids instead of capillary pressure differences between the two sides of fault and smear mud factors. There exist seal differences not only between the faults of different time but also between the sections within a fault due to the variation of depths, strata and positions. The general distribution rules of reservoirs were dominated by the seal characteristics of a fault during the time reservoirs formed. While the current features of fault seal decide the conservation of reservoirs and heights of oil accumulations. Seal or not of a fault is not absolute because the essential for fault to seal is the distribution of permeability of fault zone. Therefor, the multi cyclical activities of faults create the space-time variation of seal features of the fault. Totally, the seal extent of the faults within the area is not as perfect as to accumulate ordinary crude. Crude oil can only be sealed when it becomes viscous. Process for crude oil to become viscous and viscous happened strongly because of the fault-fold movements. Shallowly burying and even revealing of the objective layers of the reservoirs made the crude oil to be thickened by water washing biologically degradation and oxidation degradation. The northwestward deepening during or after the reservoir formation of the structural zone provided the power for oil to migrate one or more times. The main reason for oil accumulation is the formation of Lukeqin block zone during Xishanyao stage, middle Jurassic Period, Early Yanshanian Movement. While the main reason for reservoir conservation is the placidity of Triassic blocks after the formation of reservoirs. Contrasting to former opinions, it is concluded that the reservoirs in Lukeqin zone, including viscous reservoirs, were formed by one time but not more times. So the author proposes the opinion that the reservoirs of viscous oil were formed by viscous oil migration under the conditions of aptitude sets of fault seals controlled by fluid and other factors. To grope the distribution rules outside Taibei depression and discuss the formation mechanism of Anjurassic reservoirs, it is necessary to study the dominate factors for the formation of reservoirs in Lukeqin structural zone such as structural stress, fault seals and thickening mechanism of crude oil. Also, the necessary studies are the key to break through the Taibei depression and Anjurassic systems. Therefor, they are significant for the future exploration and reserve increasing of hydrocarbon within the Turpan-Hami basin. The paper studied the distribution rules of block reservoirs and forecasted the favorable zones for further exploration in Turpan-Hami basin. Conclusions can be useful for not only the exploration in the area but also the theory consult in the adjacent areas.

成年女儿照料老年父母的心理体验及关系质量作用的研究

Previous researches about family caregiving revealed that caregiving has both negative and positive effects on caregivers’ well-being. Based on Lawton’s two-factor model, this study aims at examining how caring for old parents would affect adult daughters’ psychological well-being. According to Lawton, objective stressors as caregiving would arouse two different kinds of caregivers’ subjective appraisal, i.e., negative appraisal and positive appraisal, which in turn correlate with the negative and positive dimensions of caregivers’ psychological well-being, respectively. There were two main purposes of this study: a) to verify both the negative and positive paths in the two-factor model and their relatively independence; and b) to examine the effects of relationship quality between caregiver and care-recipient on those paths. The results are as follows: 1) Caregiving stressors have significant positive predictive effect on caregivers’ negative appraisal, but have no direct effect on caregivers’ positive appraisal. 2) Caregivers’ negative appraisal has significant positive predictive effect on their negative emotional experience, while caregivers’ positive appraisal has significant positive predictive effect on their positive emotional experience. 3) Certain dimensions of relationship quality, including the Appreciation and General Appraisal, have significant negative predictive effect on caregivers’ negative appraisal, and have significant positive predictive effect on caregivers’ positive appraisal. 4) The Appreciation dimension of relationship quality moderates the path from caregiving demands to caregivers’ burden; and the General Appraisal of relationship quality moderates the path from caregivers’ positive appraisal to life satisfaction. Based on the above results, the researcher concluded that a) both the negative path and positive path exist in caregiving process, and they are relatively independent from each other; and b) relationship quality does moderate certain paths in the model. Meanwhile, the main effect of relationship quality on caregivers’ experience is also significant and more remarkable. This study attempts to explain these results in terms of coping resources. Both relationship quality and many other factors might be explained as resources that caregivers utilize to cope with stress of caregiving. With more resources, caregivers tend to appraise more positively, and less negatively, and vice versa. However, the resources which might affect caregivers’ positive appraisal, as well as the ways they work, may be different from that affect caregivers’ negative appraisal.

老化态度对老年人记忆作业影响的研究

It was a study of the effect of non-cognitive factors on cognitive aging. The purpose of the study was to research the effect of attitude towards aging on memory performance of the elderly, and trying to explore the mechanism of the attitude towards aging on memory performance. The entire study was divided into three phases and included six experiments. The questionnaire of attitude towards aging was compiled in phase 1. Then the validity and reliability of it was tested in the elderly. The self-esteem scale was also tested as above. These prepare the instruments for the study. It was a static study in phase 2. Through comparing the old groups in different environment, analyzing relationship between different non-cognitive factors and memory performance in the old subjects in identical environment, the paths of the effect of the attitude towards aging on memory performance was explored. It was a dynamic study in phase 3. By changing attitude towards aging of the elderly with direct thinking (positive, negative), the effect on memory performance was observed and the mechanism was explored. Main results: 1. The validity and reliability of the questionnaire towards aging were good. Its construct had two factors (positive, negative). Both reliability of retest and split test were over 0.90. 2. Direct thinking task was able to change the attitude towards aging effectively in old or young adults. The scores of the attitudes towards aging between positive and negative thinking subjects were different significantly. 3. The effect of attitude towards aging on memory performance was significant. The results in static or dynamic condition showed that the elderly with higher score of the questionnaire of attitude towards aging (more positive), their memory self-efficacy and memory performance were higher. In the opposite (negative), they were lower. 4. Mechanism of the effect of the attitude towards aging on memory performance: in static condition there were two paths i.e. direct path and through memory self-efficacy; in dynamic condition, besides the two paths, emotion change constituted another path. 5. Different non-cognitive factors have different effect on psychological function. The results of this study found that self-esteem has no significant correlation with memory performance, but correlates with depression; The factor of education has effect on many non-cognitive factors and memory performance.

Battling Reality

In the four years that the MIT Mobile Robot Project has benn in existence, we have built ten robots that focus research in various areas concerned with building intelligent systems. Towards this end, we have embarked on trying to build useful autonomous creatures that live and work in the real world. Many of the preconceived notions entertained before we started building our robots turned out to be misguided. Some issues we thought would be hard have worked successfully from day one and subsystems we imagined to be trivial have become tremendous time sinks. Oddly enough, one of our biggest failures has led to some of our favorite successes. This paper describes the changing paths our research has taken due to the lessons learned from the practical realities of building robots.

Transition Space

Informal causal descriptions of physical systems abound in sources such as encyclopedias, reports and user's manuals. Yet these descriptions remain largely opaque to computer processing. This paper proposes a representational framework in which such descriptions are viewed as providing partial specifications of paths in a space of possible transitions, or transition space. In this framework, the task of comprehending informal causal descriptions emerges as one of completing the specifications of paths in transition space---filling causal gaps and relating accounts of activity varied by analogy and abstraction. The use of the representation and its operations is illustrated in the context of a simple description concerning rocket propulsion.

Trajectory Mapping ("TM''): A New Non-Metric Scaling Technique

Trajectory Mapping "TM'' is a new scaling technique designed to recover the parameterizations, axes, and paths used to traverse a feature space. Unlike Multidimensional Scaling (MDS), there is no assumption that the space is homogenous or metric. Although some metric ordering information is obtained with TM, the main output is the feature parameterizations that partition the given domain of object samples into different categories. Following an introductory example, the technique is further illustrated using first a set of colors and then a collection of textures taken from Brodatz (1966).