空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

锦屏电站地下厂房卸荷岩体参数研究

It is a basic work to ascertain the parameters of rock mass for evaluation about stability of the engineering. Anisotropism、inhomogeneity and discontinuity characters of the rock mass arise from the existing of the structural plane. Subjected to water、weathering effect、off-loading, mechanical characters of the rock mass are greatly different from rock itself, Determining mechanical parameters of the rock mass becomes so difficult because of structure effect、dimension effect、rheological character, ‘Can’t give a proper parameter’ becomes one of big problems for theoretic analysis and numerical simulation. With the increment of project scale, appraising the project rock mass and ascertaining the parameters of rock mass becomes more and more important and strict. Consequently, researching the parameters of rock mass has important theoretical significance and actual meaning. The Jin-ping hydroelectric station is the first highest hyperbolic arch dam in the world under construction, the height of the dam is about 305m, it is the biggest hydroelectric station at lower reaches of Yalong river. The length of underground factory building is 204.52m, the total height of it is 68.83m, the maximum of span clearance is 28.90m. Large-scale excavation in the underground factory of Jin-ping hydroelectric station has brought many kinds of destructive phenomenon, such as relaxation、spilling, providing a precious chance for study of unloading parameter about rock mass. As we all know, Southwest is the most important hydroelectric power base in China, the construction of the hydroelectric station mostly concentrate at high mountain and gorge area, basically and importantly, we must be familiar with the physical and mechanical character of the rock mass to guarantee to exploit safely、efficiently、quickly, in other words, we must understand the strength and deformation character of the rock mass. Based on enough fieldwork of geological investigation, we study the parameter of unloading rock mass on condition that we obtain abundant information, which is not only important for the construction of Jin-ping hydroelectric station, but also for the construction of other big hydroelectric station similar with Jin-ping. This paper adopt geological analysis、test data analysis、experience analysis、theory research and Artificial Neural Networks （ANN） brainpower analysis to evaluate the mechanical parameter, the major production is as follows: (1)Through the excavation of upper 5-layer of the underground powerhouse and the statistical classification of the main joints fractures exposed, We believe that there are three sets of joints, the first group is lay fracture, the second group and the fourth group are steep fracture. These provide a strong foundation for the following calculation of and analysis; (2)According to the in-situ measurement about sound wave velocity、displacement and anchor stress, we analyses the effects of rock unloading effect，the results show a obvious time-related character and localization features of rock deformation. We determine the depth of excavation unloading of underground factory wall based on this. Determining the rock mass parameters according to the measurement about sound wave velocity with characters of low- disturbing、dynamic on the spot, the result can really reflect the original state, this chapter approximately the mechanical parameters about rock mass at each unloading area; (3)Based on Hoek-Brown experienced formula with geological strength index GSI and RMR method to evaluate the mechanical parameters of different degree weathering and unloading rock mass about underground factory, Both of evaluation result are more satisfied; (4)From the perspective of far-field stress, based on the stress field distribution ideas of two-crack at any load conditions proposed by Fazil Erdogan (1962)，using the strain energy density factor criterion (S criterion) proposed by Xue changming(1972)，we establish the corresponding relationship between far-field stress and crack tip stress field, derive the integrated intensity criterion formula under the conditions of pure tensile stress among two line coplanar intermittent jointed rock，and establish the corresponding intensity criterion for the exploratory attempt; (5)With artificial neural network, the paper focuses on the mechanical parameters of rock mass that we concerned about and the whole process of prediction of deformation parameters, discusses the prospect of applying in assessment about the parameters of rock mass，and rely on the catalog information of underground powerhouse of Jinping I Hydropower Station, identifying the rock mechanics parameters intellectually，discusses the sample selection, network design, values of basic parameters and error analysis comprehensively. There is a certain significance for us to set up a set of parameters evaluation system，which is in construction of large-scale hydropower among a group of marble mass.

地震P-S波时差耦合作用的斜坡崩滑效应研究

China locates between the circum-Pacific and the Mediterranean-Himalayan seismic belt. The seismic activities in our country are very frequent and so are the collapses and slides of slope triggered by earthquakes. Many collapses and slides of slope take place mainly in the west of China with many earthquakes and mountains, especially in Sichuan and Yunnan Provinces. When a strong earthquake happening, the damage especially in mountains area caused by geological hazards it triggered such as rock collapses, landslides and debris flows is heavier than that it caused directly. A conclusion which the number of lives lost caused by geological hazards triggered by a strong earthquake in mountains area often accounts for a half even more of the total one induced by the strong earthquake can be made by consulting the statistical loss of several representative earthquakes. As a result, geological hazards such as collapses and slides of slope triggered by strong earthquakes attract wide attention for their great costs. Based on field geological investigation, engineering geological exploration and material data analysis, chief conclusions have been drawn after systematic research on formation mechanism, key inducing factors, dynamic characteristics of geological hazards such as collapses and slides of slope triggered by strong earthquakes by means of engineering geomechanics comprehensive analysis, finite difference numerical simulation test, in-lab dynamic triaxial shear test of rock, discrete element numerical simulation. Based on research on a great number of collapses and landslides triggered by Wenchuan and Xiaonanhai Earthquake, two-set methods, i.e. the method for original topography recovering based on factors such as lithology and elevation comparing and the method for reconstructing collapsing and sliding process of slope based on characteristics of seism tectonic zone, structural fissure, diameter spatial distribution of slope debris mass, propagation direction and mechanical property of seismic wave, have been gotten. What is more, types, formation mechanism and dynamic characteristics of collapses and slides of slope induced by strong earthquakes are discussed comprehensively. Firstly, collapsed and slided accumulative mass is in a state of heavily even more broken. Secondly, dynamic process of slope collapsing and sliding consists of almost four stages, i.e. broken, thrown, crushed and river blocked. Thirdly, classified according to failure forms, there are usually four types which are made up of collapsing, land sliding, land sliding-debris flowing and vibrating liquefaction. Finally, as for key inducing factors in slope collapsing and sliding, they often include characteristics of seism tectonic belts, structure and construction of rock mass, terrain and physiognomy, weathering degree of rock mass and mechanical functions of seismic waves. Based on microscopic study on initial fracturing of slope caused by seismic effect, combined with two change trends which include ratio of vertical vs. horizontal peak ground acceleration corresponding to epicentral distance and enlarging effect of peak ground acceleration along slope, key inducing factor of initial slope fracturing in various area with different epicentral distance is obtained. In near-field area, i.e. epicentral distance being less than 30 km, tensile strength of rock mass is a key intrinsic factor inducing initial fracturing of slope undergoing seismic effect whereas shear strength of rock mass is the one when epicentral distance is more than 30 km. In the latter circumstance, research by means of finite difference numerical simulation test and in-lab dynamic triaxial shear test of rock shows that initial fracture begins always in the place of slope shoulder. The fact that fracture strain and shear strength which are proportional to buried depth of rock mass in the place of slope shoulder are less than other place and peak ground acceleration is enlarged in the place causes prior failure at slope shoulder. Key extrinsic factors inducing dynamic fracture of slope at different distances to epicenter have been obtained through discrete element numerical simulation on the total process of collapsing and sliding of slope triggered by Wenchuan Earthquake. Research shows that combined action of P and S seismic waves is the key factor inducing collapsing and sliding of slope at a distance less than 64 km to initial epicenter along earthquake-triggering structure. What is more, vertical tensile action of P seismic wave plays a leading role near epicenter, whereas vertical shear action of S seismic wave plays a leading role gradually with epicentral distance increasing in this range. On the other hand, single action of P seismic wave becomes the key factor inducing collapsing and sliding of slope at a distance between 64 km and 216 km to initial epicenter. Horizontal tensile action of P seismic wave becomes the key factor gradually from combined action between vertical and horizontal tensile action of P seismic wave with epicentral distance increasing in this distance range. In addition, initial failure triggered by strong earthquakes begins almost in the place of slope shoulder. However, initial failure beginning from toe of slope relates probably with gradient and rock occurrence. Finally, starting time of initial failure in slope increases usually with epicentral distance. It is perhaps that the starting time increasing is a result of attenuating of seismic wave from epicenter along earthquake-triggering structure. It is of great theoretical and practical significance for us to construct towns and infrastructure in fragile geological environment along seism tectonic belts and conduct risk management on earthquake-triggered geological hazards by referring to above conclusions.

开发早期油藏动态地质模型与油藏工程优化研究——以埕岛油田馆上段油藏为例

In order to developing reservoir of Upper of Ng at high-speed and high-efficient in Chengdao oilfield which is located in the bally shallow sea, the paper builds up a series of theory and means predicting and descripting reservoir in earlier period of oilfield development. There are some conclusions as follows. 1. It is the first time to form a series of technique of fine geological modeling of the channel-sandy reservoir by means of mainly seismic methods. These technique include the logging restriction seismic inversion, the whole three dimension seismic interpretation, seismic properties analysis and so on which are used to the 3-dimension distributing prediction of sandy body, structure and properties of the channel reservoir by a lot of the seismic information and a small quantity of the drilling and the logging information in the earlier stage of the oil-field development. It is the first time that these methods applied to production and the high-speed development of the shallow sea oilfield. The prediction sandy body was modified by the data of new drilling, the new reservoir prediction thinking of traced inversion is built. The applied effect of the technique was very well, according to approximately 200 wells belonging to 30 well groups in Chengdao oilfield, the drilling succeeded rate of the predicting sandy body reached 100%, the error total thickness only was 8%. 2. The author advanced the thinking and methods of the forecasting residual-oil prediction at the earlier stage of production. Based on well data and seismic data, correlation of sediment units was correlated by cycle-correlation and classification control methods, and the normalization and finely interpretation of the well logging and sedimentation micro-facies were acquired. On the region of poor well, using the logging restriction inversion technique and regarding finished drilling production well as the new restriction condition, the sand body distributing and its property were predicted again and derived 3-dimension pool geologic model including structure, reservoir, fluid, reservoir engineering parameter and producing dynamic etc. According to the reservoir geologic model, the reservoir engineering design was optimized, the tracking simulation of the reservoir numerical simulation was done by means of the dynamic data (pressure, yield and water content) of development well, the production rule and oil-water distributing rule was traced, the distributing of the remaining oil was predicted and controlled. The dynamic reservoir modeling method in metaphase of development was taken out. Based on the new drilling data, the static reservoir geologic model was momentarily modified, the research of the flow units was brought up including identifying flow units, evaluating flow units capability and establishing the fine flow units model; according to the dynamic data of production and well testing data, the dynamic tracing reservoir description was realized through the constant modification of the reservoir geologic model restricted these dynamic data by the theory of well testing and the reservoir numerical simulation. It was built the dynamic tracing reservoir model, which was used to track survey of the remaining oil on earlier period. The reservoir engineering tracking analysis technique on shallow sea oilfield was founded. After renewing the structure history since tertiary in Chengdao area by the balance section technique and estimating the activity character of the Chengbei fault by the sealing fault analysis technique, the meandering stream sediment pattern of the Upper of Ng was founded in which the meandering border was the uppermost reservoir unit. Based on the specialty of the lower rock component maturity and the structure maturity, the author founded 3 kinds of pore structure pattern in the Guanshang member of Chengdao oil-field in which the storing space mainly was primary (genetic) inter-granular pore, little was secondary solution pore and the inter-crystal pore tiny pore, and the type of throat mainly distributed as the slice shape and the contract neck shape. The positive rhythmic was briefly type included the simple positive rhythm, the complex positive rhythm and the compound rhythm. Interbed mainly is mudstone widely, the physical properties and the calcite interbed distribute localized. 5. The author synthetically analyzed the influence action of the micro-heterogeneity, the macro-heterogeneity and the structure heterogeneity to the oilfield water flood development. The efficiency of water flood is well in tiny structure of convex type or even type at top and bottom in which the water breakthrough of oil well is soon at the high part of structure when inject at the low part of structure, and the efficiency of water flood is poor in tiny structure of concave type at top and bottom. The remaining oil was controlled by sedimentary facies; the water flooding efficiency is well in the border or channel bar and is bad in the floodplain or the levee. The separation and inter layer have a little influence to the non-obvious positive rhythm reservoir, in which the remaining oil commonly locate within the 1-3 meter of the lower part of the separation and inter layer with lower water flooding efficiency.

珠江口盆地流花11-1油田储层非均质性及流体分布规律研究

In order to discover the distribution law of the remaining oil, the paper focuses on the quantitative characterization of the reservoir heterogeneity and the distribution law of the fluid barrier and interbed, based on fine geological study of the reservoir in Liuhuall-1 oil field. The refined quantitative reservoir geological model has been established by means of the study of core analysis, logging evaluation on vertical well and parallel well, and seismic interpretation and prediction. Utilizing a comprehensive technology combining dynamic data with static data, the distribution characteristics, formation condition and controlling factors of remaining oil in Liuhuall-1 oil field have been illustrated. The study plays an important role in the enrichment regions of the remaining oil and gives scientific direction for the next development of the remaining oil. Several achievements have been obtained as follows: l.On the basis of the study of reservoir division and correlation,eight lithohorizons (layer A, B_1, B_2, B_3, C, D, E, and F) from the top to the bottom of the reservoir are discriminated. The reef facies is subdivided into reef-core facies, fore-reef facies and backreef facies. These three subfacies are further subdivided into five microfacies: coral algal limestone, coralgal micrite, coral algal clastic limestone, bioclastic limestone and foraminiferal limestone. In order to illustrate the distribution law of remaining oil in high watercut period, the stratigraphic structure model and sedimentary model are reconstructed. 2.1n order to research intra-layer, inter-layer and plane reservoir heterogeneity, a new method to characterize reservoir heterogeneity by using IRH (Index of Reservoir Heterogeneity) is introduced. The result indicates that reservoir heterogeneity is medium in layer B_1 and B_3, hard in layer A, B_2, C, E, poor in layer D. 3.Based on the study of the distribution law of fluid barrier and interbed, the effect of fluid battier and interbed on fluid seepage is revealed. Fluid barrier and interbed is abundant in layer A, which control the distribution of crude oil in reservoir. Fluid barrier and interbed is abundant relatively in layer B_2,C and E, which control the spill movement of the bottom water. Layer B_1, B_3 and D tend to be waterflooded due to fluid barrier and interbed is poor. 4.Based on the analysis of reservoir heterogeneity, fluid barrier and interbed and the distribution of bottom water, four contributing regions are discovered. The main lies on the north of well LH11-1A. Two minors lie on the east of well LH11-1-3 and between well LH11-1-3 and well LH11-1-5. The last one lies in layer E in which the interbed is discontinuous. 5.The parameters of reservoir and fluid are obtained recurring to core analysis, logging evaluation on vertical well and parallel well, and seismic interpretation and prediction. Theses parameters provide data for the quantitative characterization of the reservoir heterogeneity and the distribution law of the fluid barrier and interbed. 6.1n the paper, an integrated method about the distribution prediction of remaining oil is put forward on basis of refined reservoir geological model and reservoir numerical simulation. The precision in history match and prediction of remaining oil is improved greatly. The integrated study embodies latest trend in this research field. 7.It is shown that the enrichment of the remaining oil with high watercut in Liuhua 11-1 oil field is influenced by reservoir heterogeneity, fluid barrier and interbed, sealing property of fault, driving manner of bottom water and exploitation manner of parallel well. 8.Using microfacies, IRH, reservoir structure, effective thickness, physical property of reservoir, distribution of fluid barrier and interbed, the analysis of oil and water movement and production data, twelve new sidetracked holes are proposed and demonstrated. The result is favorable to instruct oil field development and have gotten a good effect.

滑坡和泥石流致灾条件叠合技术及风险分析-以川藏公路然乌至鲁朗段为例

As one part of national road No. 318, Sichuan-Tibet (Chengdu-Lasha) Highway is one of traffic life lines connecting Tibet municipality to the inland, which is very important to the economic development of Tibet. In addition, it is still an important national defence routeway, with extremely important strategic position on maintaining the stability and solidarity of Tibet municipality and consolidating national defence. Particular geological condition, terrain and landform condition and hydrometeorological condition induce large-scale debris flows and landslides (including landslips) and the like geological hazards frequently occur along the highway. High frequency geological hazards not only result in high casualties and a great property loss, but also block traffic at every turn, obstructing the Sichuan-Tibet highway seriously. On the basis of considerable engineering geological investigation and analysis to the relative studying achievements of predecessors, it is found that one of the dominating reason incurring landslides or debris flows again and again in a place is that abundant loose materials are accumulated in valleys and slopes along the highway. Taking landslides' and debris flows along Ranwu-Lulang section of Sichuan-Tibet highway as studying objects, the sources and cause of formation of loose accumulation materials in the studying area are analyzed in detail, the major hazard-inducing conditions, hazard, dynamic risk, prediction of susceptibility degree of landslides and debris flows, and the relations between landslides and debris flows and various hazard-inducing conditions are systematically researched in this paper. All of these will provide scientific foundation for the future highway renovating and reducing and preventing geological hazards. For the purpose of quantitatively analyzing landslide and debris flow hazards, the conception of entropy and information entropy are extended, the conception of geological hazard entropy is brought forward, and relevant mathematics model is built. Additionally, a new approach for the dynamic risk analysis of landslide and debris flow is put forward based on the dynamic characteristics of the hazard of hazard-inducings and the vulnerability of hazard-bearings. The formation of landslide and debris flow is a non-linear process, which is synthetically affected by various factors, and whose formation mechanics is extremely complex. Aiming at this question, a muli-factors classifying and overlapping technique is brought forward on the basis of engineering geomechanics meta-synthesis (EGMS) thought and approach, and relevant mathematics model is also built to predict the susceptibility degree of landslide or debris flow. The example analysis result proves the validity of this thought and approach. To studying the problem that whether the formation and space distribution of landslides and debris flows are controlled by one or several hazard-inducing conditions, the theme graphics of landslides and debris flows hazard and various hazard-inducing conditions are overlapped to determine the relationship between hazard and hazard-inducing conditions. On this basis, the semi-quantitative engineering zonation of the studying area is carried out. In addition, the overlapping analysis method of the hazard-indue ing conditions of landslides and debris flows based on "digital graphics system" is advanced to orderly organize and effectively manage the spatial and attributive data of hazard and hazard-inducing conditions theme graphics, and to realize the effectively combination of graphics, images and figures.