基于双重逆极限空间的地貌信息数据挖掘方法

为了挖掘不完备的地貌数据信息,将地貌信息描述为拓扑空间的形式,从逆极限定义出发,利用映射诱导法建立与原地貌信息空间拓扑同构的双重逆极限空间作为目标信息空间,然后利用拓扑传递性分析原地貌信息空间与目标信息空间的拓扑相似度,实现不完备的地貌信息空间数据挖掘.本文利用GIS软件,结合地貌的坡度、坡向等地貌特征因子分别模拟了地貌信息数据挖掘前后的模型并进行了对比,结果表明此方法是可行的.     

反演模型分辨率的估算方法

反演问题的时空间分辨率或称时空分辨长度是评估模型精细程度的重要参数,决定了该模型应用的范围和价值,但是分辨长度估算却是比反演更复杂和麻烦的数学问题。除了层析成像中广泛利用理论模型恢复试验定性提取空间分辨长度外,通过求解分辨率矩阵可定量获得分辨长度。通过矩阵操作给出的分辨率矩阵包括三类：直接分辨率矩阵、正则化分辨率矩阵和混合分辨率矩阵。这三类矩阵包含了反演本身不同侧面的信息,因此在一个反演应用中,同时提供这三类分辨率矩阵可更全面地评估反演模型分辨率分布。最近An（2012）提出了从大量随机理论模型及其解中统计出分辨率矩阵的方法。这种分辨率矩阵是从模拟真实反演实验的输入和输出模型中通过反演得到的,因此这种分辨率矩阵更能反映整个反演所涉及到的更多因素和过程;同时由于这种分辨率矩阵计算过程无需进行矩阵操作且不依赖于具体正演和反演方法,因此可以被应用于更普遍的反演问题。实际应用证明统计分辨率分析方法适用于对二维和三维层析成像反演模型进行分辨率分析。     

曲波域不规则地震数据提高分辨率技术研究

受野外观测条件的限制,采集的地震数据体通常不规则,并缺失一部分数据道。传统的单道提高分辨率方法无法兼顾横向地震信息,处理结果存在空间一致性问题。为此,本文提出在曲波域内进行不规则地震数据,通过曲波变换实现对地震数据的稀疏表征,将提高分辨率问题转化为曲波域1-范数约束的稀疏促进求解,得到规则化的高分辨率地震数据体。该方法避免传统单道提高分辨率方法存在的局限性,在提高分辨率的同时,能够恢复缺失的地震数据、压制随机噪声,进而提高地震数据的完备性,模型和实际资料试算,验证了该方法的正确性、有效性和适用性。      

双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

工业CT系统空间分辨率两种测试方法分析与评价

空间分辨率是CT图像质量的两项重要性能指标,其测试工作也是CT设备性能测试的重要内容。本文提出工业CT系统空间分辨率常用的两种测试方法,编写基于圆盘法的空间分辨率测试软件,并结合测试软件进行空间分辨率的测试。两种测试方法的对比实验在三套工业CT系统中开展,并在相同试验条件下实施、数据采集、计算和分析。结果发现,在三套工业CT设备测试结果中,基于圆盘法的测试结果要明显低于基于标准中规定的线对卡测试方法测试结果,本文根据傅立叶变换讨论两种测试方法测试结果产生差异的原因。      

基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演研究-理论模型

使用广义最小残量方法迭代求解三维频率-空间域声波方程,反演时使用多尺度、多重网格的策略,探讨了如何快速实现高分辨率的三维频率-空间域迭代法声波全波形速度反演.通过对理论模型进行三维频率-空间域迭代法声波全波形反演数值试验,证实该方法的计算速度快、计算效率高,反演所得速度的分辨率高.从而为基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演成像打下方法基础.     

基于模型空间压缩的大地电磁三维反演研究

本文提出了一种基于模型空间压缩技术的大地电磁三维反演方法.该方法在传统大地电磁三维反演理论的基础上,通过小波变换将待反演的空间域模型参数映射到小波域进行反演,获得小波域更新模型后再通过小波逆变换得到空间域反演模型.由于小波变换具有压缩特性和多尺度分辨能力,本文反演方法可在一定程度上提高反演分辨率.为了提高反演效率,我们针对基于L_1范数的模型约束求解不易收敛的反演问题,提出了一种基于模型粗糙度的简单有效的预条件处理技术.为验证本文算法的有效性,本文首先对经典的"棋盘"模型进行三维反演测试.反演结果表明本文算法的反演效率与传统方法相当,但对于深部异常体具有更好的分辨能力.最后,我们通过对实测数据反演进一步验证本文算法的有效性.     

地震活动信息熵及时空分维在地震预报中的意义

本文选取新疆境内强震的高发区——乌恰地区作为研究区域,分别采用合成缺信量、合成信息熵、信息熵、容量维和信息维等五种不同的方法,通过时间域和空间域分析对乌恰地区1980—1988年发生的三次强震进行了动态追踪的研究。结果表明,无论是从熵的观点(均匀程度),还是从分维数的观点(有序程度)出发,强震前的减熵(降维)过程是明显的,研究表明,这种强震前的减熵(降维)过程可能成为强震震源体处于临界破裂的一种特征物理量,因而可作为强震预测以及重点危险区动态监测的一种方法。     

小波分析在遥感图像融合中的应用

利用小波分析的方法对多分辨率、多波段的遥感图像信息进行融合，原始图像经小波变换，分解成多个子图像，选择适宜的小波基，进行分块处理，再做逆变换，进行图像融合重构，使得融合后的图像最大限度地保留多波段光谱信息的同时，提高了图像的空间分辨率及局部目标的清晰度，与传统方法相比，小波分析方法具有明显的优越性。     

A geographic data model for representing ground water systems

The Arc Hydro ground water data model is a geographic data model for representing spatial and temporal ground water information within a geographic information system (GIS). The data model is a standardized representation of ground water systems within a spatial database that provides a public domain template for GIS users to store, document, and analyze commonly used spatial and temporal ground water data sets. This paper describes the data model framework, a simplified version of the complete ground water data model that includes two-dimensional and three-dimensional (3D) object classes for representing aquifers, wells, and borehole data, and the 3D geospatial context in which these data exist. The framework data model also includes tabular objects for representing temporal information such as water levels and water quality samples that are related with spatial features.     

X射线三维显微镜及其典型应用

X射线三维显微镜（显微CT）是目前发展最热门、成长最快的CT成像技术门类之一,在微纳制造技术、新材料以及电子科学等领域起着十分重要的作用并得到广泛的应用。X射线三维显微镜能够以亚微米的细节分辨能力对被检测对象内部结构进行无损的三维成像,是一种新型的检测设备。本文介绍了三英精密仪器有限公司研制的nano Voxel-2000系列X射线三维显微镜的结构和工作原理,展示了该仪器在系统设计、成像方法及应用分析方面的特色。该系列的X射线三维显微镜最大管电压为150 k V,采用5μm的微焦源实现低于500 nm的成像分辨率,能够进行不同尺寸样品、不同分辨率的成像。在此基础上,探讨了X射线显微镜的典型应用,主要对样品的内部结构、形貌进行表征。希望典型应用的实例能起到抛砖引玉的作用,吸引更多学者利用X射线三维显微成像技术开展科学研究,拓宽X射线三维显微镜的应用领域。     

中国大陆西南地区二次强震活动的信息和信息熵的分析

本文从现代统计物理学中熵的概念出发,用所定义的合成缺信量和地震活动信息熵的方法,讨论了中国西南地区永善、龙陵二次强震前后地震活动的信息和信息熵,采用时间域、空间域以及空间扫描三种不同的分析方法,结果表明,强震前缺信量的降低和信息熵的减熵过程是明显的,并且具有熵值下降过程中发震,随后熵值恢复的孕震模式。     

3D Mapping of calcite and a demonstration of its relevance to permeability evolution in reactive fractures

There is a need to better understand reaction-induced changes in fluid transport in fractured shales, caprocks and reservoirs, especially in the context of emerging energy technologies, including geologic carbon sequestration, unconventional natural gas, and enhanced geothermal systems. We developed a method for 3D calcite mapping in rock specimens. Such information is critical in reactive transport modeling, which relies on information about the locations and accessible surface area of reactive minerals. We focused on calcite because it is a mineral whose dissolution could lead to substantial pathway alteration because of its high solubility, fast reactivity, and abundance in sedimentary rocks. Our approach combines X-ray computed tomography (XCT) and scanning electron microscopy. The method was developed and demonstrated for a fractured limestone core containing about 50% calcite, which was 2.5 cm in diameter and 3.5 cm in length and had been scanned using XCT. The core was subsequently sectioned and energy dispersive X-ray spectroscopy was used to determine elemental signatures for mineral identification and mapping. Back-scattered electron microscopy was used to identify features for co-location. Finally, image analysis resulted in characteristic grayscale intensities of X-ray attenuation that identify calcite. This attenuation mapping ultimately produced a binary segmented 3D image of the spatial distribution of calcite in the entire core. To demonstrate the value of this information, permeability changes were investigated for hypothetical fractures created by eroding calcite from 2D rock surfaces. Fluid flow was simulated using a 2D steady state model. The resulting increases in permeability were profoundly influenced by the degree to which calcite is contiguous along the flow path. If there are bands of less reactive minerals perpendicular to the direction of flow, fracture permeability may be an order of magnitude smaller than when calcite is contiguous. These results emphasize the importance of characterizing spatial distribution of calcite in heterogeneous rocks that also contain a similar abundance of less reactive minerals.     

Time-space multiscale analysis by use of tensor product wavelets and its application to hydrology and GRACE data

This paper presents a wavelet analysis of temporal and spatial variations of the Earth’s gravitational potential based on tensor product wavelets. The time-space wavelet concept is realized by combining Legendre wavelets for the time domain and spherical wavelets for the space domain. In consequence, a multiresolution analysis for both temporal and spatial resolution is formulated within a unified concept. The method is then numerically realized by using first synthetically generated data and finally two real data sets.     

Model ALADIN as regional climate model for Central and Eastern Europe

Results obtained with two versions of the Limited Area Model (LAM) ALADIN over differently sized integration domains (large, intermediate and small) in the European area are presented in order to investigate both the general model performance and the influence of domain choice on the quality of obtained results. The aim is also to illustrate the issues related to the strategy of selection of the optimal integration domain. Each of these studies has been performed with two versions of the ALADIN model: the first one is ALADIN-CLIMATE developed at CNRM/Météo-France, the second one is ALADIN-CLIMATE/CZ prepared at the Czech Hydrometeorological Institute (CHMI). This leaves us with total of six experiments forced by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 reanalysis data. The west Balkan domain covering Bulgaria is used as an evaluation region for investigation of the temporal and spatial properties of simulated precipitation and temperature fields. This region has been selected for its challenging orography making the results obtained here a valuable source for studies leading to further developments in climate modeling. It was found that size of the domain strongly affects the quality of obtained results. We have found that the largest domain reproduces the spatial characteristics of climate (such as bias) very well, but its use results in a poor representation of temporal aspects, which are however captured very well in experiments over both smaller domains. Our findings suggest that there is no optimal choice of domain size, securing the best results for both spatial and temporal evaluation.     

基于迭代法的同轴法X射线相衬成像的相位复原

相位复原问题是指在只有强度可以测量的情况下,从强度数据中提取出物体的相位信息。相位复原是同轴法X射线相衬成像的关键问题之一。本文采用误差递减法,通过限制物体的空间大小,结合测量到的强度数据,在物体的空域和频域都进行约束。通过反复迭代,进行带宽外推,最终实现物面和像面的相位复原。仿真实验证明了此迭代算法的收敛性,并且以塑料管为样品,在上海光源采集到同轴法X射线相衬成像的实验数据,对上述算法进行了验证。     

OpendTect系统在层序地层研究中的应用

在等时层序格架建立的基础上,利用地震资料研究各体系域内部沉积体的空间展布、叠置样式和沉积演化过程,进而预测各体系域中砂体的分布,是层序地层学应用于勘探实践的重要内容.OpendTect系统在三维地震倾角扫描的基础上,①以地震采样点空间信息为导向,实现了层序界面控制下的地震小层自动追踪,可以精细刻画出体系域内部沉积体地震反射结构和空间展布;②通过对各层序单元内部地震小层追踪结果进行层拉平,实现了时间域到Wheeler域的自动变换,可以更加直观地分析地层沉积演化历史;③通过沿沉积体内部小层面提取地震属性,实现了等时地层切片的提取,减少了地震属性多解性,提高了预测精度.OpendTect系统提供了一种新的基于层序地层的地震综合解释技术和思路.     

X射线断层扫描投影图像背景不一致的校准

在X射线断层扫描(X-ray CT)成像过程中,经常会出现射线源、样品和探测器三者之间的相对移动,导致采集的投影图像(Projection image)以及亮场图像(flat field image)的空间位置配准关系出现偏差。这直接影响了投影图像背景的正常扣除,以及后续的CT图像重建。本文利用互相关方法,求出各投影图像背景与探测器原始采集位置之间的不一致量。然后将投影图像和亮场图像作相应变换,最后扣除背景得到校准后的图像。由于该方法对图像灰度较为敏感,所以在对实际样品校准之前,本文先对该样品投影图像及亮场图像进行了照度归一化。实验结果表明,该方法能够有效校正由射线源、样品和探测器之间相对移动造成的CT投影图像背景不一致,获得较好的CT投影图像和重建切片。      

3D wavefield extrapolation with optimal separable approximation

An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.