共查询到20条相似文献，搜索用时 98 毫秒
1.
针对煤矿地质勘探中TEM面积性探测获取的三维数据体,基于Matlab语言对其三维可视化技术进行研究与实践.阐述煤矿TEM三维数据体格式及其可视化要求,提出‘立体图’、‘切片图’和‘等值面图’三种综合性可视化方式,论述使用Matlab语言调用Surfer库函数进行二次开发实现TEM测线数据的批量网格化,解决了煤矿TEM数据体的不规则性问题.同时论述三维可视化技术中关于表面图、三维曲线图、等值线、等值面等三维图形的关键绘制技术与应用函数,实现规则TEM数据体的基本图形表达.在以上工作基础上,作者开发出相应的三维可视化应用软件,通过实例展示其丰富、实用的显示功能.该软件为煤矿TEM资料的处理与解释工作提供了新手段,帮助提高对TEM解释结果的认知能力和资料解释水平. 相似文献
2.
3.
基于CUDA的三维数据并行可视化 总被引：1，自引：0，他引：1
科学计算可视化是发达国家20世纪80年代后期提出并发展起来的一个新的研究领域,它运用计算机图形学和图像处理技术,将三维数据转换为图形及图像在屏幕上显示出来并进行交互处理,主要包括面绘制和体绘制两种方法。光线投射算法是最经典的三维数据体绘制方法,图像质量较高,但计算时间较长,基于CUDA的编程技术利用GPU的多核并行运算功能可显著提高计算速度,实现大规模三维数据的实时可视化。 相似文献
4.
运用Echarts可视化插件和全景技术,将GIS技术与地震信息服务进行有效结合,构建基于WebGIS的地震信息可视化地图显示、标注、查询等功能,实现历史大规模地震目录数据和台站信息在可视化系统中的二、三维一体化展示效果,通过发布震源深度参数的分布图、三维动态折纸模型和全景地震台站,加强公众结合区域特征观看和理解地震活动的直观性和感官效果. 相似文献
5.
6.
在心脏的三维可视化的研究中,传递函数是一个十分重要的工具.然而,心脏复杂的结构使传递函数的设计十分困难.本文介绍了基于可见性直方图的传递函数设计以用来显示心脏双源CT数据.先对心脏双源CT体数据做预处理,提出在视线投射算法中利用基于可见性直方图的心脏传递函数对心脏数据进行三维显示,并且提出半自动方法来确定传递函数.通过这些方法,用户可以通过调节传递函数突出显示感兴趣的部分,甚至可以观察到心腔的许多细节.最后,利用CUDA技术对计算及显示进行加速,对三维体数据进行实时动态显示.实验表明,利用本文提出的方法可以有效地去掉心脏双源CT中的肋骨,肺静脉等数据,突出显示心肌组织及血液,并且通过传递函数可以看到主动脉瓣等心脏内部结构,对于心外科医生临床诊断有一定的意义. 相似文献
7.
8.
运用Echarts可视化插件和全景技术,将GIS技术与地震信息服务进行有效结合,构建基于WebGIS的地震信息可视化地图显示、标注、查询等功能,实现历史大规模地震目录数据和台站信息在可视化系统中的二、三维一体化展示效果,通过发布震源深度参数的分布图、三维动态折纸模型和全景地震台站,加强公众结合区域特征观看和理解地震活动的直观性和感官效果。 相似文献
9.
针对地壳构造形成的动力学机制,在广泛收集区域地质、地质构造和地球物理等资料的基础上,着重对重力数据进行了数字化、坐标、重力公式、投影方式和比例尺的统一化处理,进而进行了网格化处理.为将区域异常和局部异常分离开来,以便以对深部构造的系统研究,笔者选取趋势分析法分别对研究区内预处理后的布格重力异常数据(Δg)进行了三维多项式迭代拟合计算,得到了区域布格重力异常二阶、五阶、十阶趋势分析结果.继而,利用二级近似公式迭代法对布格重力异常五阶趋势分析区域异常数据进行了三维运算处理,得到了黄海及周边地区的莫霍面埋深值.分析了黄海及邻区莫霍面起伏特征,并进行了深部构造区划,探讨了深部断裂构造与莫霍面起伏间的成因联系,为深部构造和活动断裂演化的地球动力学研究提供了依据. 相似文献
10.
三维可视化解释技术是一项全新的地震解释技术,与传统的地震资料解释有着本质的不同,它是从三维可视化显示出发,以地质体或三维研究区块为单元,采用点、线、面和体相结合的数据体空间可视化解释。该技术结合相干体技术可进行复杂断层的解释与组合;应用空间域层位自动追踪技术能够准确落实构造形态;采用目标雕刻技术,利用反演数据体可实现立体透视岩性体的空间展布及厚度趋势。将该技术应用在大港油田官西地区复杂断块构造解释和歧南西斜坡隐蔽油气藏储层预测两项实际工作中,取得了良好效果,表明可视化解释技术在石油勘探中具有广阔的应用前景,是提高勘探效率和勘探成功率的有效技术手段。 相似文献
11.
ICT切片重构模型的剖切及应用 总被引：3，自引：0，他引：3
在ICT切片重构模型的可视化中,为了分析重构模型的任意截面形状,并且观察该截面的内部结构,需要对重构模型进行剖切.本文针对剖切提出一种简单方法,该方法利用剖切面对重构模型所在的立方体图像空间进行剖切,生成交点集合,然后将空间剖切平面映射到二维图像平面上,最后计算空间剖切平面图像的灰度值,得到剖面图像.通过剖面图像,可以直观地看到内部结构形状,为后续无损检测奠定基础. 相似文献
12.
We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions. 相似文献
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions. 相似文献
13.
《Geofísica Internacional》2014,53(3):289-308
Edge enhancement is an element of analysis to derive the spatial structure of satellite images. Two methods to extract edges from multispectral satellite images are presented. A multispectral image is modeled as a vector field with a number of dimensions equal to the number of bands in the image. In this model, a pixel is defined as a vector formed by a number of elements equal to the number of bands. Two vector operators are applied to such vector field. In our first method, we extend the definition of the gradient. In this extension, the vector difference of the window central pixel with neighboring pixels is obtained. A multispectral image is then generated where each pixel represents the maximum change in spectral response in the image in any direction. This image is named a multispectral gradient. The other method, considers the generalization of the Laplacian by means of an η-dimensional Fourier transform. This image is named a multispectral Laplacian. The vector operators perform a simultaneous extraction of edge-content in the spectral bands of a multispectral image. Our methods are parameter-free. Our methods work for a multispectral image of any number of bands. Two examples are discussed that involve multispectral satellite images at two scales. We compare our results with widely used edge enhancement procedures. The evaluation of results shows better performance of proposed methods when compared to widely used edge operators. 相似文献
14.
Earth degassing specifically of carbon dioxide CO2 is of increasing interest with respect to the global carbon budget, related climate effects, earthquake and volcano eruption mechanisms, as well as plant physiological reactions in gas-rich environments. Investigations in all of these disciplines require the detection of surface CO2 degassing structures and quantification of their emissions. We introduce minimal thermal change detection based on infrared imaging as a new remote sensing tool for the detection of earth surface thermal anomalies suiting among others to discover earth degassing locations of any origin. The method allows for seamless areal search and monitoring of degassing structures in any terrain. As proof of concept infrared imaging measurements were performed at the Bossoleto vent on the eastern master fault of the Siena Graben (Tuscany, Italy). It is known for the migration of a large amount of CO2-rich gas from deep geothermal reservoirs. Field data acquired confirmed the qualification of the method. Detection of CO2 degassing locations from infrared image time series worked reliably and optimal detection conditions were identified (dry, calm, cloudless weather between dusk and dawn). A simple model of heat exchange processes involved and observed was developed. In a first attempt this model was applied to determine the gas exit temperature, the area of gas thermal reach and the gas flux from recorded image series. It is the first method that allows remote areal survey of mofette fields and the associated CO2 flux quantification sole from infrared image time series. 相似文献
15.
Xiaoyu Zhang Qunshan Zhang George A. McMechan Gladys Gonzalez 《Geophysical Prospecting》2017,65(3):669-686
Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region. 相似文献
16.
We propose a fast method for imaging potential field sources. The new method is a variant of the “Depth from Extreme Points,” which yields an image of a quantity proportional to the source distribution (magnetization or density). Such transformed field is here transformed into source‐density units by determining a constant with adequate physical dimension by a linear regression of the observed field versus the field computed from the “Depth from Extreme Points” image. Such source images are often smooth and too extended, reflecting the loss of spatial resolution for increasing altitudes. Consequently, they also present too low values of the source density. We here show that this initial image can be improved and made more compact to achieve a more realistic model, which reproduces a field consistent with the observed one. The new algorithm, which is called “Compact Depth from Extreme Points” iteratively produces different source distributions models, with an increasing degree of compactness and, correspondingly, increasing source‐density values. This is done through weighting the model with a compacting function. The compacting function may be conveniently expressed as a matrix that is modified at any iteration, based on the model obtained in the previous step. At any iteration step the process may be stopped when the density reaches values higher than prefixed bounds based on known or assumed geological information. As no matrix inversion is needed, the method is fast and allows analysing massive datasets. Due to the high stability of the “Depth from Extreme Points” transformation, the algorithm may be also applied to any derivatives of the measured field, thus yielding an improved resolution. The method is investigated by application to 2D and 3D synthetic gravity source distributions, and the imaged sources are a good reconstruction of the geometry and density distributions of the causative bodies. Finally, the method is applied to microgravity data to model underground crypts in St. Venceslas Church, Tovacov, Czech Republic. 相似文献
17.
With a high-resolution seismic survey, we targeted the shallow structural features of the Cheb Basin, Northwest Bohemia. The aim of this study was to evaluate these features’ relation to seismic activity below the Cheb Basin with the first reflection seismic image of parts of the basin. We recorded a densely sampled, 3200 m long seismic survey stretching from Kopanina towards SE, into the basin. From the recordings of 170 shots at 192 geophones with 384-m total spread, we obtained a reflection seismic image down to the basement at 300-m depth, supplemented by a shallow tomographic velocity model for the uppermost 40 m. Strong lateral velocity contrasts indicate Quaternary river terraces, into which the Plesná river incises. We observe reverse faults in the lower basin sediments, which we interpret as signs of a push-up structure related to the N-S oriented Po?átky-Plesná Zone (PPZ). However, we do not observe any vertical fault offsets in the younger sedimentary layers, which suggests that any normal or reverse faulting must be older than 20 Ma. The image agrees well with previous lithostratigraphic borehole data and previous sedimentological and tectonic models. Our explanation for the formation of the escarpment at the eastern bank of the Plesná valley, previously interpreted as outcrop of the PPZ, relies solely on incision of the Plesná river into weak sediments. 相似文献
18.
Improving the gradient of the image‐domain objective function using quantitative migration for a more robust migration velocity analysis 
下载免费PDF全文
下载免费PDF全文 Mark Noble 《Geophysical Prospecting》2015,63(2):391-404
Migration velocity analysis aims at determining the background velocity model. Classical artefacts, such as migration smiles, are observed on subsurface offset common image gathers, due to spatial and frequency data limitations. We analyse their impact on the differential semblance functional and on its gradient with respect to the model. In particular, the differential semblance functional is not necessarily minimum at the expected value. Tapers are classically applied on common image gathers to partly reduce these artefacts. Here, we first observe that the migrated image can be defined as the first gradient of an objective function formulated in the data‐domain. For an automatic and more robust formulation, we introduce a weight in the original data‐domain objective function. The weight is determined such that the Hessian resembles a Dirac function. In that way, we extend quantitative migration to the subsurface‐offset domain. This is an automatic way to compensate for illumination. We analyse the modified scheme on a very simple 2D case and on a more complex velocity model to show how migration velocity analysis becomes more robust. 相似文献
19.
基于Markov链模型的储层岩相随机模拟 总被引：3，自引：4，他引：3
在油气储层随机建模研究中，基于Markov链模型的方法是一类较受欢迎的技术，同时也是一类不成熟的技术，问题的症结之一在于侧向的转移概率矩阵很难求取，针对这种情况，作者在深入理解Walther相律的基础上，借鉴模拟退火算法的相应思路，提出了一种岩相模拟的新方法，该方法依据不同岩相的百分比进行随机模拟得到一幅初始图象，而后以按岩相组织剖面得到的垂向和侧向的岩相转移概率矩阵的相似性作为判别标准对图象进行扰动，直至得到满意的图象，二维模型试算结果表明了这种岩相随机模拟方法的可行性。 相似文献
20.
Paul Sava 《Geophysical Prospecting》2015,63(5):1086-1096
Waveform inversion is a velocity‐model‐building technique based on full waveforms as the input and seismic wavefields as the information carrier. Conventional waveform inversion is implemented in the data domain. However, similar techniques referred to as image‐domain wavefield tomography can be formulated in the image domain and use a seismic image as the input and seismic wavefields as the information carrier. The objective function for the image‐domain approach is designed to optimize the coherency of reflections in extended common‐image gathers. The function applies a penalty operator to the gathers, thus highlighting image inaccuracies arising from the velocity model error. Minimizing the objective function optimizes the model and improves the image quality. The gradient of the objective function is computed using the adjoint state method in a way similar to that in the analogous data‐domain implementation. We propose an image‐domain velocity‐model building method using extended common‐image‐point space‐ and time‐lag gathers constructed sparsely at reflections in the image. The gathers are effective in reconstructing the velocity model in complex geologic environments and can be used as an economical replacement for conventional common‐image gathers in wave‐equation tomography. A test on the Marmousi model illustrates successful updating of the velocity model using common‐image‐point gathers and resulting improved image quality. 相似文献