An improved shuffled frog leaping algorithm and its application in the optimization of cascade reservoir operation

ABSTRACT

A new improved shuffled frog leaping algorithm (SFLA), the chaos catfish effect SFLA (CCESFLA), is proposed by coupling a local “refine search” mechanism and a “global incentive adjustment” mechanism. Chaotic technology is introduced in the local “refine search” mechanism to improve local search ability, by implementing more refined local search around the optimal individuals. The catfish effect mechanism is adopted in the “global incentive adjustment” mechanism to improve global convergence, by motivating the frogs to “jump out” of the local steady state. The operation optimization by the CCESFLA is carried out taking the Li Xianjiang cascade reservoirs in China as an example. Compared with SFLA, particle swarm optimization, immune SFLA and cloud SFLA, the average annual power generation using the CCESFLA can be increased by 6.7, 7.5, 3.0 and 0.8%, respectively. The convergence process of the CCESFLA is more stable, and its execution time is the least of the three improved SFLAs.     

倾角分解共反射面元叠加方法

共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距剖面成像方法,传统的CRS叠加实现是以数据驱动的方式对属性参数进行自动搜索并对其进行优化合成相应的CRS叠加算子,通过该算子进行叠加能够得到信噪比和连续性更高的零偏移距剖面.但是数据驱动的实现方式带来了不可避免的“倾角歧视现象”,它造成了弱有效反射信号损失和运动学特征失真的问题.本文提出的倾角分解CRS叠加方法成功解决了上述问题,使CRS叠加方法更具实用价值.     

The variation of cosmogenic Kr and Xe in a core from Estherville mesosiderite: direct evidence that the lunar131Xe anomaly is a depth effect

Kr and Xe were measured by a stepwise heating technique in three samples of a drill core in the “Minnesota” fragment of the Estherville mesosiderite. The cosmogenic⁷⁸Kr/⁸³Kr decreased from the “top” sample to the “bottom” sample(“top” = 0.163 ± 0.005, “bottom” = 0.151 ± 0.005) while the cosmogenic¹³¹Xe/¹²⁶Xe ratio increased(“top” = 5.58 ± 0.35, “bottom” = 6.92 ± 0.17). Cosmic-ray track studies have shown that the “top” sample was indeed closer to the preatmospheric surface than the “bottom” sample by ～ 10 cm. This is the first direct evidence, in a sample of known geometry, that the cosmogenic¹³¹Xe/¹²⁶Xe ratio increases as a function of depth, and as such, confirms the hypothesis that the lunar¹³¹Xe anomaly is a bona fide depth effect due to resonance neutron capture in¹³⁰Ba.     

强烈地震震源破裂和深层过程与地震短临预测探索

强烈地震震源区和其周边地域的介质,在力源作用下开始破裂,即微破裂,而当其在应力不断积累、逐渐发展并形成具有一定规模的“破裂链”时,即应力集中达到临界状态时,震源介质所积累的巨大能量以波动的形式瞬间释放,并冲出地表发生地震,且在地表产生一系列的破坏和大型断裂.在强烈地震发生后依据设定的初始模型和远、近地震台网记录,通过震源机制解计算,可求得震源介质的破裂长度和破裂过程.事实表明,一系列大地震的“孕育”、发生和发展,由震源深处到地表均呈现出强烈的破裂效应,其破裂长度可由几公里、几十公里、乃至千余公里长.有的地震在发生前确亦存在着破裂响应和迹象.为此,若在地表和深井中能进行破裂效应的观测,以捕捉初始微破裂和其形成“破裂链”的动力过程,这对于短临地震预测可能是一个十分重要的途径.     

Estimation of source function and medium response function by autoregressive method

By modelling seismograms as “low” and “high” order autoregressive (AR) processes, the source function and the medium response function are separated from a single channel seismogram. Akaike's final prediction error is used as a statistic to select the appropriate “low” and “high” AR order of the process. Case studies of synthetic data show that the recovered source and reflectivity functions compare very well with the input functions. Using this method, arrivals of the surface reflected P phases of five explosions from the Soviet region and of two earthquakes from Kamchatka, recorded at Gauribidanur Seismic Array (GBA), India, are identified. Certain features of the source and source region of these events are also inferred.     

莺歌海盆地“暗点”型气藏地震识别技术研究

莺歌海盆地乐东区深层与东方区具有类似的成藏条件,是高温高压领域获取下一个重大勘探突破的最现实区域.但是该区已钻遇的“亮点”型储层物性差,普遍表现为低孔-特低渗特征,而区域基础规律研究表明,区域潜在甜点储层在地震上应该表现为“暗点”型特征.但是,“暗点”型储层在叠后地震资料上无法有效识别,如何准确预测“暗点”型储层是乐东区勘探新突破的重点攻关方向.本文针对“暗点”识别难的问题,提出了基于叠前地震资料的储层预测技术.从岩石物理出发系统性地梳理该区“暗点”型储层的地球物理特征,总结出其具有“叠加无反射、近道弱波峰、远道弱波谷、道集极性反转”的规律,根据该规律建立了Fstack属性重构算法,利用该算法开展甜点储层筛选并进行“暗点”储层预测,结合道集资料、反演结果等地球物理信息,最终确定潜在“暗点”型气层的分布范围.钻井结果证实该技术对乐东区深层潜在“暗点”型储层具有很好的预测效果.     

Groundwater potential zones using a combination of geospatial technology and geophysical approach: case study in Dehradun,India

ABSTRACT

A combination of geospatial, geophysical and statistical models using satellite data, the weighted index overlay (WIO) method and two-dimensional electrical resistivity tomography (2D-ERT) is applied to generate the highest potential groundwater area and to further explore the groundwater in Dehradun, India. The results show that of 19.7 km² total basin area, 0.26% falls under the “poor” category as a prospect zone for groundwater, 4.3% is “moderate”, 10.10% “moderately good”, 4.9% “good” and 0.17% “very good”. In addition, the demonstration of the geophysical survey is presented, in which Purkal Youth Society Division (PYSD) site is categorized as a shallow aquifer zone and the Guru Nanak Fifth Centenary School (GNFCS) site is a deeper aquifer zone. Our study emphasizes remote sensing and geographic information system integrated with a geophysical survey to support prospecting the most probable area and confirm the existence of groundwater.     

Magnetic polarization of the Schumann resonances: An asymptotic theory

Horizontal components of extremely low-frequency (ELF) fields in the gyrotropic Earth–ionosphere cavity are expressed through the fields of an isotropic cavity. A linear approximation of the perturbation method with respect to two small parameters is used. These are the impedance z₀ of the isotropic ionosphere and the ratio γ≡z₂/z₁ of the off-diagonal to diagonal elements of the surface impedance matrix of the gyrotropic ionosphere. It is shown that the effect of splitting of the Schumann resonance frequencies caused by the existence of the d.c. geomagnetic field results in an “apparent” angular shift of all ground-based sources of ELF radiation toward the geographic West (“refraction of ELF fields in the gyrotropic Earth–ionosphere cavity”) and excitation of an elliptical magnetic polarization. An analytical relation for the mentioned phenomena is derived. Explicit expressions for polarization characteristics of the electromagnetic ELF noise excited by an arbitrary number of equatorial thunderstorm centers are obtained. Expected values of the “apparent” angular shift for these sources are estimated. Analysis of the estimated magnitudes shows that locations of thunderstorm centers and isolated Q-bursts determined with the use of multi-component ELF measurements should be corrected to exclude the “gyrotropic refraction” effect. Equations to calculate the necessary corrections are presented. The results are generalized to the case of an arbitrary latitudinal location of the sources.     

滑坡危险性等级区划方法研究

以四川省地理环境条件为研究背景,根据相关单位提供的资料,对四川省滑坡危险性等级进行了区划研究。将地形地貌、地层岩性、断裂构造、河流水系、降雨量、地震烈度等6项作为滑坡的主要影响因素,采用模式识别方法进行滑坡危险性等级区划,对算法中所涉及到的危险系数计算公式、因子权重分析、阈值选取等进行了一系列控制试验,验证了算法的可行性以及降低算法中存在的不确定性。采用“识别率”和“改变率”2个准则来判断“分类结果作为新的训练集”,即RTS试验的收敛性,从而给出识别率高、改变率稳定的分类结果,以及能合理反映识别结果的最佳参数。通过3次逐级识别分类,将四川省滑坡危险度划分为7个等级,区划结果与实际滑坡发生情况吻合。本文方法同样适用于其它地区的滑坡危险性等级区划。     

地震预测模型优化方法研究

在分析现有基于观测指标“异常-正常”的“二态”地震前兆模式不足的基础上, 提出了异常可表现为多种状态的“多态”前兆模式, 给出“多态”前兆模式下预测效能判定指标, 进一步提出了基于预测效能最优的单项预测模型参数的选择方法。 在广义(时间、 空间、 时-空联合的预测指标)“多态”前兆模式下, 将所建立的模型应用于我国华北和南北带地区, 分别确定了两地区活断层预测地震(M_S≥6)空间分布的最优模型(分别为断层周围20 km、 30 km)及其预测效能(分别为0.42、 0.34)。     

Numerical investigation of alternative fracture stiffness measures and their respective scaling behaviours

We study the mechanical deformation of fractures under normal stress, via tangent and specific fracture stiffnesses, for different length scales using numerical simulations and analytical insights. First, we revisit an equivalent elastic layer model that leads to two expressions: the tangent stiffness is the sum of an “intrinsic” stiffness and the normal stress, and the specific stiffness is the tangent stiffness divided by the fracture aperture at current stress. Second, we simulate the deformation of rough fractures using a boundary element method where fracture surfaces represented by elastic asperities on an elastic half‐space follow a self‐affine distribution. A large number of statistically identical “parent” fractures are generated, from which sub‐fractures of smaller dimensions are extracted. The self‐affine distribution implies that the stress‐free fracture aperture increases with fracture length with a power law in agreement with the chosen Hurst exponent. All simulated fractures exhibit an increase in the specific stiffness with stress and an average decrease with increase in length consistent with field observations. The simulated specific and tangent stiffnesses are well described by the equivalent layer model provided the “intrinsic” stiffness slightly decreases with fracture length following a power law. By combining numerical simulations and the analytical model, the effect of scale and stress on fracture stiffness measures can be easily separated using the concept of “intrinsic” stiffness. We learn that the primary reason for the variability in specific stiffness with length comes from the fact that the typical aperture of the self‐affine fractures itself scales with the length of the fractures.     

Biological activity and physicochemical properties of thermally activated water types: Mollusk-based bioassay

Physicochemical and bioassay methods are used to carry out a comprehensive study of the properties of two types of activated waters: “thermally activated,” i.e., thermal water, referred to as “white spring” or “cold boiling water” in folk medicine, and “melt” water after quick freeze. The test object for bioassay were apple shell mollusks Pomacea canaliculata from a laboratory culture used to carry out patented and widely tested methods of freshwater quality bioassay. Oxygen and hydrogen peroxide concentrations and pH were determined in water. The integral biological activity of water was evaluated by a method based on the use of absorption spectra to evaluate the rate of redox transformation in the quinone-hydroquinol system. Similar studies were carried out with three types of reference samples: settled original tap water, the same water cooled to the freezing point, and water artificially enriched with hydrogen peroxide. The water activated by thermal impact has been found to have higher biological activity than the original tap water.     

Verfahren zur katalytisch-oxydativen Entsulfidierung von Abwässern

The oxidative conversion of strongly alkaline sulphide concentrations in wastewaters from the wet molybdenite enrichment is successfully performed within a relatively short time by the application of atmospheric oxygen, “Winkler adsorbent” (activated form of ash from Winkler generators) or “Winkler generator ash” and chlorinated lime at the same time. The reaction times were 60… 90 min at residual sulphide concentrations of 1… 10 mg/l. The initial concentrations of sulphide were between 400 and 1,200 mg/l. It was proved by comparative investigations that the “Winkler absorbent” or the “Winkler generator ash” have a catalytical effect during the oxidation by air. Therefore, also, e.g., at an initial concentration of sulphide of 600 mg/l only 4 g/l chlorinated lime were required, which corresponded to only about 0.13 times the stoichiometrically calculated quantity with the available chlorine content being taken into account. Flotation cells are suitable as reactors.     

ESSAI DE FILTRAGE NON-LINEAIRE APPLIQUE AUX PROFILS AEROMAGNETIQUES*

Digital aeromagnetic profiles recorded with a high accuracy can be filtered with an electronic computer. Linear filtering is not very efficient, because the frequency spectrum of a given anomaly is very wide. But other methods are possible which, through a step by step analysis of the profile, leave some categories of anomalies unaltered while they remove the others completely. Our method uses as a criterium the width of the anomalies; it distinguishes and processes differently “bell-shaped” and “multi-legged” anomalies. Some examples of the use of the method are shown.     

A NONITERATIVE PROCEDURE FOR INVERTING PLANE-WAVE REFLECTION DATA AT SEVERAL ANGLES OF INCIDENCE USING THE RICCATI EQUATION*

Various exact methods of inverting the complete waveform of vertical seismic reflection data to produce acoustic impedance profiles have been suggested. These inverse methods generally remain valid for nonvertical, plane-wave data, provided total reflection does not occur. Thus, in principle, the “seismogram” at each ray parameter in a slant stack can be interpreted separately. Rather than invert each plane-wave seismogram separately, they can all be interpreted simultaneously and an “average” model thus obtained. Inversion for both the velocity and the density also becomes possible when two or more plane-wave seismograms are simultaneously inverted. The theory for a noniterative inversion method, based on the time-domain Riccati equation, is discussed. Numerical examples of inversions using this technique on synthetic data demonstrate its numerical stability and the advantage of simultaneous inversion of several seismograms to reduce the effect of noise in the data and increase the stability of the inversion process.     

A canonical correlation approach to the determination of homogeneous regions for regional flood estimation of ungauged basins

Abstract

A canonical correlation method for determining the homogeneous regions used for estimating flood characteristics of ungauged basins is described. The method emphasizes graphical and quantitative analysis of relationships between the basin and flood variables before the data of the gauged basins are used for estimating the flood variables of the ungauged basin. The method can be used for both homogeneous regions, determined a priori by clustering algorithms in the space of the flood-related canonical variables, as well as for “regions of influence” or “neighbourhoods” centred on the point representing the estimated location of the ungauged basin in that space.     

Assessment of workplace accident risks in underground collieries by integrating a multi-goal cause-and-effect analysis method with MCDM sensitivity analysis

This study introduces a method using a multi-goal fuzzy cognitive map (FCM) and multi-criteria decision making based on sensitivity analysis to assess the risks associated with working accidents in underground collieries. Safety, stoppage in operation, and operational and capital costs are considered as the main goals during the FCM process with significant emphasis on safety. Workplace accidents data from Kerman underground collieries are statistically evaluated to find the degrees of occurrence probability, severity, and work-disability duration as the main risk factors. The causes and effects of accidents are analyzed using FCM based on three goals and the effects of risk factors. A sensitivity analysis on the weights of the goals is conducted with the aim of increasing the workplace safety in TOPSIS environment after solving the designed multi-goal FCM. Results indicate that “gas poisoning,” “roof fall,” and “debris and destruction” take the first three ranks and impose high risks to the system. By contrast, “collision, hit, and crash” presents the lowest risk among all accidents.     

A practice‐oriented estimation of the failure probability of building structures

In the presented practice‐oriented probabilistic approach for the seismic performance assessment of building structures, the SAC‐FEMA method, which is a part of the broader PEER probabilistic framework and permits probability assessment in closed form, is combined with the pushover‐based N2 method. The most demanding part of the PEER probabilistic framework, that is incremental dynamic analysis, is replaced by the much simpler N2 method, which requires considerably less input data and much less computational time, but which can, nevertheless, often provide: acceptable estimates for the mean values of the structural response. Using some additional simplifying assumptions that are consistent with seismic code procedures, an explicit equation for a quick estimation of the annual probability of “failure” (i.e. the probability of exceeding the near collapse limit state) of a structure can be derived, which is appropriate for practical applications, provided that predetermined default values for the dispersion measures are available. In the paper, this simplified approach is summarized and applied to the estimation of the “failure” probability of reinforced concrete frame buildings representing both old structures, not designed for earthquake resistance, and new structures designed according to Eurocode 8. The results of the analyses indicate a high probability of the “failure” of buildings, which have not been designed for seismic loads. For a building designed according to a modern code, the conservatively determined probability of “failure” is about 30 times less but still significant (about 1% over the lifetime of the structure). Copyright © 2011 John Wiley & Sons, Ltd.     

Optimal Estimation of Two-Dimensional Contaminant Transport

The state-space estimation technique presented herein provides a method for obtaining optimal estimates of concentrations for two-dimensional plumes in ground water. The concentration of a plume was defined as the state variable. The technique uses the Kalman filter and involves combining two independent estimates of plume concentrations. One estimate is called the process modeling data and the other is called the “measurement” data. The process modeling data is obtained from a numerical model. The “measurement” data is obtained from field measurements; however, for illustration in this paper it was generated by a different numerical model than the one used to obtain the process modeling data. The state-space technique produces a distribution of contaminant concentrations that is more accurate than either of the distributions generated by the process modeling or the “measurement” data. An example is presented to show that the technique produces significant improvements in the prediction of plume concentration distributions.     

A framework for sourcing of evaporation between saturated and unsaturated zone in bare soil condition

ABSTRACT

Sourcing subsurface evaporation (E_ss) into groundwater (E_g) and unsaturated zone (E_u) components has received little scientific attention so far, despite its importance in water management and agriculture. We propose a novel sourcing framework, with its implementation in dedicated post-processing software called SOURCE (used along with the HYDRUS1D model), to study evaporation sourcing dynamics, define quantitatively “shallow” and “deep” water table conditions and test the applicability of water table fluctuation (WTF) and “bucket” methods for estimation of E_g and E_u separately.

For the “shallow” and “deep” water table we propose E_g?>?0.95E_ss and E_g = 0 criteria, respectively. Assessment of the WTF method allowed sourcing of very small fluxes otherwise neglected by standard hydrological methods. Sourcing with SOURCE software was more accurate than the standard “bucket” method mainly because of greater flexibility in spatio-temporal discretization. This study emphasized the dry condition relevance of groundwater evaporation which should be analysed by applying coupled flow of heat, vapour and liquid water.

Editor D. Koutsoyiannis; Associate editor S. Kanae