Study on the Dissociation Behavior of Gas Hydrate in Porous Sediment Based on Fractal Theory

The dissociation process of gas hydrate was regarded as a gas-solid reaction without solid production layer when the temperature was above the zero centigrade. Based on the shrinking core model and the fractal theory, a fractional dimension dynamical model for gas hydrate dissociation in porous sediment was established. The new approach of evaluating the fractal dimension of the porous media was also presented. The fractional dimension dynamical model for gas hydrate dissociation was examined with the previous experimental data of methane hydrate and carbon dioxide hydrate dissociations, respectively. The calculated results indicate that the fractal dimensions of porous media acquired with this method agree well with the previous study. With the absolute average deviation (AAD) below 10%, the present model provided satisfactory predictions for the dissociation process of methane hydrate and carbon dioxide hydrate.

Fractal dimensions for unsaturated porous media

The analytical expressions of the fractal dimensions for wetting and non-wetting phases for unsaturated porous media are derived and are found to be a function of porosity, maximum and minimum pore sizes as well as saturation. There is no empirical constant in the proposed fractal dimensions. It is also found that the fractal dimensions increase with porosity of a medium and are meaningful only in a certain range of saturation S-w, i.e. S-w > S-min for wetting phase and S-w < S-max for non-wetting phase at a given porosity, based on real porous media for requirements from both fractal theory and experimental observations. The present analysis of the fractal dimensions is verified to be consistent with the existing experimental observations and it makes possible to analyze the transport properties such as permeability, thermal dispersion in unsaturated porous media by fractal theory and technique.

DYNAMIC SCALING AND FRACTAL BEHAVIOR OF SPINODAL PHASE-SEPARATION IN POLYMER MIXTURES OF POLY(METHYL METHACRYLATE) WITH POLY(STYRENE-CO-ACRYLONITRILE)

Dynamic scaling and fractal behaviour of spinodal phase separation is studied in a binary polymer mixture of poly(methyl methacrylate) (PMMA) and poly(styrene-co-acrylonitrile) (SAN). In the later stages of spinodal phase separation, a simple dynamic scaling law was found for the scattering function S(q,t):S(q,t) approximately q(m)-3S approximately (q/q(m)). The possibility of using fractal theory to describe the complex morphology of spinodal phase separation is discussed. In phase separation, morphology exhibits strong self-similarity. The two-dimensional image obtained by optical microscopy can be analysed within the framework of fractal concepts. The results give a fractal dimension of 1.64. This implies that the fractal structure may be the reason for the dynamic scaling behaviour of the structure function.

A comparative study on multiscale fractal dimension descriptors

Fractal theory presents a large number of applications to image and signal analysis. Although the fractal dimension can be used as an image object descriptor, a multiscale approach, such as multiscale fractal dimension (MFD), increases the amount of information extracted from an object. MFD provides a curve which describes object complexity along the scale. However, this curve presents much redundant information, which could be discarded without loss in performance. Thus, it is necessary the use of a descriptor technique to analyze this curve and also to reduce the dimensionality of these data by selecting its meaningful descriptors. This paper shows a comparative study among different techniques for MFD descriptors generation. It compares the use of well-known and state-of-the-art descriptors, such as Fourier, Wavelet, Polynomial Approximation (PA), Functional Data Analysis (FDA), Principal Component Analysis (PCA), Symbolic Aggregate Approximation (SAX), kernel PCA, Independent Component Analysis (ICA), geometrical and statistical features. The descriptors are evaluated in a classification experiment using Linear Discriminant Analysis over the descriptors computed from MFD curves from two data sets: generic shapes and rotated fish contours. Results indicate that PCA, FDA, PA and Wavelet Approximation provide the best MFD descriptors for recognition and classification tasks. (C) 2012 Elsevier B.V. All rights reserved.

Texture analysis by multi-resolution fractal descriptors

This work proposes a novel texture descriptor based on fractal theory. The method is based on the Bouligand- Minkowski descriptors. We decompose the original image recursively into four equal parts. In each recursion step, we estimate the average and the deviation of the Bouligand-Minkowski descriptors computed over each part. Thus, we extract entropy features from both average and deviation. The proposed descriptors are provided by concatenating such measures. The method is tested in a classification experiment under well known datasets, that is, Brodatz and Vistex. The results demonstrate that the novel technique achieves better results than classical and state-of-the-art texture descriptors, such as Local Binary Patterns, Gabor-wavelets and co-occurrence matrix.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.

子午岭林区典型植被下土壤结构及稳定性指标分析

运用分形理论,研究了子午岭林区5种天然次生植被(以6 a天然恢复弃耕地为对照)下土壤结构特征,分析了土壤水稳性团聚体分维、孔隙分维、平均重量直径等3个指标在描述土壤结构稳定性方面的差异。研究表明,相对于弃耕地,各个植被群落均能明显改善土壤结构,降低土壤水稳性团聚体分形维数,提高孔隙分形维数,增强土壤结构的稳定性。土壤水稳性团聚体分维、孔隙分维与>0.25 mm水稳性团聚体含量、土壤有机碳、容重的相关系数均达到了极显著水平,均能作为评价土壤结构稳定性的指标;而团聚体平均重量直径与土壤有机碳含量、容重相关性不显著,只与>5 mm团聚体含量和>0.25 mm团聚体含量有极显著正相关关系,因此,仅可作为大团聚体含量的评价指标。

黄土丘陵区不同土地利用方式土壤团粒结构分形特征

土壤团粒、水稳性团粒和微团粒状况是影响土壤结构和性质的重要因素。运用分形理论研究黄土丘陵区纸坊沟流域7种不同土地利用方式土壤团粒结构的分形特征。结果表明:表层土壤团粒的分形维数在1.641~2.114之间,其大小顺序为人工草地>果园>坡耕地>乔木林地>灌丛>天然草地>人工梯田。土壤水稳性团粒的分形维数在1.774~2.384之间,其变化顺序为果园>乔木林地>人工草地>坡耕地>天然草地>人工梯田>灌丛,二者均表现出≥0.25 mm粒级土壤团粒含量越低,分形维数越高的规律;分形维数与≥5 mm、5~2 mm和≥0.25 mm粒级的土壤团粒、水稳性团粒呈极显著负相关(P<0.01),与<0.25 mm粒级的呈极显著正相关(P<0.01);结合团粒结构体破坏率可知,灌丛土壤结构与稳定性最好,其次为天然草地,人工草地最差;表层土壤微团粒结构分形维数在2.360~2.487之间,大小顺序为天然草地>乔木林地>坡耕地>人工梯田>人工草地>果园>灌丛,灌丛有助于促进土壤微结构形成,分形维数与0.25~0.05 mm和>0.001 mm粒级土壤微团粒含量呈极显著负相关(P<0.01),与<0.001 mm粒级土壤微团粒含量呈极...

农田土壤特性的空间变异性及分形特征

对 2 4 m× 2 4 m面积内的 49个观测点进行了土壤含水量和坚实度的空间变异性研究。根据测定结果 ,用经典统计方法计算了均值、标准差、变异系数等统计特征值 ,认为该法不能定量的描述土壤特性的空间变异性。于是 ,用分形理论计算了含水量与坚实度的分形维数。结果表明 ,在一定范围内 ,土壤具有分形特征 ,分形维数可表示土壤特性在空间的不均一程度

海拉尔盆地贝尔凹陷基岩潜山特征及储层预测

The lithology of the buried hill of Triassic Budate group in Beier depression is epimetamorphic clastic rock and volcanic clastic rock stratum. Recently the favorable hydrocarbon show was discovered in buried hill of base rock, and large-duty industrial oil stream was obtained in some wells in Beier depression. Based on the information of seismos and wells, the tectonic framework, tectonic deformation times and faulted system of the Beier depression are comprehensively studied, then configuration, evolutional history, genetic type and distributed regularity of buried hill are defined. According to observing description and analysis of core sample, well logging and interpretive result of FMI, the lithological component, diagenetic type and diagenetic sequence of buried hill reservoir are confirmed, then reservoir space system of buried hill is distinguished, and vegetal feature, genetic mechanism and distributed regularity of buried hill fissure are researched, at the same time the quantitative relationship is build up between core fissures and fissures interpreted by FMI. After that fundamental supervisory action of fault is defined to the vegetal degree of fissure, and the fissure beneficial places are forecasted using fractal theory and approach. At last the beneficial areas of Budate group reservoir are forecasted by reservoir appraisal parameters optimization such as multivariate gradually regression analysis et. al. and reservoir comprehensive appraisal method such as weighing analyze and clustering procedure et. al. which can provide foundation for the next exploratory disposition. Such production and knowledge are obtained in this text as those: 1. Four structural layers and two faulting systems are developed, and four structural layers are carved up by three bed succession boundary surfaces which creates three tectonic distortional times homology. Three types of buried hill are divided, they are ancient physiognomy buried hill, epigenetic buried hill, and contemporaneous buried hill. 2. Reservoir space of Budate buried hill is mainly secondary pore space and fissure, which distributes near the unconformity and/or inside buried hill in sections. The buried hill reservoir experienced multi-type and multi-stage diagenetic reconstruction, which led to the original porosity disappeared, and multi secondary porosity was created by dissolution, superficial clastation and cataclasis et. al. in diagenetic stage, which including middle crystal pore, inter crystal pore, moldic pore, inter particle emposieu, corrosion pore space and fissure et. al. which improved distinctly reservoir capability of buried hill. 3. The inner reservoir of buried hill in Beier depression is not stratigraphic bedded construction, but is fissure developing place formed by inner fault and broken lithogenetic belt. The fissures in inner reservoir of buried hill are developed unequally with many fissure types, which mainly are high angle fissure and dictyonal fissures and its developing degree and distribution is chiefly controlled by faulting. 4. The results of reservoir comprehensive evaluate and reservoir predicting indicates that advantageous areas of reservoir of buried hill chiefly distributes in Sudeerte, Beixi and Huoduomoer, which comprehensive evaluate mainly Ⅱand Ⅲ type reservoir. The clues and results of this text have directive significance for understanding the hydrocarbon reservoir condition of buried hill in Beier depression, for studying hydrocarbon accumulated mechanism and distributed regularity, and for guiding oil and gas exploration. The results of this text also can enrich and improve nonmarine hydrocarbon accumulated theory.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

溶蚀岩体渗透水力学特性研究

Seepage control in karstic rock masses is one of the most important problems in domestic hydroelectric engineering and mining engineering as well as traffic engineering. At present permeability assessment and leakage analysis of multi-layer karstic rock masses are mainly qualitative, while seldom quantitative. Quantitative analyses of the permeability coefficient and seepage amount are conducted in this report, which will provide a theoretical basis for the study of seepage law and seepage control treatment of karstic rocks. Based on the field measurements in the horizontal grouting galleries of seepage control curtains on the left bank of the Shuibuya Hydropower Project on the Qingjiang river, a hydraulic model is established in this report, and the computation results will provide a scientific basis for optimization of grouting curtain engineering. Following issues are addressed in the report. (1) Based on the in-situ measurements of fissures and karstic cavities in grouting galleries, the characteristics of karstic rock mass is analyzed, and a stochastic structural model of karstic rock masses is set up, which will provide the basis for calculation of the permeability and leakage amount of karstic rock mass. (2) According to the distribution of the measured joints in the grouting galleries and the stochastic results obtained from the stochastic structural model of karstic rock mass between grouting galleries, a formula for computation of permeability tensor of fracturing system is set up, and a computation program is made with Visual Basic language. The computation results will be helpful for zoning of fissured rock masses and calculation of seepage amount as well as optimization of seepage control curtains. (3) Fractal theory is used to describe quantitatively the roughness of conduit walls of karstic systems and the sinuosity of karstic conduits. It is proposed that the roughness coefficient of kastic caves can be expressed by both fractal dimension Ds and Dr that represent respectively the extension sinuosity of karstic caves and the roughness of the conduit walls. The existing formula for calculating the seepage amount of karstic conduits is revised and programmed. The seepage amount of rock masses in the measured grouting galleries is estimated under the condition that no seepage control measures are taken before reservoir impoundment, and the results will be helpful for design and construction optimization of seepage curtains of the Shuibuya hydropower project. This report is one part of the subject "Karstic hydrogeology and the structural model and seepage hydraulics of karstic rock masses", a sub-program of "Study on seepage hydraulics of multi-layer karstic rock masses and its application in seepage control curtain engineering", which is financially supported by the Hubei Provincial key science and technology programme.

花岗岩风化土微观结构及其工程地质特性研究

Saprolite is the residual soil resulted from completely weathered or highly weathered granite and with corestones of parent rock. It is widely distributed in Hong Kong. Slope instability usually happens in this layer of residual soil and thus it is very important to study the engineering geological properties of Saprolite. Due to the relic granitic texture, the deformation and strength characteristics of Saprolite are very different from normal residual soils. In order to investigate the effects of the special microstructure on soil deformation and strength, a series of physical, chemical and mechanical tests were conducted on Saprolite at Kowloon, Hong Kong. The tests include chemical analysis, particle size analysis, mineral composition analysis, mercury injection, consolidation test, direct shear test, triaxial shear test, optical analysis, SEM & TEM analysis, and triaxial shear tests under real-time CT monitoring.Based on the testing results, intensity and degree of weathering were classified, factors affecting and controlling the deformation and strength of Saprolite were identified, and the interaction between those factors were analyzed.The major parameters describing soil microstructure were introduced mainly based on optical thin section analysis results. These parameters are of importance and physical meaning to describe particle shape, particle size distribution (PSD), and for numerical modeling of soil microstructure. A few parameters to depict particle geometry were proposed or improved. These parameters can be used to regenerate the particle shape and its distribution. Fractal dimension of particle shape was proposed to describe irregularity of particle shapes and capacity of space filling quantitatively. And the effect of fractal dimension of particle shape on soil strength was analyzed. At the same time, structural coefficient - a combined parameter which can quantify the overall microstructure of rock or soil was introduced to study Saprolite and the results are very positive. The study emphasized on the fractal characteristics of PSD and pore structure by applying fractal theory and method. With the results from thin section analysis and mercury injection, it was shown that at least two fractal dimensions Dfl(DB) and Df2 (Dw), exist for both PSD and pore structure. The reasons and physical meanings behind multi-fractal dimensions were analyzed. The fractal dimensions were used to calculate the formation depth and weathering rate of granite at Kowloon. As practical applications, correlations and mathematical models for fractal dimensions and engineering properties of soil were established. The correlation between fractal dimensions and mechanical properties of soil shows that the internal friction angle is mainly governed by Dfl 9 corresponding to coarse grain components, while the cohesion depends on Df2 , corresponding to fine grain components. The correlations between the fractal dimension, friction angle and cohesion are positive linear.Fractal models of PSD and pore size distribution were derived theoretically. Fragmentation mechanism of grains was also analyzed from the viewpoint of fractal. A simple function was derived to define the theoretical relationship between the water characteristic curve (WCC) and fractal dimension, based on a number of classical WCC models. This relationship provides a new analytical tool and research method for hydraulic properties in porous media and solute transportation. It also endues fractal dimensions with new physical meanings and facilitates applications of fractal dimensions in water retention characteristics, ground water movement, and environmental engineering.Based on the conclusions from the fractal characteristics of Saprolite, size effect on strength was expressed by fractal dimension. This function is in complete agreement with classical Weibull model and a simple function was derived to represent the relationship between them.In this thesis, the phenomenon of multi-fractal dimensions was theoretically analyzed and verified with WCC and saprolite PSD results, it was then concluded that multi-fractal can describe the characteristics of one object more accurately, compared to single fractal dimension. The multi-fractal of saprolite reflects its structural heterogeneity and changeable stress environment during the evolution history.