Home – Home – Jornals – Contemporary Problems of Ecology 2006 number 3

The use of integral reflection and transmission operators to account for boundary conditions follows from the theory of wave scattering at curved interfaces in layered inhomogeneous media obtained in our previous studies. In this respect it appears practical to investigate the possibility of applying these operators instead of plane-wave reflection and transmission coefficients to describe high-frequency wavefields by Kirchhoff-type surface singular integrals common in forward and inverse seismic problems. We compare the integral reflection operator and reflection coefficient approaches in numerical experiments for single scattering of a spherical wave at a curved interface. We show that the use of the integral reflection operator allows eliminating artefacts associated with the use of the reflection coefficient and ensures including head waves.

A numerical-analytical solution for seismic and acoustic-gravity wave propagation is applied to a heterogeneous `Earth-atmosphere' model. Seismic wave propagation in elastic earth is described by a system of first-order dynamic equations of the elasticity theory and propagation of acoustic-gravity waves in the atmosphere by the linearized Navier-Stokes equations. The solution algorithm combines the time-domain integral Laguerre transform, the finite integral Bessel transform along the radial coordinate, and the FD solution of the reduced problem along the vertical coordinate. The suggested algorithm is tested with numerical experiments for the heterogeneous `Earth-atmosphere' model for different source locations.

Particle motion of homogeneous and inhomogeneous time-harmonic plane waves propagating in unbounded viscoelastic anisotropic media is generally elliptical. Exception is linear polarization of P and S waves propagating along some specific directions. A typical example is a linear polarization of SH waves propagating in a plane of symmetry of a viscoelastic anisotropic medium. Two most important characteristics of the particle motion are the orientation of the axes of the polarization ellipse and its eccentricity. They both usually vary considerably with the direction of wavefront propagation, and with varying strength of inhomogeneity of the considered plane wave. The orientation of the P -wave polarization ellipse generally differs from the direction of wavefront propagation, and it is usually closer to the direction of the energy flux. The orientation of the polarization ellipses of S waves often differs from the direction perpendicular to the wavefront propagation, and it is usually closer to the direction perpendicular to the direction of the energy flux. The eccentricity of the polarization ellipse depends particularly strongly on the inhomogeneity of the plane wave. For homogeneous plane waves, the particle motion is usually nearly linear, i.e., polarization ellipses have large eccentricity, and the eccentricity decreases with increasing inhomogeneity of the wave. For strongly inhomogeneous plane waves, the polarization ellipse becomes nearly circular, eccentricity being very small. The eccentricity of the polarization ellipse usually also decreases in the vicinity of singular directions.

A new method for estimating the quality of geoelectrical data inversion implies formalization of the concepts for horizontal and vertical resolution and equivalence. Various aspects of the suggested approach are illustrated by examples.

Amplitude versus offset (AVO) or amplitude versus angle (AVA) curves are nowadays regularly extracted from seismic data for various purposes of reservoir studies and characterization. After adequate processing, these curves represent the variation of the reflection coefficients with respect of offset or angle at points of key target reflectors, such as, for example, the top of a reservoir. Besides the information that can be obtained merely from their shape, AVO/AVA curves can also be used to invert more quantitative attributes, such as the intercept and the gradient of the reflection coefficient or, even better, the elastic-parameter contrasts. For a common-midpoint gather, the AVA curve is generally derived from its AVO counterpart by means of a well-known expression that relates the reflection angle to offset. It is to be noted that a successful inversion of the sought-for attributes is strongly dependent on the approximations of the reflection coefficient that are considered.
The recently introduced reflection impedance concept provides an attractive approximation of the elastic PP-reflection coefficient as a function of the ray parameter. In this sense, that approximation can be of value when amplitude versus ray parameter (AVP) curves are available from seismic data. It is to be noted that an AVP curve tend to be more reliable that its AVA counterpart. This is because the ray parameter, as a function of offset, depends on the RMS-velocity only, as opposed to the incidence angle, which also depends on the (more unstable) interval velocity.
In this paper, we propose an algorithm to invert elastic-parameter contrasts from AVP curves using the reflection impedance approximation of the PP-reflection coefficient. First results shown on synthetic data indicate that the procedure may offer a promising alternative to existing methods of inverting reservoir attributes from AVO/AVA curves.

For a horizontally layered medium with isotropic layers with constant velocity gradient, it is possible to estimate the velocity function (gradient and velocity at the top of the layer) and thickness of each layer. From large-offset PP seismic reflections one can estimate three traveltime parameters: the zero-offset two-way traveltime, the NMO velocity and a heterogeneity coefficient, using the shifted hyperbola approximation or a fractional approximation. From the estimated traveltime parameters at the top and bottom of a layer, it is possible to compute the thickness and velocity function of the layer. It is necessary to solve a simple nonlinear equation which also can be solved approximately. There are two solutions, corresponding to a positive and negative velocity gradient. Therefore, the sign of the velocity gradient must be chosen. When there is one solution, the velocity gradient in the layer is zero, and the result is the standard Dix equations.

The paper addresses traveltime processing for 2D models with laterally inhomogeneous layers and curved interfaces in shallow subsurface. Analytical equations are derived to approximate the relation between shallow velocity anomalies and NMO velocity as a basis for traveltime inversion to obtain an initial-approximation velocity model. The new approach is applied to real field data from regions with shallow velocity anomalies.

A model composed of silicate sand, crushed granite, and cement was subjected to biaxial compression, and the resulting deformations have been described. During the experiment lasting for several months, with a quasi-static level of applied stresses, water was repeatedly injected into the model. It has been established that when the volume of water injected into a zone of active fractures is small as compared with the volume of the model, the acoustic emission drastically increases. The shape of seismograms and curves of repeatedness of the acoustic signals recorded before and after water injection do not differ significantly, which implies that the effect is of trigger nature. The time succession of injection-caused acoustic events is qualitatively similar to the swarm and, in some case, aftershock activity of earthquakes. Hence, this factor may have an effect on seismicity. Different types of time succession of induced acoustic activity are described by an equation known from the kinetic theory of strength provided that the parameters of activation energy and working stresses depend on time.

Upper mantle anisotropy in Asia was investigated by inversion of Love and Rayleigh wave group velocity dispersion curves obtained by the FTAN procedure at periods between 10 and 150 s along paths that traverse Central Asia and Siberia. Our data were supplemented with data from the Center for Imaging the Earth's Interior (Boulder, Colorado), and 2500 to 4000 paths, depending on period, were used altogether. Locally averaged group velocity dispersion curves for Love and Rayleigh waves obtained by 2D seismic tomography for different periods were inverted to velocity-depth profiles, jointly and separately for SV (Rayleigh) and SH (Love) waves. SH waves have higher velocities than SV waves in the crust and upper mantle to a depth of ~300 km. The velocity profiles were used to estimate the mean anisotropy coefficient over depth intervals from the Moho to 100 km, from 100 to 200 km, and from 200 to 300 km. Anisotropy is the most prominent between 100 and 200 km, i.e., in the asthenosphere, and in anticorrelation with asthenospheric velocity: the lower the velocity the higher the anisotropy coefficient. Anisotropy is the greatest beneath active deforming regions of orogeny and vanishing in stable cratonic areas.

Several long-range profiles were shot with nuclear explosions within the Siberian Platform. The seismic records of both chemical- and nuclear-explosion profiles have revealed new peculiarities of the upper-mantle structure in the region. Seismic sections were constructed by ray tracing. To build the most informative starting model, the τ   (V,  x) method was applied, which implies construction of τ-time sections using the refraction and wide-angle reflection travel times for several velocity lines V = const. The modeling has revealed abnormally high (up to 8.4-8.5 km/s) upper-mantle velocities beneath the ancient Siberian Platform and low ones beneath the young West Siberian Plate. The upper mantle is shown to be of layered structure, with reflecting boundaries at depths of about 100, 150, 240, and 320 km. The boundary N1 at depths of 70-120 km is the most persistent; it often underlies low-velocity zones. The existence of velocity inversion zones and reflecting boundaries might be accounted for by fluid concentration at particular depths of the lithosphere, which suggests its rheological layering. The asthenosphere-lithosphere boundary has not been recognized in the velocity cross-sections.

A universal theory that accounts for the structure of fractured crust is suggested using data on strength and finite failure of rocks in high-pressure and high-temperature triaxial deformation tests. Correlation of the experimental results to the seismic cross section of the crust and mathematical modeling data indicates that the conditions at the Moho correspond to closure of the cracks network and annihilation of its hydraulic permeability. Therefore, upper mantle is rather dry than water-saturated, which means that the existing views of rock physics need a revision. Variations in crustal thickness are controlled by the trend of the geotherm and, besides, by olivine-to-serpentinite transition by the Hess reaction in oceanic crust. The Conrad discontinuity fits the horizontal foot of listric faults, waveguides (low-velocity zones) correspond to the cracked-porous inclusions and the lower crust to cataclastic state of rocks. The crust cut into fault-bounded blocks is modeled in the context of rotation tectonic waves that follow the sin-Gordon equation.

The variation principle was applied to the bifurcated deformation of an inelastic body to clear up under what conditions a regular system of shear and compaction bands can form. It has been established that a regular system of bands where plastic deformations are localized should form most likely in materials with a distinct yield point. Internal friction and dilatancy promoting the development of instability of the material in the regime of strain-hardening contribute much to this phenomenon.

Results of laboratory experiments and observation data on postseismic deformations in different regions were used to detect the regularities of formation of slow relative block displacements provoked by dynamic events. It is shown that under a gradual change of the stress-strain state of a rock massif, the asymptotic form of displacement with time is close to the law of a quasi-static flow, and under a drastic change in deformation rate, the relative block displacements are controlled mainly by the laws of friction force change during shear. The results obtained show that the dynamics of forces resisting to shears along block boundaries must be taken into account when constructing geomechanical models of different scales.

The paper presents an approach implying constructing continuum images of liquid- or gas-saturated porous and cracked solids. The hypothesis of classical continuum mechanics that closely positioned points of a solid occur in similar physical conditions breaks down in the case of porous media in which fluid and matrix material have strongly different properties. We suggest to represent a real microheterogeneous body by its continuum image using continuum operators, whereby the conservation laws are applied rather to the image than to its prototype.
The continuum image is described by the standard equation of motion and the respective real solid by infinite-order equations. The latter include velocities of elastic waves as well as infinitely low velocities of anomalies due to a great number of the degrees of freedom in microheterogeneous media consisting of separate elementary blocks.

We discuss a mechanism of nonhysteretic strain-amplitude dependent dissipation of elastic waves in microinhomogeneous media containing linearly dissipative and nonlinear elastic soft defects. The combined effect of these factors can cause well-pronounced essentially dissipative nonlinearity of elastic waves. The mechanism presumably works in rocks and other materials with microstructure where cracks and grain contacts act as soft inclusions. Unlike homogeneous materials where strain dependent variations of dissipation and elasticity are commonly of the same order, dissipation variations in microinhomogeneous materials can be times as great as the respective elasticity variations. A slight (a few percent) change in Young's modulus can correspond to a manifold dissipation change at moderate strain typical of acoustic and seismic applications. Moreover, in the case of variable frequency, this nonlinear mechanism can increase or decrease Young's modulus or reduce its variation to zero at certain frequencies while dissipation variations remain quite large. The nonhysteretic mechanism may act concurrently with hysteretic mechanisms, which are due to friction and adhesion effects and are commonly invoked for materials with microstructure, and influence the relationship between strain dependent dissipation and elasticity observed in experiments on nonlinear dissipative effects.

A method for finding the mutual correspondence between stereopair frames and geometrical transforming a videoimage sequence is presented. The method is adapted to parallelization and implementation on a cluster computer system. Two schemes for parallelization of the videoframe transformation process are considered. Efficiency of the proposed parallelization algorithms is analyzed a priori.

Models based on Lattice Gas Cellular Automata (LGCA) can be considered as an alternative to the traditional approach to spatial dynamics simulation. The possibilities of using LGCA automata to simulate a 2D wave are investigated. Results of some experiments (simulation of moving boundaries, effects of quiescent particles on the wave propagation speed, and wave simulation in inhomogeneous media) are discussed.

the problem of cellular automaton simulating fluid flows in porous media is considered. A 3D cellular automaton model RD-1 was used. the model structure is briefly described. The boundary conditions and the time averaging procedure are considered in detail. A series of computational experiments concerned with simulating a fluid flow through a porous medium is presented. The flow characteristics are compared with the Darcy law and other physical laws.Tthe velocity field of a model flow in a porous medium is illustrated.

An algorithm for 3D gyrodynamics simulation of a dust cloud by the particle-in-cell method is considered. Approaches to parallelization of the methods are presented. Three versions of parallel realizations are compared. The best parallel program version gave a 7-times speedup on 32 processors.

A neural network approach to constructing finite-difference adaptive meshes is proposed. It is based on the learning algorithm for Kohonen self-organizing maps. An algorithm for constructing meshes illustrated by an example of 2D domains on a plane is presented. It is shown that the proposed approach allows parallel mesh construction with a more than 90 % efficiency.

An experimental Parallel Library for Vertical Image Processing (PLVIP) based on vertical data representation is considered. The main characteristics of the library (data format, computing process management, and implemented programs) and examples of its application are presented.

A parallel method for solving the Poisson equation as a result of using the fast Fourier transform is constructed. The main feature of the method is increasing speedup of a parallel program due to applying a multigrid method to calculating the zero potential harmonic. A multigrid method for cylindrical geometry is implemented.

A parallel substitution algorithm was the base for developing a simple procedure of constructing a self-reproducing cellular structure of arbitrary size.

A method for replacing the problem of nodes of the percolation theory by the problem of links is described. It is hypothesized that the correlations arising between the states of grid elements do not change the values of critical percolation probabilities. the hypothesis is partially substantiated.

A definition of a parallel automaton, that is a P-automaton, is given. Notions of atomic, partial, and complete states of a P-automaton, a parallel input and output are introduced. Parallel functions of transitions and outputs, which are given by transition and output matrices, are defined on these sets. Conditions to which the matrices should correspond in order to ensure a correct and unique P-automation description are formulated. Three modes of its functioning, namely, synchronous, ordinary, and mixed, are considered. Statements about the conditions of reducing a P-automaton to a parallel-sequential composition of sequential automata and on sufficiency of the ordinary P-automaton model are formulated.

A new approach to integration of computing clusters for higher total performance is presented. the limits arising from local addressing of nodes in the clusters are overcome. the architecture of inter-cluster communication service, the functions of its components, and a model of application management based on the Service are considered. the dependence of the number of inter-cluster communication channels on the number of clusters and processes of a parallel application is shown. Results of performance tests of the Service are illustrated by a model and real systems of two closed clusters. Efficiency of utilizing the Service is shown in comparison with utilizing low-level functions of the protocol TCP/IP.

Problems of optimizing the structure and organizing the exchanges in circulant networks used as communication networks of parallel computer systems are considered and results obtained in this field are reviewed. Triple analytically defined circulant networks of any diameter with the maximum number of vertices are obtained. An effective dynamic two-terminal routing algorithm and a broadcasting algorithm that minimizes the time of execution and the load of messages in the network are developed for them.

Embedding of a diagnostic circulant graph into a graph of computing system intermodule connections is considered. the host graph is a graph of group of direct product of cyclic subgroups. The embedding is bases on partitioning the group into cossets via a cyclic subgroup of the chosen element degrees. Dependence of the introduced robustness index of the embedded circulant on the relative position of initial vertices in cossets of group is studied. It is found that the structural robustness index value depends substantially on the mentioned factor. Embedding parameters for which the embedded graph structural robustness does not depend on the relative position of initial vertices in its cossets are found.

The paper explores the concept of efficiency of the constitutional control, the criteria and evaluations of its efficiency. The author analyses the problems of researching the efficiency of the constitutional control and execution of the Constitutional court decisions. The author also provides the justification of the value characteristics of the criterial approach to the efficiency of constitutional control.

The aim of the paper is to define the state of legal education under the conditions of radical transformations of the state and the economy and to single out the main factors hindering its development. The expansion of legal education, change of its content and forms is the result of fundamental social changes. The strategy of this process is due to objective factors: the systemic crisis in the country, the task of forming of a new community of lawyers, the gap between the real functioning of educational institutions and their ideal models, the low level of legal education. The article is based on the analysis of the legislature on legal education and the practice of its application, personal experience of teaching legal subjects.

It is noted that alongside legitimate interests which have substantive character there can be specific legitimate interests. These interests are of a procedural character which is more clearly manifested in the course of civil proceedings. By analyzing civil procedural rules and subject composition of civil procedural legal relationships the author studies and classifies procedural legitimate interests and offers their definitions.

The paper contains a comparative analysis of the rights to things and intellectual property. Since the Roman law, the civil rights are subdivided into rights to things and obligations. However, this dichotomy does not correspond to the modern development of the law, in which the division of rights into absolute and relative becomes more important. Absolute rights are subdivided into personal non-property rights, rights to things (concerning things

The author aims to demonstrate the need in a special method of securing a preliminary contract. Concluding a preliminary contract the parties provide for transference of money to secure the performance of the contract. Usually this money is termed a deposit or an advance. But these categories can

The article argues for the necessity of transition from the statement "the law gives the owner freedom to act" to the status of the law as a mechanism of restrictions defining the limits of freedom in property relations. The author proposes to consider the limits of the owner

A continuous-discrete model of time course of a phytocenosis is considered. The properties of the model are studied, and particularly for the case of two-age cenosis it is demonstrated that cyclic rhythms of only length 2 can be realized. Private cases of the model with a large number of age classes are studied numerically, and it is shown that in the cases under study cyclic regimes of various lengths, including chaotic regimes, are possible.

A new method for calculation of paleoclimatic parameters by palinological data based on geometrical interpretation of relations between the elements of spore-pollen spectra and climatic systems is proposed. Examples of reconstruction of paleoclimate in Holocene are presented.

A series of parameters used in modern ecology as indicators of state or "health" of ecological systems is considered. Among the goal functions, exergy is chosen as a theoretically well substantiated (in ther

An attempt is made to compare local faunas of the Northern Hemosphere on the example of the Arctiinae subfamily of lepidoptera at the level of genera with the help of cluster an

An experimental and theoretical study of a closed aquatic ecosystem "producer - consumer" (Chlorella vulgaris - Paramecium caudatum) has been cattied out. The experiment was made on an installation with spatially divided units which consisted of a reactor for Chlorella cultivation (the "producer" unit) and a reactor for paramecia (the " consumer" unit). A direct relation between the paramecia reproduction with chlorella consumption and nitrogen release in ammonium form which is the most preferable for the growth of chlorella has been found. As a result of the theoretical study of the aquatic biotic "producer - consumer" cycle with spatially divided units, a contribution of paramecia to the nitrogen circulation is confirmed. It is demonstrated that simultaneously with the increase of concentration of nitrogen released due to metabolism of paramecia, an increase of productivity of microalgae takes place.

Peat deposits of the hydromorphic swamp complex on the poorly explored right-hand bank of the Yenisei river in the subzone of northern taiga have been studied. A brief geobotanic characteristic of the vegetable cover of the swamp massif, stratigraphy and some general technical properties of peat deposits are presented. It is found that the changes of humidity, degree of decomposition and ash content of peat have an oscillatory character, which is due to re-deposition and multiple mixture of peat layers under the influence of thermokarst processes. A conclusion is made of degradation of a flat-frost mound bog peat deposits and their evolution to pool hummock ridge bog complexes.