A parallel dynamic programming algorithm for multi-reservoir system optimization     

《Advances in water resources》2014

This paper develops a parallel dynamic programming algorithm to optimize the joint operation of a multi-reservoir system. First, a multi-dimensional dynamic programming (DP) model is formulated for a multi-reservoir system. Second, the DP algorithm is parallelized using a peer-to-peer parallel paradigm. The parallelization is based on the distributed memory architecture and the message passing interface (MPI) protocol. We consider both the distributed computing and distributed computer memory in the parallelization. The parallel paradigm aims at reducing the computation time as well as alleviating the computer memory requirement associated with running a multi-dimensional DP model. Next, we test the parallel DP algorithm on the classic, benchmark four-reservoir problem on a high-performance computing (HPC) system with up to 350 cores. Results indicate that the parallel DP algorithm exhibits good performance in parallel efficiency; the parallel DP algorithm is scalable and will not be restricted by the number of cores. Finally, the parallel DP algorithm is applied to a real-world, five-reservoir system in China. The results demonstrate the parallel efficiency and practical utility of the proposed methodology.  相似文献   

基于CUDA的地震相干体并行算法   总被引：1，自引：0，他引：1  

张全 《地质与勘探》2020,56(1):147-153

相干体技术在地震勘探资料解释方面得到了广泛的应用,由于相干体技术处理的对象是三维地震数据体,所以算法运算时间较长。为了缩短解释周期,本文充分发挥GPU并行计算优势,对C3相干体算法进行并行化分析。从硬盘读取数据到GPU上计算相干值并写入硬盘的整个过程进行分析,剔除了冗余数据的读取,完成了C3相干体算法的并行化设计与实现。最后分别对串行算法与并行算法进行性能测试,结果表明本文设计的并行算法在保证精度的前提下达到了16倍左右的加速比,对加快地震资料解释具有重要意义。  相似文献   

基于二次场的可控源电磁法三维有限元-无限元数值模拟   总被引：2，自引：2，他引：0       下载免费PDF全文

张林成  汤井田  任政勇  肖晓 《地球物理学报》2017,60(9):3655-3666

本文提出了一种新的混合有限元-无限元三维可控源电磁法(CSEM)问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.  相似文献   

The varve-related ice-dammed lake events in Glen Roy and vicinity: a new interpretation     

《Proceedings of the Geologists' Association. Geologists' Association》2017,128(1):146-150

This paper presents a new interpretation of the sequence of events in Glen Roy and vicinity during the Loch Lomond Stadial that can be inferred from a detailed varve record constructed by Palmer et al. (2010). 300 years of Younger Dryas glacier advance in the Scottish Highlands are recorded by very thin varves formed in an ice-dammed lake up to 35 km long. At a varve site now occupied by Loch Laggan the lake stood permanently at 260 m, but in Glen Roy varves were also laid down in a lake at 325 m and, later, 350 m caused by glacier advance. Initial ice retreat recorded by a gradual increase in varve thickness was soon followed by much thicker varves. The varve sequences are interrupted by a sand bed caused by sudden drainage of the 350 m lake. The major varves of the Glen Roy sequence show that storminess was still increasing in intensity at least 160 years after glacier retreat had begun. At the Loch Laggan site 15 cm of deformed sediments register an earthquake that produced 3 m faulted uplift of all three Glen Roy shorelines, a response to the abrupt removal of 5 km³ of water when the 260 m lake was catastrophically drained by jökulhlaup. The deformed sediments are immediately followed by varves deposited in a local lake, ice-dammed lake sedimentation now having ceased, having lasted more than 460 years.  相似文献   

并行计算技术及其在勘探地球物理学中的现状与展望   总被引：3，自引：2，他引：1       下载免费PDF全文

黄易  师学明  范建柯  胡文宝 《地球物理学进展》2010,25(2):642-649

本文详细阐述了并行计算技术及其在地球物理勘探数据处理中的发展现状和发展趋势，分析了几个代表性的并行算法实例.这些结果表明，在拥有强大的并行机的基础上，基于并行计算开发环境(MPI和PVM等)设计高效的并行算法，通过分配合理的并行粒度、通信开销、负载平衡等执行高效的并行计算，可以有效加快处理速度、降低成本.目前，并行算法在地震数据处理中应用已较为成熟，近年来向更实用的基于PC机群的并行技术发展.然而，在非地震方法中，并行算法应用较少见文献报道，研究尚处于初级研究阶段.在大地电磁的二维和三维正、反演问题上，并行计算技术逐渐得到越来越多关注和重视.随着资源和能源需求的增长，地球物理勘探向深度和广度快速发展，大幅增长的数据量使得高性能并行计算机和高效的并行算法在勘探地球物理学中的发展和应用将占据愈来愈重要的地位.  相似文献   

Exploring uncertainty in remotely sensed data with parallel coordinate plots     

Yong Ge  Sanping Li  V. Chris Lakhan  Arko Lucieer 《International Journal of Applied Earth Observation and Geoinformation》2009

The existence of uncertainty in classified remotely sensed data necessitates the application of enhanced techniques for identifying and visualizing the various degrees of uncertainty. This paper, therefore, applies the multidimensional graphical data analysis technique of parallel coordinate plots (PCP) to visualize the uncertainty in Landsat Thematic Mapper (TM) data classified by the Maximum Likelihood Classifier (MLC) and Fuzzy C-Means (FCM). The Landsat TM data are from the Yellow River Delta, Shandong Province, China. Image classification with MLC and FCM provides the probability vector and fuzzy membership vector of each pixel. Based on these vectors, the Shannon's entropy (S.E.) of each pixel is calculated. PCPs are then produced for each classification output. The PCP axes denote the posterior probability vector and fuzzy membership vector and two additional axes represent S.E. and the associated degree of uncertainty. The PCPs highlight the distribution of probability values of different land cover types for each pixel, and also reflect the status of pixels with different degrees of uncertainty. Brushing functionality is then added to PCP visualization in order to highlight selected pixels of interest. This not only reduces the visualization uncertainty, but also provides invaluable information on the positional and spectral characteristics of targeted pixels.  相似文献   

A parallel evolutionary strategy based simulation–optimization approach for solving groundwater source identification problems     

Baha Y. Mirghani  Kumar G. Mahinthakumar  Michael E. Tryby  Ranji S. Ranjithan  Emily M. Zechman 《Advances in water resources》2009

Groundwater characterization involves the resolution of unknown system characteristics from observation data, and is often classified as an inverse problem. Inverse problems are difficult to solve due to natural ill-posedness and computational intractability. Here we adopt the use of a simulation–optimization approach that couples a numerical pollutant-transport simulation model with evolutionary search algorithms for solution of the inverse problem. In this approach, the numerical transport model is solved iteratively during the evolutionary search. This process can be computationally intensive since several hundreds to thousands of forward model evaluations are typically required for solution. Given the potential computational intractability of such a simulation–optimization approach, parallel computation is employed to ease and enable the solution of such problems. In this paper, several variations of a groundwater source identification problem is examined in terms of solution quality and computational performance. The computational experiments were performed on the TeraGrid cluster available at the National Center for Supercomputing Applications. The results demonstrate the performance of the parallel simulation–optimization approach in terms of solution quality and computational performance.  相似文献   

A parallel mesh-free contaminant transport model based on the Analytic Element and Streamline Methods     

Karl W. Bandilla  Alan J. RabideauIgor Janković 《Advances in water resources》2009

This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources.  相似文献   

Effect of parallel electric field on electrostatic ion-cyclotron instability in anisotropic plasma in the presence of ion beam and general distribution function—Particle aspect analysis     

Ruchi Mishra 《Planetary and Space Science》2006,54(2):188-199

Dispersion relation, resonant energy transferred, growth rate and marginal instability criteria for the electrostatic ion-cyclotron wave with general loss-cone distribution in low-β anisotropic, homogeneous plasma in the auroral acceleration region are discussed by investigating the trajectories of the charged particles. Effects of the parallel electric field, ion beam velocity, steepness of the loss-cone distribution and temperature anisotropy on resonant energy transferred and growth rate of the instability are discussed. It is found that the effect of the parallel electric field is to stabilize the wave and enhance the transverse acceleration of ions whereas the effect of steepness of loss-cone, ion beam velocity and the temperature anisotropy is to enhance the growth rate and decrease the transverse acceleration of ions. The steepness of the loss-cone also introduces a peak in the growth rate which shifts towards the lower side of the perpendicular wave number with the increasing steepness of the loss-cone.  相似文献   

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition   总被引：3，自引：0，他引：3  

Oliver Baur  Gerrit Austen  Jürgen Kusche 《Journal of Geodesy》2008,82(4-5):207-221

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.  相似文献