CN106368689B - 一种快速生成测井成像图的方法 - Google Patents
一种快速生成测井成像图的方法 Download PDFInfo
- Publication number
- CN106368689B CN106368689B CN201510434941.0A CN201510434941A CN106368689B CN 106368689 B CN106368689 B CN 106368689B CN 201510434941 A CN201510434941 A CN 201510434941A CN 106368689 B CN106368689 B CN 106368689B
- Authority
- CN
- China
- Prior art keywords
- logging
- data
- imaging graph
- sampling
- line number
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
Abstract
本发明公开了一种快速生成测井成像图的方法，涉及成像测井技术领域。该方法，在保证生成的测井成像图符合使用要求的前提下，通过固定绘图区高度方向的像素行数，保证生成测井图像时，即便采样数据有数G大小，每次加载的磁盘数据也不会超过绘图区高度(单位是像素)*每行数据数目，而一般情况下绘图区的高度仅仅是几百像素，所以可以保证每次放大、缩小或滚动浏览操作，只提取几百行数据，提高测井成像图的生成速度，提高系统读取数据的性能。
Description
技术领域
本发明涉及成像测井技术领域，尤其涉及一种快速生成测井成像图的方法。
背景技术
成像测井(imaging logging)是根据钻孔中地球物理场的观测，对井壁和井周围物体进行物理参数成像的技术，是地理学的一门基础学科。
当今油田勘探中，测井成像数据非常大，一般会达到数千兆字节规模。用户在浏览测井成像图时，经常需要随时对成像图进行放大、缩小或滚动浏览，而每进行一次放大、缩小或滚动操作，就需要读取一次测井成像数据并重新生成测井成像图，由于测井成像数据非常大，所以，每次数据读取并生成图像，都需要较长的时间，因此，用户每次放大、缩小或滚动操作时，不仅测井成像图生成速度慢，而且会严重占用设备的资源，影响设备的性能。
发明内容
本发明的目的在于提供一种快速生成测井成像图的方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明采用的技术方案如下：
一种快速生成测井成像图的方法，包括如下步骤：
S1，获取测井成像图的采样行数；
S2，确定绘图区高度的像素行数；
S3，利用如下公式，计算每读取一行测井数据需要间隔的采样行数：
C＝A/B，
式中，
A，测井成像图的采样行数；
B，绘图区高度的像素行数；
C，每读取一行测井数据需要间隔的采样行数；
S4，根据每读取一行测井数据需要间隔的采样行数，间隔读取测井数据，根据测井数据生成测井成像图。
优选地，S1包括如下步骤：
A1，获取待测井筒的范围；
A2，获取测井的采样间隔；
A3，根据井筒范围和采样间隔，计算测井成像图的采样行数。
优选地，A3采用如下公式进行计算：
测井成像图的采样行数＝待测井筒的范围/测井的采样间隔。
优选地，所述采样间隔为0.001-0.01米。
更优选地，所述采样间隔为0.001米。
本发明的有益效果是：本发明实施例提供的快速生成测井成像图的方法，在保证生成的测井成像图符合使用要求的前提下，通过固定绘图区高度方向的像素行数，保证生成测井图像时，即便采样数据有数G大小，每次加载的磁盘数据也不会超过绘图区高度(单位是像素)*每行数据数目，而一般情况下绘图区的高度仅仅是几百像素，所以可以保证每次放大、缩小或滚动浏览操作，只提取几百行数据，提高测井成像图的生成速度，提高系统读取数据的性能。
附图说明
图1是本发明实施例提供的快速生成测井成像图的方法流程示意图；
图2a是原始测井成像图；
图2b是采用本实施例提供的方法得到的一个测井成像图；
图2c是采用本实施例提供的方法得到的另一个测井成像图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
如图1所示，本发明实施例提供了一种快速生成测井成像图的方法，包括如下步骤：
S1，获取测井成像图的采样行数；
S2，确定绘图区高度的像素行数；
S3，利用如下公式，计算每读取一行测井数据需要间隔的采样行数：
C＝A/B，
式中，
A，测井成像图的采样行数；
B，绘图区高度的像素行数；
C，每读取一行测井数据需要间隔的采样行数；
S4，根据每读取一行测井数据需要间隔的采样行数，间隔读取测井数据，根据测井数据生成测井成像图。
上述方法可以在初始的测井成像图生成时使用，即，在实际测井生成成像图时，确定采样行数。比如，要生成地下2000米到3000米范围的成像图，测井数据是每隔0.001米采样一行数据，那么2000米到3000米就会有(3000-2000)/0.001＝1000000行数据，绘图区只有500行像素，因此读取数据时就每隔1000000/500＝2000行读一行数据，共读取500行数据，将该500行数据生成测井成像图。从而极大的降低了图像的数据大小，则用户在后续的放大、缩小或滚动浏览中，每次操作只需要读取500行数据即可，与原始得到的采样数据相比，极大的降低了读取数据量，提高了测井成像图的生成速度，以及系统读取数据的性能。
上述方法还可以在已经采用现有技术生成了测井成像图，而在用户进行浏览操作时使用，比如，采用现有技术生成的原有的测井成像图的高度是1000像素(或者说原图的数据有1000行采样数据)，如果按照现有技术进行图像解析并重新生成图像，则每次都需要读取1000行数据，而如果使用本实施例提供的方法，则由于绘图区高度是100像素，则每次只需要读取100像素的数据，其他的900像素的数据被丢弃掉，具体地，按照本实施例的方法，每隔10个像素采样一行图像，其余9行图像丢弃掉，从而不仅保证了生成的新的图像符合测井的使用要求，而且极大的降低了读取数据量，提高了测井成像图的生成速度，以及系统读取数据的性能。
上述所述的成像方案，不是压缩图像，而是丢弃某些像素，因此，采用上述方法生成的图像势必会所有失真。例如一个1000*1000像素的图像，在一个100*100像素的设备上显示，势必会产生失真。但是，本发明实施例提供的方法，在丢弃测井采样数据时，采用的方法为：按照相同间隔选择读取的数据行，将其余的数据丢弃。采用上述方法得到的图像，虽然有所失真，但是，可以保证生成的测井成像图符合测井的使用要求。如图2a-图2c所示，图2a为没有丢弃数据时的原始图像，图2b为丢弃二分之一的数据时生成的图像，图2c为丢弃三分之二的数据时生成的图像，将图2a-2c进行比较，可以看出，虽然丢弃数据后，图像有所失真，但是依然在可接受范围之内,依然可以用于测井分析。
丢失数据是因为在进行概要浏览，能起到模糊效果实现大概的相貌。这就像从高空俯视大地一样，距离越高，细节越是看不到。但总体的地貌却能大体确定。
因此，按照本发明实施例提供的方法得到的图像，会提供放大缩小功能，放大或者缩小都会重新生成一次图像，放大到1:1的比例就是真实图像。就像你从高空看大地一样，距离越近，看的越清楚越真实，所以，虽然有所失真，但是，放大后还是很清晰的。
由于一般成像图的数据组织，每行的数据是连续存储在一块的，因此，上述方法在实际使用过程中，如果需要间隔读取N行数据，磁盘指针需要移动跳跃N次，完成数据的读取。
本发明的一个优选实施例中，S1包括如下步骤：
A1，获取待测井筒的范围；
A2，获取测井的采样间隔；
A3，根据井筒范围和采样间隔，计算测井成像图的采样行数。
其中，A3可以采用如下公式进行计算：
测井成像图的采样行数＝待测井筒的范围/测井的采样间隔。
比如，要生成地下2000米到3000米范围的成像图，测井数据是每隔0.001米采样一行数据，那么2000米到3000米就会有(3000-2000)/0.001＝1000000行数据。
本发明的一个实施例中，所述采样间隔可以为0.001-0.01米。
在本发明的一个优选实施例中，所述采样间隔为0.001米。
通过采用本发明公开的上述技术方案，得到了如下有益的效果：本发明实施例提供的快速生成测井成像图的方法，在保证生成的测井成像图符合使用要求的前提下，通过固定绘图区高度方向的像素行数，保证生成测井图像时，即便采样数据有数G大小，每次加载的磁盘数据也不会超过绘图区高度(单位是像素)*每行数据数目，而一般情况下绘图区的高度仅仅是几百像素，所以可以保证每次放大、缩小或滚动浏览操作，只提取几百行数据，提高测井成像图的生成速度，提高系统读取数据的性能。
本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
本领域人员应该理解的是，上述实施例提供的方法步骤的时序可根据实际情况进行适应性调整，也可根据实际情况并发进行。
上述实施例涉及的方法中的全部或部分步骤可以通过程序来指令相关的硬件来完成，所述的程序可以存储于计算机设备可读取的存储介质中，用于执行上述各实施例方法所述的全部或部分步骤。所述计算机设备，例如：个人计算机、服务器、网络设备、智能移动终端、智能家居设备、穿戴式智能设备、车载智能设备等；所述的存储介质，例如：RAM、ROM、磁碟、磁带、光盘、闪存、U盘、移动硬盘、存储卡、记忆棒、网络服务器存储、网络云存储等。
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。
Claims (1)
1.一种快速生成测井成像图的方法，其特征在于，包括如下步骤：
S1，获取测井成像图的采样行数；
S2，确定绘图区高度的像素行数；
S3，利用如下公式，计算每读取一行测井数据需要间隔的采样行数：
C=A/B，
式中，
A，测井成像图的采样行数；
B，绘图区高度的像素行数；
C，每读取一行测井数据需要间隔的采样行数；
S4，根据每读取一行测井数据需要间隔的采样行数，间隔读取测井数据，根据测井数据生成测井成像图；按照相同间隔选择数据行，将其余的数据丢弃；
S1包括如下步骤：
A1，获取待测井筒的范围；
A2，获取测井的采样间隔；
A3，根据井筒范围和采样间隔，计算测井成像图的采样行数；
测井成像图的采样行数=待测井筒的范围/测井的采样间隔；
所述采样间隔为0.001-0.01米。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510434941.0A CN106368689B (zh) | 2015-07-22 | 2015-07-22 | 一种快速生成测井成像图的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510434941.0A CN106368689B (zh) | 2015-07-22 | 2015-07-22 | 一种快速生成测井成像图的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106368689A CN106368689A (zh) | 2017-02-01 |
| CN106368689B true CN106368689B (zh) | 2020-02-18 |
Family
ID=57879893
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510434941.0A Active CN106368689B (zh) | 2015-07-22 | 2015-07-22 | 一种快速生成测井成像图的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106368689B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843935B (zh) * | 2017-10-24 | 2019-12-13 | 北京吉奥特能源科技有限责任公司 | Xrmi电成像图像自动生成方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5502686A (en) * | 1994-08-01 | 1996-03-26 | Western Atlas International | Method and apparatus for imaging a borehole sidewall |
| CN1361475A (zh) * | 2000-12-29 | 2002-07-31 | 中国科学院水生生物研究所 | 一种轮纹图象的处理方法 |
| US20090290814A1 (en) * | 2008-05-21 | 2009-11-26 | Yoo Sin Sung | Image compensation apparatus and method |
| CN102044058B (zh) * | 2010-11-22 | 2012-11-21 | Tcl集团股份有限公司 | 一种在电子设备中处理图片的方法、系统及电子设备 |
-
2015
- 2015-07-22 CN CN201510434941.0A patent/CN106368689B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106368689A (zh) | 2017-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schaff et al. | Improvements in magnitude precision, using the statistics of relative amplitudes measured by cross correlation | |
| Sun et al. | Mapping floods due to Hurricane Sandy using NPP VIIRS and ATMS data and geotagged Flickr imagery | |
| KR101015365B1 (ko) | 시각적으로 대표되는 비디오 썸네일을 생성하는 컴퓨터 구현 방법, 컴퓨팅 장치, 컴퓨터, 비디오 장치, 비디오 기록 장치 및 컴퓨터 판독가능 기록매체 | |
| US11042991B2 (en) | Determining multiple camera positions from multiple videos | |
| Zhu et al. | A study of National Lightning Detection Network responses to natural lightning based on ground truth data acquired at LOG with emphasis on cloud discharge activity | |
| CN103871058B (zh) | 基于压缩采样矩阵分解的红外小目标检测方法 | |
| Hermance | Stabilizing high‐order, non‐classical harmonic analysis of NDVI data for average annual models by damping model roughness | |
| Wang et al. | Lightning attachment processes of three natural lightning discharges | |
| CN107179550B (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| Byrne et al. | Estimating badger social-group abundance in the Republic of Ireland using cross-validated species distribution modelling | |
| Scaioni et al. | Analysis of spatial sensor network observations during landslide simulation experiments | |
| Hillier et al. | Production and preservation of the smallest drumlins | |
| CN106368689B (zh) | 一种快速生成测井成像图的方法 | |
| Dong et al. | The S-STK/LTK algorithm for arrival time picking of microseismic signals | |
| Gao et al. | A modified extreme event-based synchronicity measure for climate time series | |
| Mariethoz et al. | Feature-preserving interpolation and filtering of environmental time series | |
| Shirer et al. | Estimating the correlation dimension of atmospheric time series | |
| CN114200522B (zh) | 深度域地震子波提取方法和装置、存储介质及电子设备 | |
| CN109522520B (zh) | 地下水水位波动与多元因素的多重小波相干性分析方法 | |
| CN109636727B (zh) | 一种超分辨率重建影像空间分辨率评价方法 | |
| CN110991267A (zh) | 基于图像或视频人群计数的密度图生成方法和装置 | |
| Peng et al. | Early magnitude estimation for the MW7. 9 Wenchuan earthquake using progressively expanded P-wave time window | |
| Blázquez et al. | On the sensitivity of the fractal dimension of the skeleton of sea surface temperature signatures from upwelling events on spatial resolution and threshold: application on imagery of the Caribbean Sea | |
| CN106204501B (zh) | 一种压缩感知恢复方法 | |
| CN111309093B (zh) | 一种理想波形模型的建立方法、装置、设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 834000 Jiyun Road, Karamay District, Karamay City, Xinjiang Uygur Autonomous Region, 197-2 Patentee after: HONGYOU SOFTWARE Co.,Ltd. Address before: 834000 No. 22, long march road, the Xinjiang Uygur Autonomous Region, Karamay Patentee before: KARAMAY HONGYOU SOFTWARE Co.,Ltd. |