基于大地电磁二维反演的MPI和CUDA并行算法研究

大地电磁测深法属于天然场源的电磁勘探方法,是以岩石的电性差异为基础和前提的勘探方法。所采用的大地电磁二维反演方法为共轭梯度法,该方法避免了求解雅可比矩阵,效率较高,但是在将模型剖分较细时,多频率进行计算时效率有待提高。基于大地电磁频率依次独立处理数据的特点,采用了MPI的并行算法来提高效率,为了使得计算效率更高,在MPI的基础上增加了CUDA并行运算方法,用多个进程同时来计算各频率数据,在求解方程的过程中采用CUDA进行加速,得到计算结果。通过对正演和反演图的比较,验证了程序的正确性。对并行算法的效率进行了统计,进程数为2~4时,加速比能达到2.15~3.09,比单一的MPI并行算法的加速比要高,验证了程序的有效性。     

可控源音频大地电磁法在大平山地区矿产勘探中的应用

大平山地区位于长江中下游成矿带宁芜矿集区的北段,矿产资源丰富。但该地区的地形、地质等条件复杂,准确地划定该区的断裂构造走向、埋深及向地下延伸情况对进一步查明区内的铁、铜等矿产资源情况至关重要。以研究区的区域地质概况、地球物理特征为基础,首先对该区的重力场数据的向上延拓处理结果进行分析,科学合理地采用可控源音频大地电磁法(CSAMT)勘探剖面,然后采用CSAMT方法对其进行野外的数据采集。通过对CSAMT数据的视电阻率、相位和二维反演结果进行综合分析。结果表明,大平山地区的断裂分布主要有3条,断裂构造的走向主要呈现为NE、NW。最后,根据钻井资料对CSAMT结果进行验证,发现两者具有较好的一致性。由此可说明CSAMT在地形、地质等复杂条件地区进行断裂构造勘查应用的可行性和有效性,可为大平山及其邻区的下一步矿产勘查和靶区预测等工作指明方向。     

CSAMT法在新疆焉耆盆地铀矿勘查中的应用

根据新疆焉耆盆地铀矿勘查中的勘探深度要求和测区的电阻率及地形特征,采用了可控源音频大地电磁法(CSAMT)进行勘查.对野外采集的数据进行静态校正等处理后,以卡尼亚电阻率和阻抗相位断面图为辅助资料,利用1D和2D反演电阻率断面图进行了地质推断解释.最后根据对断裂构造、构造单元划分等的推断解释资料,制成了新疆焉耆盆地CSAMT测量综合成果图,为在新疆焉耆盆地进一步的铀矿勘查提供了重要依据.     

我国陆上油气田勘探效率评价方法及应用研究

应用统计分析方法对我国石油天然气勘探效率的单年值计算方法进行了改进,分别对探井成功率、每米进尺探明储量、每口探井探明储量和每公里测线探明储量等指标进行了拓展,并应用改进后的勘探效率评价方法对SP石油集团及其所属各油田的总体和老区勘探数据进行了重新计算,并通过勘探效率指标之间的相互联系,对计算结果进行了解释,最终较为科学地获得了有关油气田勘探效率的排序.研究结果表明,在对油气探区的勘探效率进行评价时,应把握油气勘探活动的风险性和较长周期性特征,注意实践中使用的各项指标间的联系,并选择累计平均或移动平均指标作为评价经营业绩的依据.     

CSAMT在新疆库车坳陷铀矿勘探中的应用

为了对新疆库车坳陷地区进行铀矿勘探,采用了勘探深度大、抗干扰能力强的CSAMT法进行铀成矿地质环境的研究.通过野外采集数据,对数据进行静校正处理,钻孔验证了CSAMT测量资料与推断解释地层、岩性的测井资料基本一致,达到了物探电磁测量资料指导钻探的目的.并根据CSAMT测量资料,按照电性层与地层的对应关系,勾划了侏罗系顶板埋深平面图,为进一步勘探工作奠定了基础.     

基于MPI并行及非结构化网格的2.5维CSAMT正演研究

可控源音频大地电磁CSAMT观测数据受发射源的影响,长期以来该方法仅远区信号较弱的数据能被利用,造成近区及过渡区信号较强的资料严重浪费。为避免数据浪费,在Windows平台使用Fortran语言和MPICH2工具,开发了考虑源的CSAMT2.5D正演程序,同时采用了MPI并行技术加速计算。由于实际地下介质边界复杂,程序不同于前人矩形网格,直接采用非结构化三角网格能快速适应复杂边界。通过与一维CSAMT程序对比,验证了2.5D程序的正确性;并通过复杂模型试算,表明考虑源的2.5DCSAMT实用性,为下一步反演工作奠定了理论基础。     

可控源音频大地电磁法(CSAMT)一维全区Levenberg-Marquardt反演方法

将Levenberg-Marquardt这一全新反演方法引入可控源音频大地电磁法(CSAMT)勘探的全区数据反演,提高了CSAMT勘探的深度,避免了测深数据的浪费。Levenberg-Marquardt反演的核心是引入变阻尼的思想,在反演的过程中根据向量各个分量的收敛速度的不同,给定不同的阻尼因子,实现反演的自适应化,这样就大大增加了整个反演的收敛速度和反演精度。反演中采用改进的Bostick变换作为初始模型,减少了反演的迭代次数。     

基于多种位场数据及3D反演结果的综合解释方法

近期位场梯度数据在勘探地球物理观测数据解释中的作用日益增加,随着计算机等硬件设施的进步,处理大型数据能力的增强,重力数据3D反演方法也越来越多的应用到重力数据解释工作中。本文提出结合重力数据、梯度数据及重力数据3D反演结果的综合解释方法,并将此解释方法应用到虎林盆地观测数据、梯度数据及3D反演结果的综合解释工作中。本文首先讨论单个棱柱体的重力数据、梯度数据的计算方法以及重力数据3D反演方法;然后,编制了matlab程序计算重力数据、梯度张量数据并处理重力数据的3D反演;最后,根据计算的数据详细讨论了梯度九分量之间的关系及重力数据反演结果的数值特征,提出结合梯度九分量数据及重力数据3D反演方法联合解释方法。模型试算及解释结果表明重力梯度数据对进一步提高勘探地球物理观测数据的解释精度及确定埋深较浅的异常源位置及水平尺寸等信息起到重要的作用,将其与重力数据3D反演方法联合解释可更准确的刻画异常源的分布特征。最后经过此综合解释方法将虎林盆地划分为8个构造单元。     

基于可控源音频大地电磁法的场源影响校正新方法

基于一维正演模拟,对可控源音频大地电磁法（CSAMT）实测资料处理的场源影响校正问题进行了研究。通过计算远场因子与发射频率、收发距和大地电磁阻抗倒数模三者乘积（综合因数）的关系,实现了简洁有效的查表场源影响校正方法。同时,为提高校正方法解决实际生产问题的能力,又开展了视电阻率数据场源影响校正与反演处理协同联进策略（即反演递归校正法）的研究工作,从而进一步提高了数据处理效果;考虑到反演递归校正处理方法,受初步校正结果影响较大,本文又提出了设计层状背景模型的校正方案,即在工作区前期地质与物性等资料综合研究成果的基础上,通过对工区地质—地电模型的构建,来进一步提高CSAMT实测资料的场源影响校正效果,并经由理论模型算例和实测数据处理验证了方法的可行性和有效性。     

CSAMT在珊瑚钨锡矿勘查中的应用

珊瑚矿区构造复杂,多组断裂交汇发育控制着钨锡矿田（或矿床）的空间分布。为了达到理想的勘探效果,我们采用勘探深度大、分辨率强、观测效率高的CSAMT对测区进行勘探。通过野外采集数据,对数据进行静校正处理,以及和已知的地质资料比较,验证勘探效果,为进一步的勘探工作奠定基础。     

Characterization of Ordovician carbonate reservoirs, southeastern Saskatchewan, Canada

The discovery of the prolific Ordovician Red River reservoirs in 1995 in southeastern Saskatchewan was the catalyst for extensive exploration activity which resulted in the discovery of more than 15 new Red River pools. The best yields of Red River production to date have been from dolomite reservoirs. Understanding the processes of dolomitization is, therefore, crucial for the prediction of the connectivity, spatial distribution and heterogeneity of dolomite reservoirs.The Red River reservoirs in the Midale area consist of 3～4 thin dolomitized zones, with a total thickness of about 20 m, which occur at the top of the Yeoman Formation. Two types of replacement dolomite were recognized in the Red River reservoir: dolomitized burrow infills and dolomitized host matrix. The spatial distribution of dolomite suggests that burrowing organisms played an important role in facilitating the fluid flow in the backfilled sediments. This resulted in penecontemporaneous dolomitization of burrow infills by normal seawater. The dolomite in the host matrix is interpreted as having occurred at shallow burial by evaporitic seawater during precipitation of Lake Almar anhydrite that immediately overlies the Yeoman Formation. However, the low δ18O values of dolomited burrow infills (-5.9‰～ -7.8‰, PDB) and matrix dolomites (-6.6‰～ -8.1‰, avg. -7.4‰ PDB) compared to the estimated values for the late Ordovician marine dolomite could be attributed to modification and alteration of dolomite at higher temperatures during deeper burial, which could also be responsible for its 87Sr/86Sr ratios (0.7084～0.7088) that are higher than suggested for the late Ordovician seawaters (0.7078～0.7080). The trace amounts of saddle dolomite cement in the Red River carbonates are probably related to "cannibalization" of earlier replacement dolomite during the chemical compaction.     

A computer generator for randomly layered structures

AcomputergeneratorforrandomlylayeredstructuresYUJia shun1,2,HEZhen hua2(1.TheInstituteofGeologicalandNuclearSciences,NewZealand;2.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,China)Abstract:Analgorithmisintrod…     

海南岛地体及其毗邻陆缘晚中生代—新生代古地磁研究和构造演化

本文叙述了对海南岛及其毗邻大陆边缘白垩纪到第四纪地层岩石进行古地磁研究的全部工作过程。通过分析岩石中剩余磁矢量的磁偏角及磁倾角的变化,提出海南岛白垩纪以来经历的构造演化模式如下:早期伴随顺时针旋转而向南迁移,后期伴随逆时针转动并向北运移。联系该地区及邻区的地质、地球物理资料,对海南岛上述的构造地体运动提出以下认识:北部湾内早期有一拉张作用,主要是该作用使湾内地壳显著伸长减薄,形成北部湾盆地。从而导致了海南岛的早期构造运动,而海南岛后期的构造运动则主要是受南海海底扩张的影响。海南地体运动规律的阐明对于了解北部湾油气盆地的形成演化有重要的理论和实际意义。     

Exciton Seebeck in Molecular Systems

Various applications relevant to the exciton dynamics,such as the organic solar cell,the large-area organic light-emitting diodes and the thermoelectricity,are operating under temperature gradient.The potential abnormal behavior of the exicton dynamics driven by the temperature difference may affect the efficiency and performance of the corresponding devices.In the above situations,the exciton dynamics under temperature difference is mixed with     

Achieving Unusual States of Matter Under High Pressure

The elongation method,originally proposed by Imamura was further developed for many years in our group.As a method towards O(N)with high efficiency and high accuracy for any dimensional systems.This treatment designed for one-dimensional(ID)polymers is now available for three-dimensional(3D)systems,but geometry optimization is now possible only for 1D-systems.As an approach toward post-Hartree-Fock,it was also extended to     

《Chinese Science Bulletin》SUBJECT INDEX TO VOLUME 59(2014)

正~~     

An on-line scaling method for improving scalability of a database cluster

The explosive growth of the Internet and database applications has driven database to be more scalable and available, and able to support on-line scaling without interrupting service. To support more client's queries without downtime and degrading the response time, more nodes have to be scaled up while the database is running. This paper presents the overview of scalable and available database that satisfies the above characteristics. And we propose a novel on-line scaling method. Our method improves the existing on-line scaling method for fast response time and higher throughputs. Our proposed method reduces unnecessary network use, i.e. , we decrease the number of data copy by reusing the backup data. Also, our on-line scaling operation can be processed parallel by selecting adequate nodes as new node. Our performance study shows that our method results in significant reduction in data copy time.     

An improved R-tree based on childnode''s probability

R-Tree is a good structure for spatial searching. But in this indexing structure,either the sequence of nodes in the same level or sequence of traveling these nodes when queries are made is random. Since the possibility that the object appears in different MBR which have the same parents node is different, if we make the subnode who has the most possibility be traveled first, the time cost will be decreased in most of the cases. In some case, the possibility of a point belong to a rectangle will shows direct proportion with the size of the rectangle. But this conclusion is based on an assumption that the objects are symmetrically distributing in the area and this assumption is not always coming into existence. Now we found a more direct parameter to scale the possibility and made a little change on the structure of R-tree, to increase the possibility of founding the satisfying answer in the front sub trees. We names this structure probability based arranged R-tree (PBAR-tree).     

Integrating GIS Web services based on mediating architecture

The geographic information service is enabled by the advancements in general Web service technology and the focused efforts of the OGC in defining XML-based Web GIS service. Based on these models, this paper addresses the issue of services chaining,the process of combining or pipelining results from several interoperable GIS Web Services to create a customized solution. This paper presents a mediated chaining architecture in which a specific service takes responsibility for performing the process that describes a service chain. We designed the Spatial Information Process Language (SIPL) for dynamic modeling and describing the service chain, also a prototype of the Spatial Information Process Execution Engine (SIPEE) is implemented for executing processes written in SIPL. Discussion of measures to improve the functionality and performance of such system will be included.     

Tracking data management using topology of moving object database system for moving object constrained trajectory in LBS application

Advances in wireless technologies and positioning technologies and spread of wireless devices, an interest in LBS (Location Based Service) is arising. To provide location based service, tracking data should have been stored in moving object database management system (called MODBMS) with proper policies and managed efficiently. So the methods which acquire the location information at regular time intervals then, store and manage have been studied. In this paper, we suggest tracking data management techniques using topology that is corresponding to the moving path of moving object. In our techniques, we update the MODBMS when moving object arrived at a street intersection or a curved road which is represented as the node in topology and predict the location at past and future with attribute of topology and linear function. In this technique, location data that are corresponding to the node in topology are stored, thus reduce the number of update and amount of data. Also in case predicting the location,because topology are used as well as existing location information, accuracy for prediction is increased than applying linear function or spline function.