-
Date submitted2023-09-05
-
Date accepted2024-11-07
-
Date published2025-02-25
Investigation of the accuracy of constructing digital elevation models of technogenic massifs based on satellite coordinate determinations
At all stages of the life cycle of buildings and structures, geodetic support is provided by electronic measuring instruments – a laser scanning system, unmanned aerial vehicles, and satellite equipment. In this context, a set of geospatial data is obtained that can be presented as a digital model. The relevance of this work is practical recommendations for constructing a local quasigeoid model and a digital elevation model (DEM) of a certain accuracy. A local quasigeoid model and a DEM were selected as the study objects. It is noted that a DEM is often produced for vast areas, and, therefore, it is necessary to build a local quasigeoid model for such models. The task of assessing the accuracy of constructing such models is considered; its solution will allow obtaining a better approximation to real data on preassigned sets of field materials. A general algorithm for creating both DEM and local quasigeoid models in the Golden Software Surfer is presented. The constructions were accomplished using spatial interpolation methods. When building a local quasigeoid model for an area project, the following methods were used: triangulation with linear interpolation (the least value of the root mean square error (RMSE) of interpolation was 0.003 m) and kriging (0.003 m). The least RMSE value for determining the heights by control points for an area project was obtained using the natural neighbour (0.004 m) and kriging (0.004 m) methods. To construct a local quasigeoid model for a linear project, the following methods were applied: kriging (0.006 m) and triangulation with linear interpolation (0.006 m). Construction of the digital elevation model resulted in the least aggregate value of the estimated parameters: on a flat plot of the earth’s surface – the natural neighbour method, for a mountainous plot with anthropogenic topography – the quadric kriging method, for a mountainous plot – quadric kriging.
-
Date submitted2023-10-31
-
Date accepted2024-11-07
-
Date published2025-01-24
Specifics of magnetotelluric studies in Antarctica
One of the priority areas of scientific research in Antarctica is the study of its deep structure. Most of the continent is covered with a thick ice sheet, so the main geoscientific data are acquired using geophysical methods, among which magnetotelluric (MT) ones have the greatest penetration depth and insignificant environmental impact. The possibility of acquiring high-quality MT data in the conditions of the sixth continent has long been questioned. The work is aimed at studying the specifics of magnetotelluric survey in Antarctica. The following tasks were set: to summarize the world experience of studying Antarctica using MT sounding methods; to identify factors that negatively affect the high-quality data acquisition; to determine methods for minimizing the influence of these factors. The article analyses geophysical studies conducted by the magnetotelluric sounding method in the Antarctic region from 1964 to the present. The application of the method is complicated by the following: extremely low temperature affects the drop in the batteries capacity, freezing of the non-polarizing electrodes solution, and changes in the strength properties of materials. Electromagnetic noise occurs during strong winds; proximity to the magnetotelluric field source can violate the plane wave principle on which the method is based. The ice sheet covering most of Antarctica does not allow acquiring optimal values of the contact resistance of the electrode grounding; the extended coastline distorts the acquired data. Studies of the influence of factors complicating the MT sounding method in the coastal and central parts of Antarctica made it possible to formulate recommendations for preparing equipment and adapting the work procedure, modifying the processing flow and a set of measures to ensure safety, the implementation of which will both allow safe performance of geophysical investigations and high-quality data acquisition.
-
Date submitted2024-04-11
-
Date accepted2024-06-03
-
Date published2024-07-04
Assessment of the efficiency of acid mine drainage purification (using the example of copper-pyrite mines in the Middle Urals)
According to the results of the anti-rating of regions with extreme pollution of watercourses in the Sverdlovsk region, the largest number of polluted rivers has been recorded in recent years – more than a quarter of all high and extremely high pollution. One of the sources of pollution of natural water bodies in the Middle Urals are closed and flooded copper-pyrite mines, where acidic mine drainage continue to form and unload to the surface. Several of them have organized collection and a two-stage acidic drainage purification system, including neutralization with lime milk and settling in a clarifier pond. Despite the identical schemes, different indicators of pollutants are recorded during discharge into water bodies. The aim of the work is to evaluate the effectiveness of the applied acid mine drainage purification system and identify the parameters affecting the quality of treated mine water. Laboratory studies were performed using methods of flame emission spectrometry, flame atomic absorption, atomic absorption spectrometry, mass spectrometry with ionization in inductively coupled plasma, potentiometric, etc. It has been established that the existing mine drainage purification system at the Degtyarskii mine makes it possible to significantly reduce the concentrations of most toxic components of mine waters to almost standard values. At the Levikhinskii mine, the multiplicity of exceeding the maximum permissible concentrations reaches hundreds and thousands of times. To achieve a higher degree of purification, it is necessary that the duration of passive purification is sufficient for the reactant to interact with acidic waters. However, to ensure this possibility, it will require the creation of a cascade of ponds with an area of several thousand hectares. If the current two-stage system is quite effective for the Degtyarskii mine, then for Levikhinskii it is necessary to switch to the use of more modern systems, including three stages of purification.
-
Date submitted2021-10-27
-
Date accepted2023-06-20
-
Date published2023-12-25
Geomechanical analysis of the impact of the new tunnels construction in the vicinity of existing underground subway structures on the state of the soil massif
The specificity of the behavior of the soil massif near the tunnel under construction in difficult mining conditions is considered. It was revealed through the joint interpretation of the results of geophysical measurements in the tunnel and computer simulation. The results of field geophysical studies to identify areas of decompacted soil behind the lining in two existing tunnels during successive drilling of two new tunnels under them are described. A method to analyze the response of the lining to impact was used providing for the calculation of its energy. It has been established that the decompaction zones are mainly located in the lateral lower areas of the tunnel. To substantiate the mechanism of formation of cavities, computer simulations were carried out using the finite element method with the COMSOL Multiphysics software. The finite element model is built on the Drucker – Prager criterion in the variant of a two-dimensional problem statement. It is shown that at the initial position of two old tunnels, the areas of decompaction can develop mainly on the sides. The position of the zones changes significantly when excavating two new tunnels. Soil decompaction zones appear between the tunnels and there is a tendency for the areas to spread to the upper point of the tunnel. According to geophysical data time delays in the impact of new tunnels on the existing line are noted, as well as a decrease in the size of decompacted soil areas over time. There is a satisfactory agreement between the positions of the decompaction areas and voids obtained by the geophysical method and the results of numerical simulation.
-
Date submitted2022-08-01
-
Date accepted2022-11-17
-
Date published2023-02-27
Use of machine learning technology to model the distribution of lithotypes in the Permo-Carboniferous oil deposit of the Usinskoye field
- Authors:
- Denis V. Potekhin
- Sergei V. Galkin
Permo-Carboniferous oil deposit of the Usinskoye field is characterized by an extremely complex type of the void space with intense cross-sectional distribution of cavernous and fractured rock. In this study, for this production site, the process of 3D geological modeling has been implemented. At the first stage, it provided for automated identification of reservoir volumes by comparing the data of core and well logging surveys; at the second stage, identification of rock lithotypes according to Dunham classification is performed on the basis of comparison of thin sections examination and well logging data. A large array of factual information enables the use of machine learning technology on the basis of Levenberg – Marquardt neural network apparatus toward achievement of our research goals. The prediction algorithms of reservoir and rock lithotype identification using well logging methods obtained on the basis of the training samples are applied to the wells without core sampling. The implemented approach enabled complementing the 3D geological model with information about rock permeability and porosity, taking into account the structural features of the identified lithotypes. For the Permo-Carboniferous oil deposit of the Usinskoye field, the volumetric zoning of the distribution of different rock lithotypes has been established. Taking into account the lithotypes identified based on machine learning algorithms, density and openness of fractures were determined, and fracture permeability in the deposit volume was calculated. In general, during the implementation, the machine learning errors remained within 3-5 %, which suggests reliability of the obtained predictive solutions. The results of the research are incorporated in the existing 3D digital geological and process model of the deposit under study.
-
Date submitted2022-05-26
-
Date accepted2022-11-17
-
Date published2022-12-29
Renovation method of restoring well productivity using wavefields
A stagewise theoretical substantiation of the renovation vibrowave method of influencing the near-wellbore zone of reservoir for restoring well productivity is presented. The area of treatment by the proposed method covers the reservoir with a heterogeneous permeability with fractures formed by fracking. In this method a decrease in concentration of colmatants occurs due to a change in direction of contaminants migration. Under the influence of pressure pulses, they move deep into the reservoir and disperse through the proppant pack. The results of mathematical modelling of the propagation of pressure wave and velocity wave and the calculations of particles entrainment in wave motion are presented.
-
Date submitted2021-12-21
-
Date accepted2022-06-20
-
Date published2022-11-10
Scientific justification of the perforation methods for Famennian deposits in the southeast of the Perm Region based on geomechanical modelling
The article presents the results of analysing geological structure of the Famennian deposits (Devonian) in the Perm Region. Numerical modelling of the distribution of inhomogeneous stress field near the well was performed for the two considered types of perforation. With regard for the geometry of the forming perforation channels, numerical finite element models of near-wellbore zones were created considering slotted and cumulative perforation. It is ascertained that in the course of slotted perforation, conditions are created for a significant restoration of effective stresses and, as a result, restoration of reservoir rock permeability. Stress recovery area lies near the well within a radius equal to the length of the slots, and depends on the drawdown, with its increase, the area decreases. From the assessment of failure areas, it was found that in case of slotted perforation, the reservoir in near-wellbore zone remains stable, and failure zones can appear only at drawdowns of 10 MPa and more. The opposite situation was recorded for cumulative perforation; failure zones near the holes appear even at a drawdown of 2 MPa. In general, the analysis of results of numerical simulation of the stress state for two simulated types of perforation suggests that slotted perforation is more efficient than cumulative perforation. At the same time, the final conclusion could be drawn after determining the patterns of changes in permeability of the considered rocks under the influence of changing effective stresses and performing calculations of well flow rates after making the considered types of perforation channels.
-
Date submitted2022-03-17
-
Date accepted2022-10-04
-
Date published2022-11-10
Improving the reliability of 3D modelling of a landslide slope based on engineering geophysics data
Landslides are among the most dangerous geological processes, posing a threat to all engineering structures. In order to assess the stability of slopes, complex engineering surveys are used, the results of which are necessary to perform computations of the stability of soil masses and assess the risks of landslide development. The results of integ-rated geological and geophysical studies of a typical landslide slope in the North-Western Caucasus spurs, composed of clayey soils, are presented. The purpose of the work is to increase the reliability of assessing the stability of a landslide mass by constructing a 3D model of the slope, including its main structural elements, identified using modern methods of engineering geophysics. Accounting for geophysical data in the formation of the computed 3D model of the slope made it possible to identify important structural elements of the landslide, which significantly affected the correct computation of its stability.
-
Date submitted2021-03-31
-
Date accepted2021-09-29
-
Date published2021-10-21
Methodology of modeling nonlinear geomechanical processes in blocky and layered rock masses on models made of equivalent materials
- Authors:
- Boris Yu. Zuev
The research purpose is to develop a methodology that increases the reliability of reproduction and research on models made of equivalent materials of complex nonlinear processes of deformation and destruction of structured rock masses under the influence of underground mining operations to provide a more accurate prediction of the occurrence of dangerous phenomena and assessment of their consequences. New approaches to similarity criterion based on the fundamental laws of thermodynamics; new types of equivalent materials that meet these criteria; systems for the formation of various initial and boundary conditions regulated by specially developed computer programs; new technical means for more reliable determination of stresses in models; new methods for solving inverse geomechanical problems in the absence of the necessary initial field data have been developed. Using the developed methodology, a number of complex nonlinear problems have been solved related to estimates of the oscillatory nature of changes in the bearing pressure during dynamic roof collapse processes; ranges of changes in the frequency of processes during deformation and destruction of rock mass elements, ranges of changes in their accelerations; parameters of shifts with a violation of the continuity of the rock mass under the influence of mining: secant cracks, delaminations, gaping voids, accounting for which is necessary to assess the danger of the formation of continuous water supply canals in the water-protection layer.
-
Date submitted2021-03-16
-
Date accepted2021-07-27
-
Date published2021-10-21
Allocation of a deep-lying brine aquifer in the rocks of a chemogenic section based on the data of geophysical well logging and 2D seismic exploration
Advancement in the production of potassium fertilizers is an important strategic task of Russian agricultural industry. Given annually growing production rates, the reserves of discovered potassium-magnesium salt deposits are noticeably decreasing, which creates the need to ensure stable replenishment of the resource base through both the discovery of new deposits and the exploitation of deep-lying production horizons of the deposits that are already under development. In most cases, deposits of potassium-magnesium salts are developed by underground mining. The main problem for any salt deposit is water. Dry salt workings do not require any additional reinforcement and can easily withstand rock pressure, but with an inflow of water they begin to collapse intensively – hence, special attention is paid to mine waterproofing. Determination of spatial location, physical and mechanical properties of the aquifer and water-blocking stratum in the geological section represent an important stage in the exploration of a salt deposit. The results of these studies allow to validate an optimal system of deposit development that will minimize environmental and economic risks. On the territory of Russia, there is a deposit of potassium-magnesium salts with a unique geological structure – its production horizon lies at a considerable depth and is capped by a regional aquifer, which imposes significant limitations on the development process. To estimate parameters of the studied object, we analyzed the data from CDP seismic reflection survey and a suite of methods of radioactive and acoustic well logging, supplemented with high-frequency induction logging isoparametric sounding (VIKIZ) data. As a result of performed analysis, we identified location of the water-bearing stratum, estimated average thickness of the aquifers and possible water-blocking strata. Based on research results, we proposed methods for increasing operational reliability of the main shaft in the designed mine that will minimize the risks of water breakthrough into the mine shaft.
-
Date submitted2021-02-05
-
Date accepted2021-03-30
-
Date published2021-04-26
Assessment of negative infrastructural externalities when determining the land value
- Authors:
- Elena N. Bykova
The work forms and substantiates the concept of land value, based on a new institutional theory. The infrastructural component of the cost of land in the presented concept determines, on the one hand, the efficiency of the use of natural resources, properties, demand for land on the market, on the other hand, the costs, which are determined not only by capital investments in construction of engineering infrastructure, but also by losses associated with restrictions on activities within zones with special conditions for territory use, creation of unfavorable conditions for economic activity, small contours, irregularities and others on a specific land plot, which are external negative infrastructural externalities that create losses of rights holders of land plots that are not compensated by the market, falling within the boundaries of these zones. Methods for assessing the impact of such negative infrastructural externalities on the cost of land encumbered by zones in different conditions of land market activity have been developed and tested, based on an expert-analytical approach (depressed market); the ratio of market values of land plots encumbered and unencumbered by a specific zone, and qualimetric modeling (inactive market); modeling by introducing into the model the factor of presence of zones with special conditions for territory use, based on the grouping of zones according to similar regulations for use, or by introducing the parameters of this factor (active market). Methods for taking into account spatial deficiencies and compensating for restrictions and prohibitions on activities on the territory of land plots with an individual market assessment are proposed.
-
Date submitted2020-08-28
-
Date accepted2020-12-16
-
Date published2020-12-29
Metrological support of equipment for geophysical research
The article discusses the problems associated with metrological support of equipment for geophysical research, issues of ensuring the uniformity of well measurements, the creation of Russian standards for calibrating well equipment when determining the porosity coefficients and oil, gas, water saturation, determining the parameters of defects during well cementing and technical condition of casing strings and the water-oil-gas flow. The problems of creating full-fledged methods for measuring the parameters of oil and gas fields with conventional and hard-to-recover reserves have been investigated. The key directions of development of the Russian metrological support of geophysical well measurements were determined. The tasks that need to be solved to create metrological support for geophysical well logging as an industry that meet international standards are indicated. The expediency of creating a Russian Geophysical Center for Metrology and Certification, the need to develop a new and update the existing regulatory framework, which will allow Russian geophysics to reach the level of world leaders in the field of geophysical research, are substantiated.
-
Date submitted2019-06-22
-
Date accepted2019-09-11
-
Date published2020-04-24
Effective capacity building by empowerment teaching in the field of occupational safety and health management in mining
The paper is dealing with a developed concept named Empowerment Teaching, which is based on practical teaching experience gained in various mining universities. It is demonstrated that this concept can be used to increase the effectiveness of knowledge transfer to mining countries in the world, as well as to overcome cultural barriers between lecturers and their students. The two models of participatory training, which are proposed to be named “physical” and “emotional” models, are portrayed. The authors are convinced that participatory training methods can be an ideal answer to a challenge associated with workers’ competencies in mining, namely – the potential of highly motivated and well-educated young academics is often diminished by a lack of ability to apply their knowledge. A special emphasis is made on the possible application of empowerment teaching for educational and training activities in the field of occupational safety and health (OSH), which is a matter of utmost importance for the mining industry. Several benchmarking initiatives in the field of OSH (“safety culture”, zero-accident vision) are underlined to be encouraged and promoted by means of new teaching methods. The examples of successful international cooperation among universities are given, as well.
-
Date submitted2019-03-24
-
Date accepted2019-05-13
-
Date published2019-08-23
Calculation of Oil-saturated Sand Soils’ Heat Conductivity
- Authors:
- J. Sobota
- V. I. Malarev
- A. V. Kopteva
Nowadays, there are significant heavy high-viscosity oil reserves in the Russian Federation with oil recovery coefficient not higher than 0.25-0.29 even with applying modern and efficient methods of oil fields development. Thermal methods are the most promising out of the existing ways of development, main disadvantage of which is large material costs, leading to the significant rise in the cost of extracted oil. Thus, creating more efficient thermal methods and improving the existing ones, is the task of great importance in oil production. One of the promising trends in enhancing thermal methods of oil recovery is the development of bottomhole electric steam generators. Compared to the traditional methods of thermal-steam formation treatment, which involve steam injection from surface, well electrothermal devices can reduce energy losses and improve the quality of steam injected into the formation. For successful and efficient organization of oil production and rational development of high-viscosity oil fields using well electrothermal equipment, it is necessary to take into account the pattern of heat propagation, both in the reservoir and in the surrounding space, including the top and bottom. One of the main values characterizing this process is the heat conductivity λ of oil-bearing rocks. The article describes composition of typical oil-saturated sand soils, presents studies of heat and mass transfer in oil-saturated soils, reveals the effect of various parameters on the heat conductivity of a heterogeneous system, proposes a method for calculating the heat conductivity of oil-bearing soils by sequential reduction of a multicomponent system to a two-component system and proves the validity of the proposed approach by comparing acquired calculated dependencies and experimental data.
-
Date submitted2017-10-29
-
Date accepted2017-12-31
-
Date published2018-04-24
Risk assessment of accidents due to natural factors at the Pascuales – Cuenca multiple-use pipeline (Ecuador)
- Authors:
- Dzh. Zambrano
- S. V. Kovshov
- E. A. Lyubin
The natural aspects of the accident risk at the Pascuales – Cuenca multiple-use pipeline (Ecuador) are analysed in the paper. The Russian Methodological recommendations for the quantitative analysis of accident risks at hazardous production plants of oil trunk pipelines and oil product trunk pipelines issued in 2016 are used as a methodological framework due to relatively poorly defined evaluation mechanism for natural factors of accidents at oil trunk pipelines in the most widespread international accident risk assessment methodologies. The methodological recommendations were updated to meet the environmental conditions of oil pipelines of Latin America. It was found that the accidents due to natural factors make up approximately 15 % of cases at oil trunk pipelines in Ecuador. Natural geographical features of the areas surrounding the main Ecuadorian Pascuales–Cuenca oil trunk pipeline and its relatively short length allow defining three zones along the line in terms of the accident risk: lowland coastlines, high plateaus, and foothills. Calculations and analysis revealed that the maximum predicted specific frequency of accidents is characteristic of the lowland seaside area. The evidence showed that physical and chemical properties of soils and significant seismic activity are the root causes of failures.
-
Date submitted2014-12-19
-
Date accepted2015-02-27
-
Date published2015-12-25
Geological and mathematical model of secondary sulphidisation zone polarizability as a function of oil and gas deposit depth
- Authors:
- O. F. Putikov
- S. A. Ivanov
The differential-normalized method of electrical exploration (DNME), by means of polarization study of the secondary sulphidation zone (depth ~ 0,4-0,6 km), allows making conclusion about existence and parameters of an oil and gas deposit (depth up to ~ 2-5 km). The approximate solution of the system of nonlinear differential equations in partial derivatives for concentration of sulfur-hydrogen and pyrite (polarizability) as the function of oil-gas deposit depth has been received. The results are confirmed by experimental investigations in the North Sea.
-
Date submitted2013-07-11
-
Date accepted2013-09-05
-
Date published2014-03-17
Mathematical modelling of applied problems of rock mechanics and rock massifs
- Authors:
- A. P. Gospodarikov
- M. A. Zatsepin
The variety of the mining and geological conditions with further increasing in depth of the development of bedded deposits leads to necessity for the analysis of stress and strain state near different types of excavations.
-
Date submitted2010-07-28
-
Date accepted2010-09-07
-
Date published2011-03-21
Modelling cavities nearwell areas for studyind of possibilities neitron well methods
- Authors:
- N. A. Poznyakova
Look into the possibility neutron method geophysical logging for determination zones of cavity nearwell areas gas reservoirs and determine the filler cavities. Relevant practicability using these methods on underground gas storage. Give an example these research on model and real conditions.
-
Date submitted2009-10-18
-
Date accepted2009-12-09
-
Date published2010-09-22
Influence of blocky-hierarchy structure of geological medium on the forming of natural and technogenic geodynamic fields in the North-West part of the Amur geoblock
- Authors:
- I. Yu. Rasskazov
- B. G. Saksin
On the example of the South-Argun uranium ore district of Zabaikalie consideration is given to the interaction of regional and local fields of stresses. The district includes a number of deposits prone to rock bursts. The performed analysis contains data on recent geodynamics of the explored part of the Amur geoblock as well as on geomechanical and geologo-structural materials of Tulukuev caldera and of Antey deposit. It is shown that with the use of the M.A. Sadovsky relationship the indices of local structural forms of a row are subjected to unified laws of deformation and fracturing. Conclusion was made that at the beginning of deposit exploitation the direction of action of stress field within rock mass conformed to regional one, and then it considerably changed under influence of large-scale mining operations.
-
Date submitted2009-09-10
-
Date accepted2009-11-11
-
Date published2010-06-25
Stages of training activation methods development
- Authors:
- M. A. Perfilova
The article is dedicated to the problem of developing methods of training activation. Historical and pedagogical background for the development of these methods and distinguishes the stages of training activation methods development from the ancient times till nowadays are analysed. The brief chronological survey is presented tabular. The main features of the methods of active training are pointed out.
-
Date submitted2009-09-19
-
Date accepted2009-11-24
-
Date published2010-06-25
The forecast of the stress-strain state of massifs rocks in the bed deposits
- Authors:
- A. P. Gospodarikov
- M. A. Zatsepin
The variety of the mining and geological conditions with further increasing in depth of the development of bedded deposits leads to necessity for the analysis of stress-strain state near different types of excavations.
-
Date submitted2009-09-15
-
Date accepted2009-11-02
-
Date published2010-06-25
Mathematical modeling of stress-strain state of the mined seam deposits
- Authors:
- A. P. Gospodarikov
- M. A. Zatsepin
Development of bedded deposits is associated with man-caused distortion of specific environment – the rock massifs, which are very complicated in their composition, can vary significantly in mechanical properties and is characterized with a wide variety of laws and techniques to assess its stress-strain state.
-
Date submitted2009-08-29
-
Date accepted2009-10-30
-
Date published2010-02-01
The unit for gravity concentration of fine particles
- Authors:
- Ya. V. Kuskova
It was made new modification of device for allows increasing recovery, to increase productivity and to lower size of concentrated particles. The device has a deck executed as a circle, which divides into two sectors; each sector has circular rifles with progressively growing height from the centre to periphery of a deck. The table is equipped with the high-frequency step-by-step engine ensuring continuous rotation of a table and submission counter impulses for shear of particles in a direction, opposite to rotation of a table. The quantity of angle sectors can be various. Testing of prototype of a table has confirmed an opportunity of increase of recovery, decrease size of extracted dense minerals and increase of throughput.
-
Date submitted2009-07-22
-
Date accepted2009-09-22
-
Date published2010-04-22
Surveying methods for providing geodynamical safety of mining operations
- Authors:
- E. N. Kuznetsova
It is impossible to solve the main task of geodynamics without using technologies including traditional measurements of bench marks excursion, as well as new technologies of the end of XX – beginning of XXI century, such as photogrammetry, global positioning system, laser measurements and technologies of geophysical well logging. The paper concerns with the analysis of the mine surveying methods facilitating the solution of tasks related to the geodynamical safety assurance.
-
Date submitted2009-07-16
-
Date accepted2009-09-28
-
Date published2010-04-22
Efficient numerical methods for geotechnical problems
- Authors:
- M. A. Karasev
This article give a general overview of methods to solve large scale geotechnical problems, nonlinear high deformation and plastic problems, self contact of strata and application of high performance computing.