第20卷 第1期地 球 物 理 学 进 展Vol .20 N o .12005年3月(页码:71～77)
PROGRE SS I N GEOPHYS ICS
M a rc h 2005
大数据量地球物理测井绘图关键技术研究
李 河, 王祝文*, 李舟波, 王朝辉
(吉林大学地球探测科学与技术学院,长春130026)
摘 要 在开发地球物理测井解释软件过程中,绘图所涉及的数据量很大,在绘图过程中要进行大量的比例运算、坐标转换及坐标映射等,特别是绘制巨幅、需要连续滚屏显示的声电成像测井图像,数据量和运算量更大,因此,就会经常出现内存不足、图形刷新速度过慢、图形重叠和闪烁等现象.在开发测井解释系统G eo W ell 过程中,采用内存映射文件、内存DC 绘图、局部重绘等技术实现了大数据量巨幅测井图形的快速度、高质量绘制以及图形的连续滚动显示,提高了绘图效率,为微机环境下开发大数据量复杂测井绘图软件提供了技术参考.关键词 测井软件,内存映射文件,内存DC 绘图,局部重绘,连续滚动显示
中图分类号 TP311;P 631 文献标识码 A 文章编号 1004-2903(2005)01-0071-07
Th e k ey te chn o lo gy re search o f geoph y s ica l
w e ll lo gg ing p lo t w ith m as s da ta
L I He , WA NG Zhu -w en , L I Zhou -bo , WANG Zhao-hu i
(C olle ge of G eo-e x plor a ti on Science &Tec hnology ,Jili n Un i vers it y ,C hangchun 130026,C hi na )
Ab s tra c t In the process of deve lop ing the soft wa re of geophys ica l we ll logg ing interpreta tion ,a great volu m e of da ta is involved in p lotting .And in the course of p lotting ,m uch sca le co m pu ting ,coord ina te transfor m a tion and coord ina te m app ing is to be under way.Pa rticularly when p lotting huge acoustic i m ag ing logg ing and electric i m ag ing logg ing i m -age wh ich de m ands continuous scrolling d is p lay ,there are m uch m ore vol um es of da ta and co m pu ting work.There-f ore ,it frequently turns out to be so m e pheno m ena ,s uch as m em ory insufficiency ,s peed decreas ing of p lot refresh ,p lot overlap and screen flicker etc .W h ile deve lop ing G eo W ell of the we ll logg ing interpreta tion syste m,by us ing tec hnology of m em ory m app ing file ,m e m ory DC p lott ing and pa rt redraw ing ,the huge well logg ing i m age w ith m ass da ta can be drawn fleetly in h igh qua lity ,and the p lot can be scrolled continuously.It a lso i m proves p lott ing efficiency and offers techn ica l reference f or deve lop ing the soft wa re of co m p lica ted w ell logg ing p lot w ith m ass da ta in the PC env iron m ent .
K e yw ord s w ell logg ing soft wa re ,m e m ory m app ing file ,m e m ory DC p lotting ,part redraw ing ,con tinuous scrolling d is-p lay
收稿日期 2004-07-10; 修回日期 2004-08-20.
基金项目 国家自然科学基金"九五"重大项目和大庆石油管理局项目(49894190-42)联合资助.
作者简介 李河(1973-),男,博士研究生,主要从事地球物理测井软件开发、科学计算可视化及图像识别等方面的研究.
(E -m a il :L i hestrong @sina.co m )*
通讯作者:王祝文(1961-),男,教授,博士生导师,主要从事地球物理测井方法研究.
0 引 言
计算机硬件的发展速度很快,其运行速度更是越来越快,但是在很多复杂的图形系统软件中,如G IS 系统,CAD 图形系统,地球物理测井绘图系统等,由于要绘制大量的图形元素,图形显示一次需要大量的时间
[1—8]
.在测井绘图系统中,特别是声电成
像图像的绘制,所涉及的数据量更大,绘制的图形幅度也很大,而且所绘测井图经常是经过大量的数值计算、比例运算和坐标转换及坐标映射之后才画到屏幕上或打印输出的.在开发测井绘图软件时,当要绘制图形的井段很长的时候,如在深度比例为1:200的情况下绘制3000m 的图,考虑多道,只绘常规曲线数据量就非常庞大,如果绘制成像图像,数据
地球物理学进展20卷
量就更大.如此大的数据量,如果想一次把所有数据读入内存然后绘图,会由于占用内存太大而不可能;同时,如果采用常规绘图和滚屏方法,一次把整幅图全部绘出,必然会出现图形刷新速度过慢、图形闪烁和图形重叠等现象,因此必须采取措施,提高图形的绘制和显示速度[9—15].针对上述问题,采用内存映射文件、内存DC、局部重绘等技术,获得了良好的绘图效果.
1利用内存映射文件实现大数据量测井文件的快速读取
在测井数据处理和绘图过程中,会遇到读取大数据量测井数据文件的问题.在对数据进行处理或绘图时,最常用的方法是定义数组来临时保存数据,但是当数据量很大时,要想一次读取所有数据,数组也必须定义得很大,这样就会发生内存不足的问题,如果用传统的通过I/O文件读取技术进行数据读取的方式,会因为数据量大,速度慢不能满足要求[16,17].本文采用内存映射文件技术来解决这些问题[18].
1.1内存映射文件的工作原理
内存映射文件允许应用程序与通过指针访问动态内存相同的方式访问磁盘上的文件[19].因此,可以在进程地址空间中将磁盘上的文件部分或者全部映射到特定地址范围,然后通过指针访问内存映射文件中的内容.一旦文件被映射,就可以访问它,就像整个文件已经加载到内存一样,从而可以不必对文件执行I/O操作,这对大数据量文件来说存取效率较高.内存映射文件分三种情况:
(1)可执行文件的内存映射,主要由操作系统自身使用;
(2)数据文件的内存映射;
(3)借助于系统页面交换文件的共享内存映射.应用程序可以使用后面两种内存映射文件.系统在进行工作时,首先把数据文件的一部分映射到虚拟地址空间内,但不提交RA M,存取这段内存指令时会产生一个页面异常,系统捕获到这个异常后,分配一页RA M,并把它映射到当前进程发生异常的地址处,然后系统把文件中相应的数据读到这个页面中,继续执行刚才产生异常的指令.这就是应用程序自己不需要调用文件I/O函数的原因,也是内存映射文件技术的工作机理.
1.2使用内存映射文件的一般步骤
用内存映射文件进行文件输入输出操作需要下
面几个步骤[20]:
(1)调用C rea teF ile函数,以适当的方式创建或打开一个文件核心对象.
(2)把C rea teF ile函数返回的文件句柄作为参数传给C rea teF ile M app ing函数,C rea teF ile M ap-p i ng函数创建一个文件映射核心对象的适当属性.
(3)创建了文件映射核心对象后,调用M ap-V ie wO f F ile函数,将整个文件的部分或全部映射到内存.
(4)利用M apV ie wO f F ile函数返回的指针来读写文件数据.
(5)用Un M ap V ie wO f F ile函数撤销对文件映射核心对象的映射.
(6)使用C lose Hand le函数关闭文件映射核心对象.
(7)用C lose Hand le函数关闭文件核心对象.
1.3基于内存映射文件的XT F格式测井数据文件
的快速读取
在开发测井解释处理系统时,了解测井数据文件的记录格式非常重要,只有知道所处理的测井数据的格式,才能正确地进行数据的读取、格式转换、处理等工作.在开发G eo W ell测井解释系统时,系统要处理ECLIPS5700数控测井系统XTF格式文件,并进行绘图.当XTF文件中含有各类成像测井数据时,文件可能很大,必须采取有效的方法来读取数据,才能够解决内存不足、数据读取和处理速度慢的问题.
XTF文件包含的信息量非常大.它由标题块和数据块两大部分组成,最小的组成单元是记录,每个记录的长度均为4096字节,标题块通常包括8个记录,数据块包含的记录个数由曲线的深度范围来定, XTF文件允许不同特性的曲线并存,例如曲线的起始深度、结束深度、采样率等都可以不同[21].XTF 文件的标题块后紧跟数据块,数据块用于存放XTF 文件的数据.每条曲线的数据存放都包括两部分,即曲线头记录(4096字节)和数据记录(其长度根据每条曲线的起始深度、结束深度、采样率以及曲线类型等确定).曲线头记录的是每条曲线的信息,包括曲线名、起始和结束深度、采样率、维数等,有些内容与标题块中的相同;数据记录中存放曲线的数据,如果曲线的字节数不是4096的整数倍,不够4096的根据需要用0或其它空值补齐.XTF文件整体结构如图1所示.
27
1期李 河,等:
大数据量地球物理测井绘图关键技术研究
图1 XTF 格式文件结构
F ig 1 The structure of XTF for m a t file
在实际应用中,数据处理和绘图都存在任意时刻读取任意井段的任意井信息和任意曲线信息及曲线数据的需要.因此,基于内存映射文件,采用面向对象方法,设计了类,能够快速方便地读取当前打开井的任意深度段的曲线数据,由于篇幅所限,类的详细设计作者另文介绍.一些文献中只是介绍内存映射文件的原理和使用中的几个步骤,具体如何使用的介绍并不多见,这里举一个具体例子:假设文件已经映射到内存,需要读取XTF 文件标题信息的8个标题记录,具体如下:
BOOL CP lot D a taF ile O bject ::ReadXTF8Head(){
int Te m pF ro m Position ;//指定要读取数据在文件中的起始位置
Te m pF ro mPosition =0;
XTF_HEAD 1*Te m pXTF_Head1;//XTF_HEAD 1为XTF 文件标题记录1的结构
//读标题记录1
Te m p XTF_Head1=(XTF_HEAD 1*)((char *)m_lpvF ile +Te m pFro mPosition);//m_lpvF ile 为M ap V ie wO f F ile()的返回值,这里指经过内存映射的XTF 文件
m _XTF_Head1=*Te m pXTF_Head1;Te m pF ro mPosition +
=XTF _RECORD -LENGT H;//下一个标题记录的起始位置
//读标题记录2
XTF_HEAD 2*Te m pXTF_Head2;
Te m p XTF_Head2=(XTF_HEAD 2*)((char *)m _lpvF ile +Te m pFro mPosition);
m _XTF_Head2=*Te m pXTF_Head2;
……
return TRUE ;
}
经过实验证明,当进行了内存映射后,读取大数据量文件是非常快速方便的,不需要考虑更多细节,就象利用数组一样简单.
2 利用内存DC 实现大数据量复杂测井图
形的无闪烁绘制
在开发测井解释系统时,图形元素的类型很多,经常要添加新的图形元素或者删除旧的图形元素,并且经常需要改变图形元素的属性.例如:添加或删除曲线,添加或删除道,添加或删除解释结论等,修改曲线的线型、颜色,改变道的宽度等.为了及时反映这些变化,一般要通过快速地反复刷新屏幕作图区来实现.如果采用常规绘图方法而不加特殊处理的话,屏幕上常常会出现图形闪烁的现象.由于测井图形经常是巨幅大图,即使采用分屏显示,仍然存在屏幕闪烁、停顿等现象,这种现象会极大的影响程序的运行效果,所以必须加以解决.实际软件开发中采取了先用内存DC 绘图,然后拷贝到屏幕的方法来解决这些问题.
2.1 W i n dow s 的屏幕刷新机制
在W indo w s 中,应用程序在下面情况下会发出WM _PAI NT 消息请求绘图
[22]
:
(1)当创建窗体或改变窗体大小时;
(2)在程序中使用Inva lidate 等函数来指定更新区域时;
(3)执行其它会影响客户区的操作时;
W ind w os 接收WM _PAI NT 消息后,就会执行OnPa int 函数进行窗口的重绘工作,W indo w s 自身
3
7
地球物理学进展20卷
所提供的绘图机制中,重绘窗口分两步进行:首先, W indo w s要发送WM_ERASEBKGND消息到窗口来通知窗口擦除背景,当窗口接收到WM_ERASE-BKGND消息后,会调用OnE raseBkgnd函数来进行处理,这会导致用窗口的背景色填充整个窗口.上述背景重绘之后,接着实现窗口图形的重绘.如果不断有WM_PAI NT消息请求绘图,就总会有窗口背景重绘和图形重绘交替出现.
2.2屏幕图形闪烁产生的原因分析
首先,屏幕背景刷新会引起闪烁现象[22—24].当创建窗体或改变窗体大小时,窗体背景会被擦除并重画窗体上的图形元素.此时由系统指定作图区域,并且应用程序发出WM_PAI NT消息,在该消息所对应的窗口处理函数中,先由指定的背景刷子来擦除背景,再重画窗体上的图形元素.有些时候,可能反复触发WM_PAI NT消息.例如当测井绘图井段很长,连续滚动测井图时,就会反复触发WM_ PA I NT消息,通知窗体擦除窗体背景并重画窗体上的测井图形元素.一般情况下,缺省的背景刷子的颜色不是透明的而是有色的,例如白色,那么在一次刷新过程中,用户就会看到绘图区先被白色刷子刷为白色,然后又被画上了测井图形.由于要连续滚动显示测井图,那么就要反复进行刷新,所以用户就会看到屏幕上交替显示着擦除背景后的白色背景和程序绘制出的测井图形.由于这种交替显示的频率较高,从而使用户感觉屏幕闪烁.
其次,测井图形刷新造成屏幕闪烁现象[22—24].即使克服了屏幕背景刷新引起的闪烁现象,屏幕闪烁也仍然可能存在.测井图形复杂,并且图幅巨大,当要在窗体的客户区画多个测井图形元素时,一般的绘图方法是把图形元素一个一个地依次画在客户区上,这种方法是以图形元素为单位来画的.由于后画的图形元素会覆盖先画的元素,所以在用户看来,先画的图形元素处于下方,后画的图形元素在上方.对于一些操作,需要反复刷新客户区.例如当需要连续滚动观看测井图形时,需要反复刷新客户区以跟踪反映图形状态的变化,给用户以连续的感觉,特别当测井图形曲线较多,属性复杂,并且采用直接向屏幕绘图,一次绘完整个井段时,就会出现闪烁现象.
2.3屏幕图形闪烁解决方法
在开发测井解释系统中进行绘图时,上面两类闪烁现象是同时存在的.即:要想完全解决测井绘图屏幕闪烁问题,需要同时解决屏幕背景刷新引起的闪烁现象和测井图形刷新造成屏幕闪烁现象.单独
解决其中任何一种,都不能完全解决测井图形的闪烁问题.
(1)针对背景刷新引起的闪烁,采用了不进行背景重绘的方法来消除.
由于测井绘图本身有背景颜色,所以没有必要让系统每次刷新都要按系统默认方式再去擦除和重绘背景,也就避免了擦除背景引起的闪烁.至于测井图形本身背景颜色的重绘,可以把其看作图形元素来绘制,在解决了重绘图形元素引起的闪烁时,自然就不存在绘制测井图背景引起的闪烁了.针对前面的分析,可以采取如下方法实现不执行系统默认的擦除背景操作:
①把系统缺省的背景刷子设置为NULL,那么窗体就不会在刷新屏幕时擦除背景.可以在重载后窗体的Pre C rea te W indo w函数中实现[25].
②重载O nE raseB kgnd函数,并在实现该函数时不调用基类函数.
③把默认的背景刷修改为NULL,也可以实现不擦除背景,具体如下:
void CG eo W ell V ie w::O nIn itia l Update()
{
CS croll V iew::O nIn itia l Upda te();
SetCurren tScrollS ize();
//把系统默认的背景刷子修改为NULL
SetC la ssLong(th is->m_h W nd,GCL_HBR-BACKGROUND,(LONG)((HBR US H) (NULL)));
}
在开发测井解释系统G eo W ell中,采用了第三种方式,效果良好.
(2)针对测井图形刷新造成的屏幕闪烁现象,采用内存DC绘图来消除屏幕闪烁.
这类闪烁现象是由于测井图形刷新时屏幕是以元素为单位,把测井图形元素逐个画出来而造成的.当以屏幕为单位进行刷新,并且速度较快时,就可以消除闪烁,而在内存DC中绘图,然后拷贝到屏幕的方法恰恰可以实现这样的功能.具体是:把测井图形先用内存设备环境DC绘制在与显示兼容的位图中,然后从内存环境把测井图形复制到屏幕客户区,这样,在测井图形绘制到屏幕之前,已经将图形绘制在位图中,然后直接复制到屏幕上,实现以屏幕为单位绘图,从而消除了闪烁.在一些文献中介绍了这种方法[22-24],但是,由于测井绘图的特殊性,使得在具
47
1期李河,等:大数据量地球物理测井绘图关键技术研究
体开发中可能遇到问题.很多公开的文献中,屏幕设备环境和内存设备环境都采用设备坐标就能够满足要求,但是在测井绘图中,需要采用逻辑坐标,如果不反复实验,可能会因为所采用的映射方式不同于普通映射方式,在一些函数中不能正确传递参数,而使程序无法运行.下面给出一些重要函数的参数传递的实例:
/*p S creen DC为屏幕DC,p M e m ory DC为兼容的内存DC,Current C lien t S creen V ie w Rect 为当前的客户区坐标(设备坐标,像素),Current S croll Position Poin t为当前滚动条位置(当前映射模式下滚动条位置的逻辑坐标)*/
BOOLC D raw P lot::D raw A ll Tracks(CDC*p S creen DC,CDC*p M e m ory DC,C R ect Current C lien t Screen V ie w Rect,C Point C urrent S croll Position Poin t)
{
……
//准备内存D C和屏幕D C
C B it m ap M e m ory
D C B it m ap;//位图对象
p S creen D C->Set M ap M ode(MM_LO-METR I C);//设置屏幕坐标映射模式
p M e m ory D C->C rea te C o m pa tible DC(p S creen DC);//设置与屏幕DC兼容的内存DC p M e m ory D C->Set M ap M ode(MM_LO-METR I C);//设置内存DC坐标映射模式
//创建与屏幕DC兼容的位图对象
M e m ory DC B it m ap.C rea te C o m patible B it m-ap(p S creen DC,C urrent C lient Screen V iew Re-ct.W idth(),Current C lient S creen V ie w R ect. He ight());
//把新位图选入内存DC
C B it m ap*PO ld B it m ap=p M e m ory DC-> Se lect Ob ject(&M e m ory DC B it m ap);
C Point Current C lient Screen V ie w R ight C orner Point(Curren t C lient Screen V ie w R ect. righ t,Cu rrent C lien t S creen V ie w R ect.botto m);
//设备坐标转换为逻辑坐标
p S creen DC->D P to LP(&Current C lient S creen V ie w R ight Corner Point);
floa t Curren t C lient S creen V ie w W idth;
floa t Curren t C lient S creen V ie w He ight;
Current C lient Screen V ie w W idth=(floa t)Current C lient S creen V ie w R ight C orner Po int.x;
Current C lient Screen V ie w He ight=(floa t) (-C urrent C lien t Screen V ie w R ight C orner Po i n t.y);
//填充测井图形背景
p M e m ory DC->F ill Solid Rect(0,Current C lien t Screen V ie w R ight C orner Po int.y,Current C lien t S creen V ie w R ight C orner Point.x,-Cur-rent C lien t S creen V ie w R igh t C orner Po int.y,m_ P lot Pa ram eter.B ack C olor);
//准备完毕,可以开始在内存DC中绘图
……
//在内存中把测井图绘制完毕,拷贝到屏幕
p S creen DC->B it B lt(
Current S croll Position Poin t.x,
Current S croll Position Po int.y+Curren t C li-ent S creen V iew R igh t Corner Poin t.y,
Current C lient Screen V iew R ight C orner Po i n t.x,
-Current C lient S creen V iew R ight C orner Po i n t.y,
p M e m ory DC,0,Current C lien t S creen V ie w R ight C orner Po int.y,SRCCOPY);
p M e m ory DC->Se lect Ob ject(PO ld B it m-ap);
p M e m ory DC->D elete DC();
M e m ory DC B it m ap.D elete Ob ject();
return TRUE;
}
实际开发中,两类闪烁现象同时存在,同时采用上面两种技术,解决了测井绘图中的屏幕闪烁问题,获得了很好的效果.
3利用局部重绘法实现大数据量测井图形的快速连续滚动显示
测井图形一般幅度很大,当道很多并且绘图深度很长时,横向和纵向都不可能在当前窗口一次全部显示出来.如果采用通常的绘图方法,一次把整幅测井图全部绘制出来,由于绘图数据量非常大,绘图速度很慢,甚至用户无法接受.所以必须采用好的绘图技术,提高绘图效率.连续滚动显示是常用的方法,但是必须解决巨幅测井图形的绘图速度问题. 3.1整幅图形全刷新的方法实现滚动存在的问题
分析
57
地 球 物 理 学 进 展
20卷
为了解决屏幕闪烁问题,前面采用了内存DC
绘图的技术.那么,如果采用一般的、也是最容易实现的滚动显示方案,就是每次滚动都重新绘制整幅图形,即:全刷新的方法
[26]
.具体是:先根据实际测
井图形的大小,在内存中把整幅图都绘制出来,每次滚动时,都重新绘制整幅图形,然后根据采用的映射方式,把整幅图拷贝到屏幕.现假设当前映射方式下X 坐标向右为正,Y 坐标向上为正,那么不管当前滚动条的位置在哪儿,只要把内存中的图形左上角对应屏幕坐标的(0,0)点进行拷贝,就可以实现滚动显示.如果绘制的图形幅度较小,这种方法实现比较容易,也能满足要求.但是测井图形一般很大,每次都绘整幅图形,绘图速度会很慢,而且由于图形大,开辟的内存也大,会出现严重的内存不足问题.所以实际开发中一般不采用.
3.2 只刷新当前屏幕的方法实现滚动显示
只刷新当前屏幕的方法是指每次滚动都根据当前滚动条的位置,计算屏幕应该显示图形的哪部分,然后在内存中绘出图形,并拷贝到当前屏幕客户区.这种方法每次只需要绘制当前屏幕客户区一样大小的图形,这样一方面减少了每次滚动时绘图的数据量,提高了绘图速度,同时由于只绘制和屏幕一样大小的图形,所以占用的内存也小了,提高了绘图效率.具体如图2
所示.
图2 屏幕滚动原理F ig 2 S creen scrolling theory
图2(a)为滚动前屏幕的状态,屏幕图形由B 和C 两部分组成,其中B 部分的高度为h 1,C 部分的高度为h 2.如果每次滚动的纵向高度为h 2,现在向下滚动,则滚动一次后屏幕状态如图2(b),屏幕图形由A 和B 两部分组成,原来的C 部分滚出了屏幕.
由于采用刷新当前屏幕的方法实现滚动,则A 和B 两部分作为一个整体在内存中绘制好后,拷贝到屏幕.其它方向的滚动实现类似.如果连续滚动,则重复上述操作即可实现.这种方法比整幅图形刷新效率高很多,但是,每次刷新必须准确计算当前屏幕图形在整幅图形中的位置.
3.3 局部重绘的方法实现快速连续滚动显示
只刷新当前屏幕的方法在一些情况下可以应用,但是当测井图形复杂,数据量很大时,仍然存在绘图速度慢的问题,如绘制声电成像测井图.而采用局部重绘的方法每次只绘制新滚入屏幕的部分,可以进一步减少绘图工作量,提高绘图速度,这种方法由于每次刷新重绘的面积小,绘图占用的内存也少,具体如图2所示.
当滚动一次以后,B 部分的内容并未变化,只是在窗口中的显示位置发生了改变.因此,可以不必重新绘制B 部分,只要把B 部分拷贝后向下移动h 2距离,然后在内存中绘制A 部分,并拷贝到屏幕上,具体位置如图2(b)所示.此时原来屏幕的C 部分已经移出屏幕,不必再绘了.这样绘图数据量和内存占用都少了,提高了绘图效率.其它方向的滚动原理类似.重复上述工作,即可以实现连续快速滚动显示.开发G eo W e ll 系统时,采用了这种方法.
4 结 论
针对开发测井解释系统绘图过程中遇到的主要问题,本文给出了相应的解决办法.当然,在开发测井绘图程序过程中,不仅涉及上面几个主要技术,还有其它方面,例如:在曲线出道处理中,采用的处理技术不同,绘图的效率也就不同.我们在开发中采用了计算机图形学所涉及的剪裁技术,通过设置剪裁区域,实现曲线的出道处理.在常规曲线岩性填充中,采用了区域技术,通过区域运算,确定实际需要填充的区域,实现了各类岩性填充.在通过内存DC 实现绘图时,也利用了M FC 本身的优势,只在内存中绘制剪裁区内需要更新的区域,大大提高了绘图
效率,特别是在声电成像图像的绘制中,效果更明显.这些技术作者另有文章详细介绍.本文介绍的几个技术是作者在开发中采用的几个关键技术,希望对今后测井绘图有所帮助.参 考 文 献(R eferences):
[1] 金勇,唐文清,陈福利,等.石油测井地质综合应用网络平台
For https://www.doczj.com/doc/402244372.html,[J].石油勘探与开发,2004,31(3):92～96.
6
7
1期李河,等:大数据量地球物理测井绘图关键技术研究
[2]赵连功,刘洪.地震勘探数据资料处理软件集成化研究现状和
发展趋势[J].地球物理学进展,2003,18(4):598～601.
[3]刘光鼎,祝靓谊.近期油气勘探地球物理的一些新进展[J].地
球物理学进展,2003,18(3):363～367.
[4]张建中,陈世军,徐初伟.动态网络最短路径射线追踪[J].地球
物理学报,2004,47(5):899～904.
[5]张荣,张迁.一种新的地磁图纸数字化软件GDDS[J].地球物
理学进展,2004,19(1):183～185.
[6]崔若飞,李晋平.地震资料小波剖面制作系统[J].地球物理学
进展,2003,18(3):562～566.
[7]门少华,王祝文,梅忠武.单机测井解释系统的开发[J].长春科
技大学学报,2001,31(4):398～403.
[8]陶果,多雪峰.我国地球物理测井技术的发展与战略初探[J].
地球物理学进展,2001,16(3):47～50.
[9]胡恒山,刘家琦,王洪滨,等.基于简化的Pride理论模拟声电
效应测井响应[J].地球物理学报,2003,46(2):259～264. [10]刘四新,佟文琪.电磁波测井的现状和发展趋势[J].地球物理
学进展,2004,19(2):235～237.
[11]赵爱华,张中杰.三维复杂介质中转换波走时快速计算[J].地
球物理学报,2004,47(4):702～707.
[12]熊晓东,郭麦成,袁洪清,等.多功能测井数据处理软件的开发
[J].江汉石油学院学报,2002,24(3):76～77.
[13]黄建国,陈东.卫星介质深层充电的计算机模拟研究[J].地球
物理学报,2004,47(3):392～397.
[14]白彦彬,杨长春,袁金敏.面向对象的测井解释软件平台
(L I SP)开发[J].计算机工程与应用,2002,(8):243～245.[15]范翔宇,夏宏泉,陈平,等.井眼轨迹和地层剖面及测井曲线综
合成图研究[J].西南石油学院学报,2004,26(2):1～3. [16]吕京国,黄国满,杨明辉.用V isua l C++实现大数据量的快
速存取[J].测绘科学,2002,27(3):29～32.
[17]魏代永,郝向阳,巩学美.地图矢量数据计算机内存管理方法
[J].解放军测绘学院学,1999,16(2):109～112.
[18]邵文秀,程明霄.基于内存映射文件的数据交换技术及其在测
控软件中的应用[J].现代电子技术,2004,(1):36～38
[19]徐光伟,顾景文,阎福礼,等.基于内存映像文件的网络信息安
全性存取的探讨[J].计算机应用研究,2002,(7):48～53. [20]杨宁学,诸昌钤,聂爱丽.内存映射文件及其在大数据量文件
快速存取中的应用[J].计算机应用研究,2004,(8):187～
188.
[21]马玲华,杨劲松,李传伟,等.ECL I PS5700测井系统XTF文
件格式分析[J].测井技术,2001,25(3):225～230.
[22]柏瑞,周赣.图形系统屏幕刷新技术及快速重画方法研究[J].
计算机工程与应用,2003,(16):136～138.
[23]胡还刚,龚晖.多缓冲技术的一种实现方法及其应用[J].计算
机应用研究,2003,(4):20～21.
[24]齐邦强.W indo w s编程中的图形刷新问题[J].计算机应用,
2002,22(4):103～104.
[25]王洁宁,侯启真.应用离屏DC实现窗口无闪烁图像显示技术
[J].中国民航学院学报,2001,(19)3:38～41.
[26]张福明,邵才瑞,李洪奇.X-W indo w环境下图形连续滚动显
示方法设计[J].计算机工程,1998,24(9):68～70.
77
测井曲线代码一览表 测井类资料２0０9－0８—0７16:0１阅读437 评论0 字号: 大大中中小小 from石油科技论坛 常用测井曲线名称 测井符号英文名称中文名称 Rt tｒue formatiｏn resistivitｙ。地层真电阻率 Rxo flushed ｚｏne ｆoｒｍａｔiｏn rｅｓiｓｔivｉty 冲洗带地层电阻率 Ｉld ｄeep investigate induｃtiｏn log 深探测感应测井 Ｉlｍ meｄium investigaｔe inｄｕction log 中探测感应测井 Ils shallｏｗinｖｅstｉgatｅinducｔioｎｌog 浅探测感应测井 Rd deｅｐiｎveｓtigate douｂle lateral reｓｉstiｖｉty log 深双侧向电阻率测井 Ｒs sｈaｌｌow inｖesｔｉgate double ｌaｔerａl ｒesistｉvity ｌog 浅双侧向电阻率测井 ＲＭLＬ micｒｏlateral rｅsistｉviｔyｌog 微侧向电阻率测井COＮ induｃtｉｏn ｌoｇ感应测井 ＡC ａcoustiｃ声波时差 DEN density 密度 CN neutron 中子 GＲ natｕral ｇａmmａｒaｙ自然伽马 SP spontaneous poｔｅntial 自然电位 CAL ｂorehole ｄiａmeter 井径 K potassium 钾 TH thoriuｍ钍 U ｕraniuｍ铀 KTＨ gammａray wiｔhout uraｎiuｍ无铀伽马 NGR neｕtroｎgaｍma ｒａy 中子伽马 5700系列得测井项目及曲线名称 StaｒIｍａｇeｒ微电阻率扫描成像 CBIL 井周声波成像 MAC 多极阵列声波成像 MRIＬ核磁共振成像 TBRT 薄层电阻率 DAC 阵列声波 DVRT 数字垂直测井 HDIP 六臂倾角 MＰHI 核磁共振有效孔隙度 MBVM 可动流体体积 ＭBＶI 束缚流体体积 MPERM 核磁共振渗透率 Ｅchoｅs 标准回波数据
第五单元精确绘图 说明：为了快速地精确地绘制图形，本单元在绘图中均使用“命令化名”来扫行绘图和编辑命令，同时也希望同学们能快速掌握执行命令来画图。 5-1 第1题 【操作要求】 1，建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区内。 2．绘图：按图5-1规定的尺寸绘图，中心线型为acad iso12W120，调整线型比例。3．保存：将完成的图形以KSCAD5-1.DWG为文件名保存在考生文件夹中。 图5-1 绘制图形 【操作步骤】 第1步：绘图前的设置工作。 1．设置图形界限。选择【格式】/【图形界限】菜单命令，回车（接受默认的起始点）；输入：@150，150回车，完成图形界限的设置。 提示：如果单击【视图】工具条上的“俯视图”按钮，可将所设置的图形界限最大化。这样，可确保所绘图形在所设置的图形界限之内。 2．设置图形及线型。选择【格式】/【图层】菜单命令，在打开的【图层特性管理器】对话框中，单击按钮，新建二个图层。第一个图层，图层名为：中心线；颜色设置为“红色”,线宽默认；单击加载acad iso12W120（双点划线）线型，线宽为默认；单击按钮,将”全局比例因子”设置为:0.2。第二个图层，图层名为：尺寸标注，颜色、线型和线宽均默认。将“0层”的线宽设置为：0.2后，单击按钮，完成设置。如图5-1-1所示。 图5-1-1图层特性的设置 第2步：绘制定位线。 画垂直中心线。将“中心线层”置为：“当前层”，绘制二条相互垂直的中心线；2．画定位线。在命令行输入：O（“偏移”命令），输入：56，回车；选择水平中心线向下单击。回车继续执行“O”命令，输入：10，回车，选择垂直中心线，分别向左右单击，绘制的定位线，如图5-1-2所示。 图5-1-2 画定位中心线图5-1-3 画圆
测井曲线代码 RD、RS—深、浅侧向电阻率 RDC、RSC—环境校正后的深、浅侧向电阻率VRD、VRS—垂直校正后的深、浅侧向电阻率DEN—密度 DENC—环境校正后的密度 VDEN—垂直校正后的密度 CNL—补偿中子 CNC—环境校正后的补偿中子 VCNL—垂直校正后的补偿中子 GR—自然伽马 GRC—环境校正后的自然伽马 VGR—垂直校正后的自然伽马 AC—声波 V AC—垂直校正后声波 PE—有效光电吸收截面指数 VPE—垂直校正后的有效光电吸收截面指数SP—自然电位 VSP—垂直校正后的自然电位 CAL—井径 VCAL—垂直校正后井径 KTh—无铀伽马 GRSL—能谱自然伽马 U—铀 Th—钍 K—钾 WCCL—磁性定位 TGCN—套管中子 TGGR—套管伽马 R25—2.5米底部梯度电阻率 VR25—环境校正后的2.5米底部梯度电阻率DEV—井斜角 AZIM—井斜方位角 TEM—井温 RM—井筒钻井液电阻率 POR2—次生孔隙度 POR—孔隙度 PORW—含水孔隙度 PORF—冲洗带含水孔隙度 PORT—总孔隙度 PERM—渗透率 SW－含水饱和度 SXO—冲洗带含水饱和度
SH—泥质含量 CAL0—井径差值 HF—累计烃米数 PF—累计孔隙米数 DGA—视颗粒密度 SAND，LIME，DOLM，OTHR—分别为砂岩，石灰岩，白云岩，硬石膏含量 VPO2—垂直校正次生孔隙度 VPOR—垂直校正孔隙度 VPOW—垂直校正含水孔隙度 VPOF—垂直校正冲洗带含水孔隙度 VPOT—垂直校正总孔隙度 VPEM—垂直校正渗透率 VSW－垂直校正含水饱和度 VSXO—垂直校正冲洗带含水饱和度 VSH—垂直校正泥质含量 VCAO—垂直校正井径差值 VDGA—垂直校正视颗粒密度 VSAN，VLIM，VDOL，VOTH—分别为垂直校正砂岩，石灰岩，白云岩，硬石膏含量岩石力学参数 PFD1—破裂压力梯度 POFG—上覆压力梯度 PORG—地层压力梯度 POIS—泊松比 TOUR—固有剪切强度 UR—单轴抗压强度 YMOD—杨氏模量 SMOD—切变模量 BMOD—体积弹性模量 CB—体积压缩系数 BULK—出砂指数 MAC MAC—偶极子阵列声波 XMAC-Ⅱ—交叉偶极子阵列声波 DTC1—纵波时差 DTS1—横波时差 DTST1—斯通利波时差 DTSDTC-纵横波速度比 TFWV10-单极子全波列波形 TXXWV10-XX偶极子波形 TXYWV10- XY偶极子波形 TYXWV10- YX偶极子波形 TYYWV10- YY偶极子波形 WDST-计算各向异性开窗时间 WEND-计算各向异性关窗时间
测井符号英文名称中文名称 Rt trueformation resistivity. 地层真电阻率 Rxo flushed zone formation resistivity 冲洗带地层电阻率 Ild deep investigate induction log 深探测感应测井 Ilm medium investigate induction log 中探测感应测井 Ils shallow investigate induction log 浅探测感应测井 Rd deep investigate double lateral resistivity log深双侧向电阻率测井Rs shallow investigate double lateral resistivity log浅双侧向电阻率测井 RMLL micro lateral resistivity log 微侧向电阻率测井 CON induction log 感应测井 AC acoustic 声波时差 DEN density 密度CN neutron 中子GR natural gamma ray 自然伽马SP spontaneous potential 自然电位CAL borehole diameter 井径K potassium 钾 TH thorium 钍U uranium 铀KTH gamma ray without uranium 无铀伽马NGR neutron gamma ray 中子伽马 5700系列的测井项目及曲线名称 Star Imager 微电阻率扫描成像 CBIL 井周声波成像 MAC 多极阵列声波成像 MRIL 核磁共振成像 TBRT 薄层电阻率 DAC 阵列声波 DVRT 数字垂直测井 HDIP 六臂倾角
第五单元精确绘图 5.1第1题（机械类） 【操作要求】 1．建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区域内。 2．绘图：按图5-1规定的尺寸绘图，中心线线型为acad_iso12w100，调整线型比例。 3．保存：将完成的图形以KSCAD5-1.DWG为文件名保存在考生文件夹中。
5.2第2题（机械类） 【操作要求】 1．建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区域内。 2．绘图：按图5-2规定的尺寸绘图，中心线线型为acad_iso10w100，调整线型比例。 3．保存：将完成的图形以KSCAD5-2.DWG为文件名保存在考生文件夹中。
5.3第3题（机械类） 【操作要求】 1．建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区域内。 2．绘图：按图5-3规定的尺寸绘图，中心线线型为CENTER，点线圆的线型为acad_iso07w100，调整线型比例。 3．保存：将完成的图形以KSCAD5-3.DWG为文件名保存在考生文件夹中。
【操作要求】 1．建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区域内。 2．绘图：按图5-4规定的尺寸绘图，中心线线型为CENTER，调整线型比例。 3．保存：将完成的图形以KSCAD5-4.DWG为文件名保存在考生文件夹中。
【操作要求】 1．建立绘图区域：建立合适的绘图区域，图形必须在设置的绘图区域内。 2．绘图：按图5-5规定的尺寸绘图，中心线线型为CENTER，调整线型比例。 3．保存：将完成的图形以KSCAD5-5.DWG为文件名保存在考生文件夹中。
各种测井曲线代码 附录33 测井曲线名称代码 名称代码名称代码名称代码 0、4米电位电阻率 R04 井径1 C1 阵列感应4英尺分辨率及60英寸探测深度电阻率 AF60 0、45米电位电阻率 R045 井径2 C2 阵列感应4英尺分辨率及90英寸探测深度电阻率 AF90 0、5米电位电阻率 R05 井径3 C3 阵列感应4英尺分辨率侵入带真电阻率 AFRX 1米底部梯度电阻率 R1 井斜 DEV 补偿声波时差 AC 2、5米底部梯度电阻率 R25 井斜方位 AZIM 井径CAL 4米底部梯度电阻率 R4 高分辨率侧向电阻率 LLHR 长源距声波时差 DT 6米底部梯度电阻率 R6 方位电阻率曲线1 ARO1 纵横波速度比 VPVS 8米底部梯度电阻率 R8 方位电阻率曲线10 AR10 纵横波方式单极横波时差 DT4S 深侧向电阻率 RD 方位电阻率曲线11 AR11 纵横波方式单极纵波时差 DT4P 浅侧向电阻率 RS 方位电阻率曲线12 AR12 泊松比PR 邻近侧向电阻率 RPRX 方位电阻率曲线2 ARO2 上偶极横波时差 DT2 微侧向电阻率 RMLL 方位电阻率曲线3 ARO3 下偶极横波时差 DT1 微球型聚焦电阻率 MSFL 方位电阻率曲线4 ARO4 斯通利波时差 DTST 深感应电阻率 RILD 方位电阻率曲线5 ARO5 全波列波形 WF
中感应电阻率 RILM 方位电阻率曲线6 ARO6 声波成象ACI 八侧向电阻率 RFOC 方位电阻率曲线7 ARO7 自然伽马GR 球型聚焦电阻率 SFLU 方位电阻率曲线8 ARO8 无铀自然伽马 CGR 数字聚焦电阻率 DFL 方位电阻率曲线9 ARO9 钾 K 感应电导率 COND 阵列感应1英尺分辨率地层真电阻率AORT 钍 TH 微电位电阻率 ML1 阵列感应1英尺分辨率及10英寸探测深度电阻率 AO10 铀 U 微梯度电阻率 ML2 阵列感应1英尺分辨率及20英寸探测深度电阻率 AO20 补偿中子 CNL 钻井液电阻率 RM 阵列感应1英尺分辨率及30英寸探测深度电阻率 AO30 井壁中子 SNL 井温 TEMP 阵列感应1英尺分辨率及60英寸探测深度电阻率AO60 中子伽马 NGR 钻头直径 BS 阵列感应1英尺分辨率及90英寸探测深度电阻率 AO90 补偿密度 DEN 200兆赫兹电阻率 R4SL 阵列感应1英尺分辨率侵入带真电阻率 AORX 岩性密度 LDL 200兆赫兹幅度比 R4AT 阵列感应2英尺分辨率地层真电阻率ATRT 密度校正值 DRH 200兆赫兹介电常数 D2EC 阵列感应2英尺分辨率及10英寸探测深度电阻率 AT10 光电吸收截面指数 PE 200兆赫兹相位角 P2HS 阵列感应2英尺分辨率及20英寸探测深度电阻率 AT20 核磁共振总孔隙度 TPOR 47兆赫兹电阻率 R4SL 阵列感应2英尺分辨率及30英寸探测深度电阻率 AT30 核磁共振渗透率 KCMR 47兆赫兹幅度比 R4AT 阵列感应2英尺分辨率及60英寸探测深度电阻率 AT60 核磁共振束缚流体体积 MBVI
本次课主要讲两部分内容： 第四章精确绘图 第六章图形显示控制 教学目的： 1、掌握对象捕捉功能和绘图辅助工具，以便按尺寸精确作图。 2、掌握图形显示控制功能。 教学重点： 二维绘图坐标系、对象捕捉、图形显示控制。 教学难点： 对象捕捉；显示缩放命令的参数。 第四章精确绘图 一、二维绘图坐标系 AutoCAD中的二维坐标系有世界坐标系（WCS）和用户坐标系（UCS）两种形式。它们是根据右手笛卡尔坐标系定义的。 1、种类： （1）世界坐标系WCS： WCS是AutoCAD默认的绘图坐标系。 WCS的坐标原点位于图纸左下角，X轴为水平轴，向右为正； Y轴为垂直轴，向上为正；Z轴垂直于XY平面，指向绘图者为正。 WCS不能被改变。
WCS图标在原点WCS图标不在原点（2）用户坐标系UCS： UCS是用户自己相对于WCS建立起来的坐标系。 用户坐标系可用UCS命令来创建。 UCS图标在原点UCS图标不在原点 2、UCS命令 功能：改变坐标系的原点和方向，方便作图。 调用： （1）定义UCS的原点 * 单击“UCS”工具栏中图标； * 下拉菜单：工具→新建UCS→原点； * 在命令行中输入“UCS” 命令,选择M参数。 （2）原点不变，UCS绕Z轴旋转 * 单击“UCS”工具栏中图标； * 下拉菜单：工具→新建UCS→ Z； * 在命令行中输入“UCS” 命令, 选择N参数，再选择Z参数,指 定旋转角度。 （3）选定一个对象来定义新的坐标系 * 单击“UCS”工具栏中图标； * 下拉菜单：工具→新建UCS→对象； * 在命令行中输入“UCS” 命令,选择N参数，指定对象。 新建UCS的原点和X轴方向根据对象和拾取点决定。 （4）恢复世界坐标系 * 单击“UCS”工具栏中图标； * 下拉菜单：工具→新建UCS→世界； * 在命令行中输入“UCS” 命令,直接回车。 二、对象捕捉 1、基本设置 调出SP78图4-20“选项”对话框： * 在状态栏的“Osnap”按钮上单击鼠标右键，弹出快捷菜单，选“Settings…”→单击“选项”（Options…）按钮→Options对话框→Drafting 选项卡进行基本设置。
原始测井曲线代码 代码名称 A1R1 T1R1声波幅度 A1R2 T1R2声波幅度 A2R1 T2R1声波幅度 A2R2 T2R2声波幅度AAC 声波附加值 AA VG 第一扇区平均值AC 声波时差 AF10 阵列感应电阻率AF20 阵列感应电阻率AF30 阵列感应电阻率AF60 阵列感应电阻率AF90 阵列感应电阻率AFRT 阵列感应电阻率AFRX 阵列感应电阻率AIMP 声阻抗 AIPD 密度孔隙度 AIPN 中子孔隙度 AMA V 声幅 AMAX 最大声幅 AMIN 最小声幅 AMP1 第一扇区的声幅值AMP2 第二扇区的声幅值AMP3 第三扇区的声幅值AMP4 第四扇区的声幅值AMP5 第五扇区的声幅值AMP6 第六扇区的声幅值AMVG 平均声幅 AO10 阵列感应电阻率AO20 阵列感应电阻率AO30 阵列感应电阻率AO60 阵列感应电阻率AO90 阵列感应电阻率AOFF 截止值 AORT 阵列感应电阻率AORX 阵列感应电阻率APLC 补偿中子 AR10 方位电阻率 AR11 方位电阻率 AR12 方位电阻率 ARO1 方位电阻率 ARO2 方位电阻率 ARO3 方位电阻率ARO4 方位电阻率 ARO5 方位电阻率 ARO6 方位电阻率 ARO7 方位电阻率 ARO8 方位电阻率 ARO9 方位电阻率 AT10 阵列感应电阻率 AT20 阵列感应电阻率 AT30 阵列感应电阻率 AT60 阵列感应电阻率 AT90 阵列感应电阻率 ATA V 平均衰减率 ATC1 声波衰减率 ATC2 声波衰减率 ATC3 声波衰减率 ATC4 声波衰减率 ATC5 声波衰减率 ATC6 声波衰减率 ATMN 最小衰减率 ATR T 阵列感应电阻率 ATRX 阵列感应电阻率 AZ 1号极板方位 AZ1 1号极板方位 AZI 1号极板方位 AZIM 井斜方位 BGF 远探头背景计数率 BGN 近探头背景计数率 BHTA 声波传播时间数据 BHTT 声波幅度数据 BLKC 块数 BS 钻头直径 BTNS 极板原始数据 C1 井径 C2 井径 C3 井径 CAL 井径 CAL1 井径 CAL2 井径 CALI 井径 CALS 井径 CASI 钙硅比 CBL 声波幅度 CCL 磁性定位 CEMC 水泥图 CGR 自然伽马 CI 总能谱比 CMFF 核磁共振自由流体体积 CMRP 核磁共振有效孔隙度 CN 补偿中子 CNL 补偿中子 CO 碳氧比 CON1 感应电导率 COND 感应电导率 CORR 密度校正值 D2EC 200兆赫兹介电常数 D4EC 47兆赫兹介电常数 DAZ 井斜方位 DCNT 数据计数 DEN 补偿密度 DEN_1 岩性密度 DEPTH 测量深度 DEV 井斜 DEVI 井斜 DFL 数字聚焦电阻率 DIA1 井径 DIA2 井径 DIA3 井径 DIFF 核磁差谱 DIP1 地层倾角微电导率曲线1 DIP1_1 极板倾角曲线 DIP2 地层倾角微电导率曲线2 DIP2_1 极板倾角曲线 DIP3 地层倾角微电导率曲线3 DIP3_1 极板倾角曲线 DIP4 地层倾角微电导率曲线4 DIP4_1 极板倾角曲线 DIP5 极板倾角曲线 DIP6 极板倾角曲线 DRH 密度校正值 DRHO 密度校正值 DT 声波时差 DT1 下偶极横波时差 DT2 上偶极横波时差 DT4P 纵横波方式单极纵波时 差 DT4S 纵横波方式单极横波时 差 DTL 声波时差
第一章自然电位测井（SP） 一.分析自然电位的成因，写出扩散电动势、扩散吸附电动势、总电动势表达式。 答：（1）井内自然电位产生的原因：对于油井来说，主要有两个原因，1）地层水含盐浓度和泥浆含盐浓度不同，引起离子的扩散作用和岩石颗粒对离子的吸附作用；2）地层压力与泥浆柱压力不同时，在地层孔隙中产生过滤作用。在扩散过程中，正、负离子迁移率（速度）不同，通常是负离子快，这样在某一时刻通过同一截面的正离子数与负离子数不同，结果是浓度低的一侧形成了负离子（电荷）的富集，而浓度高的一侧形成了正 离子（电荷）的富集，从而产生了扩散电位。由于扩散吸附作用，其结果是浓度高的一侧形成了负离子（电荷）的富集，而浓度低的一侧形成了正离子（电荷）的富集，从而产生了扩散吸附电位。 （2）扩散电动势：Ed=Kd*lgCw/Cmf=Kd*lgRmf/Rm Kd- 扩散电位系数 Cw-地层水 的电化学活度 Cmf- 泥浆滤液的电化学活度 Rmf— 泥浆滤液的电阻率 Rm—地层水的电 阻率； （3）扩散吸附电动势：Ea=Ka*lgCw/Cmf=Ka*lgRmf/Rm Ka-扩散吸附电位系数 Cw- 地层水的电化学活度 Cmf- 泥浆滤液的电化学活度 Rmf— 泥浆滤液的电阻率 Rm—地 层水的电阻率； （4）总电动势：Eda=Kda*lgCw/Cmf=Kda*lgRmf/Rm Ka-扩散- 吸附电位系数 Cw- 地层水的电化学活度 Cmf- 泥浆滤液的电化学活度 Rmf— 泥浆滤液的电阻率 Rm—地层 水的电阻率： 二.不同Cw、Cmf 情况下自然电位测井曲线有哪些特征？ 答：在砂泥岩剖面中渗透层通常有自然电曲线异常现象:
浅谈CCL测井深度电缆磁记号标定方法 摘要：由于测井电缆记号标注校深操作的准确性直接关系到测井数据的准确性，我们通过对新旧标注方法的工作原理对比，总结对比新旧标定法特点。CCL电缆标定仪就是以新的方法为理论依据研制产生的，它所采用的原理是依据标准井井下套管序列为参考深度对测井电缆进行加注磁性记号。在新科技不断发展的情况下，新技术在测井校深领域的应用前景会越来越广阔。 关键词：测井电缆记号标定误差修正CCL 测井深度测量是石油测井中最为重要的参数之一，他对测井安全、测井解释、油气田后期开发都有着举足轻重的作用。在处理多井的测井数据时，校对深度的精度尤为重要，其误差对绘图、流体界面深度估算及构造解释均有较大影响。 一、问题的提出 随着油田开发的不断深化，为了得到详实有效、准确无误的井下数据，石油企业对测井作业环节提出了越来越高的要求。要求测井数据必须精准高效，从而为油田后续开发判断油层位臵、砂岩孔隙以及含油性提供科学的依据。在测井作业中，不论是岩性曲线测量、孔隙度曲线测量还是对电阻率曲线的测量，校深这一不可缺少的测井工艺发挥着重要的，不容忽视的作用，校深的准确与否直接关系着数据采集的准确性和有效性，校深的精度从一定意义上决定着一切测井数据的精度。基于以上因素，亟待需要一种测井操作简便、深度值控制
精准、数据接口兼容性好的校深仪器。测井深度电缆磁记号标定仪随即产生，由于近年被测井企业的广泛应用，其科学性和准确性也逐渐显露出来。 二、传统测井校深方法分析 传统测井电缆深度记号标注法主要有以下两种：一种是人为把记号注磁器和记号接收器之间的距离设定为25米，并以此为标尺对测井电缆进行深度记号的标定。另一种方法是以编码器待脉冲数来确定25米的距离并进行深度记号标定。 以上两种方法所标定的深度记号都会直接造成累计误差和记号大小不一的误测现象，对测井测量段的准确性的影响很大，更会将直接影响到规定测量段测井数据的采集和分析，对油气区块的分析以及后期开发十分不利。 传统的设定记号注磁器和记号接收器之间的距离为25米这种方法比较落后，它的原理等同于，用长卷尺人为的测量25米然后注一个磁记号一样。用这种方法所标定的深度记号存在着深度记号间距大小不一、误差偏大等不准确因素。而另一种通过改变编码脉冲数来确定25米长度，并以此为标准触发注磁器对测井电缆进行注磁的测深方法，普遍存在着受滑轮磨损和电缆磨损现象，造成电缆磨损粗细不均、电缆拉伸张力发生变化等不利因素造成的校深不准甚至校深失效现象。所标定的深度记号同样存在着深度记号间距大小不一，误差偏大的重大缺陷。CCL电缆深度标定仪就能很好的避免这些问题，充分保证数据的准确性，将误差率降到最小。
测井曲线代码一览表 测井类资料2009-08-07 16:01 阅读437 评论0 字号：大大中中小小 from 石油科技论坛 常用测井曲线名称 测井符号英文名称 中文名称 Rt true formation resistivity. 地层真电阻率 Rxo flushed zone formation resistivity 冲洗带地层电阻率 Ild deep investigate induction log 深探测感应测井 Ilm medium investigate induction log 中探测感应测井 Ils shallow investigate induction log 浅探测感应测井 Rd deep investigate double lateral resistivity log 深双侧向电阻率测井 Rs shallow investigate double lateral resistivity log 浅双侧向电阻率测井 RMLL micro lateral resistivity log 微侧向电阻率测井CON induction log 感应测井 AC acoustic 声波时差 DEN density 密度 CN neutron 中子 GR natural gamma ray 自然伽马 SP spontaneous potential 自然电位
CAL borehole diameter 井径 K potassium 钾 TH thorium 钍 U uranium 铀
KTH gamma ray without uranium 无铀伽马 NGR neutron gamma ray 中子伽马 5700系列的测井项目及曲线名称 Star Imager 微电阻率扫描成像 CBIL 井周声波成像 MAC 多极阵列声波成像 MRIL 核磁共振成像 TBRT 薄层电阻率 DAC 阵列声波 DVRT 数字垂直测井 HDIP 六臂倾角 MPHI 核磁共振有效孔隙度 MBVM 可动流体体积 MBVI 束缚流体体积 MPERM 核磁共振渗透率 Echoes 标准回波数据 T2 Dist T2分布数据 TPOR 总孔隙度 BHTA 声波幅度 BHTT 声波返回时间 Image DIP 图像的倾角 COMP AMP 纵波幅度 Shear AMP 横波幅度 COMP ATTN 纵波衰减 Shear ATTN 横波衰减 RADOUTR 井眼的椭圆度 Dev 井斜 原始测井曲线代码 AMP5 第五扇区的声幅值 AMP6 第六扇区的声幅值 AMVG 平均声幅 AO10 阵列感应电阻率
Z：Elevation，海拔。所有的海拔都是从平均海平面（MSL，Mean Sea Level）算起的，向上海拔增加，向下海拔减小，MSL处为0。所以就很明显，MSL以下的Z值都是负值（图2）。 KB：Kelly Bushing，补心海拔。其实，Kelly Bushing是方钻杆补心，在Petrel中就用KB表示方钻杆补心的海拔。很明显，这个值也是以平均海平面MSL为基准的（见图1）。 MD：Measured Depth，测量深度。从方钻杆补心（KB）开始，沿着井眼测量的深度（图2）。很明显，测量，包括测井的时候，都是从方钻杆补心开始。 TVD：True Vertical Depth，垂直深度。从方钻杆补心（KB）开始，井眼的真实垂向深度（图2）！老外很严谨，为什么要加一个“True”啊，因为对应MD不“True”。直井MD 和TVD是一样的，斜井，水平井就不同了，MD更多的体现了“长度”的概念，TVD才是True 的“垂向深度”。 TVDSS：True Vertical Depth Sub-sea。平均海平面下垂直深度（图2）。这个说的很明白，Sub-sea，从海平面往下真实的垂向深度。不过，这个sea是平均海平面。 补心高（度）：补心的高度，指方钻杆补心至地面的距离，也就是方钻杆补心的地面高度。 补心海拔=地面海拔+补心高（度）。 个人认为，“测井都是从地面为基准面开始测试的”这句话是不对的。 如果一个油田单位的测井资料里没有做特别地说明其测井深度是从哪一个位置为起算点的，那么石油行业标准的通常做法是测井的基准面是该井自身的补心海拔高度，即地面海拔再加上补心高的位置，也就是钻井平台的位置（也就是KB）。 下表层套管、技术套管和油层套管时，都是从补心海拔为起算点的，实际测井时，也可能表层有一段不去测量，但测量仪器依然是从钻井平台放下去的，所以深度起算点依然是该井的补心海拔。由于钻井及下套管的深度是绝对准确的，而表层套管末端与与技术套管或油层套管有一个连接处，通常叫做表套鞋或套管脚，这个深度是准确的，所以一般是用这个深度较准测井深度，而这个深度就是从钻井平台即补心海拔为起算点的。 所以，石油行业早就规定，最后提供的测井曲线深度一般必须是从该井的补心海拔为起算点的，使用时如果需要做补心海拔校正，则必须将测井深度减去补心海拔，使得测井数据均从大地水准面为起算点。 在做井震标定时，必须充分清楚地震数据起算点与测井数据起算点是否一致，不一致的必须先做到一致，否则后面的井震标定准确性就无从谈起.虽然单纯地通过识别地震标志层来标定是业界常用的做法，但是很明显，在井震数据起算点都没有校正到统一的基准面上来的前提下，这样做显然是不可取的、有很大风险的、有时甚至会产生严重的井震不一致的错误，因此，建立地震工区时，必须清楚地震数据体的基准面是多少，表层替换速度是多少，这在地震采集及地震处理时就已经是确定知道的，而且，静校正的重要内容之一就是基准面校正，力图消除地表起伏不一致带来的时差，在确定了地震基准面和表层替换速度后，在加载测井曲线数据和钻测井分层数据时，井上的深度由于是从该井补心海拔为起算点的，所以就必须将井深度校正到地震基准面上来，这里有两种做法，针对有些软件会自动做这样的校正，所以此时应该带着补心海拔数据加载井数据，而针对有些软件不自动做这种基准面校正，则需要事先将井数据深度校正到地震基准面上来，校正的方法很简单，先用
第一讲测井曲线的识别及应用 钻井取芯、岩屑录井、地球物理测井是目前比较普及的三种认识了解地层的方法。钻井获取的岩芯资料直观、准确，但成本高、效率低。岩屑录井简便、及时，但干扰因素多，深度有误差，岩屑易失真。测井是一种间接的录井手段，它是应用地球物理方法，连续地测定岩石的物理参数，以不同的岩石存在着一定物性差别，在测井曲线上有不同的变化特征为基础，利用各种测井曲线显示的特征、变化规律来划分钻井地质剖面、认识研究储层的一种录井方法；具有经济实用、收获率高、易保存的优势，是目前我们认识地层的主要途径。 鄂尔多斯盆地常规测井系列分为综合测井和标准测井两种。 综合测井系列：重点反映目的层段钻井剖面的地层特征。测量井段由井底到直罗组底部，比例尺1：200。由感应、八侧向、四米电阻、微电极、声速、井径、自然电位、自然咖玛八种测井方法组成。探井、评价井为了提高储层物性解释精度，加测密度和补偿中子两条曲线。 标准测井系列：全面反映钻井剖面地层特征，测量井段由井底到井口（黄土层底部），比例尺1：500，多用于盆地宏观地质研究。过去标准测井系列较单一，仅有视电阻率、自然咖玛测井等两三条曲线。近几年完钻井的标准测井系列曲线较完善，只比综合测井系列少了微电极测井一项。 一、测井曲线的识别 微电极系测井、四米电阻测井、感应—八侧向测井、都是以测定岩石的电阻率为物理前提，但曲线的指向意义各异。微电极常用于判断砂岩渗透性和薄层划分。感应—八侧向测井用于判定砂岩的含油水层性能。四米电阻、声速、井径、自然电位、自然咖玛用于砂泥岩性划分。它们各有特定含义，又互相印证，互为补充，所以，我们使用时必须综合考虑。 1、微电极测井 大家知道，油井完钻后由井眼向外围依次是：泥饼、冲洗带、侵入带、地层。泥饼是泥浆中的水分进入地层后，吸附、残留在砂岩壁上的泥浆颗粒物。冲洗带是紧靠井壁附近，地层中的流体几乎被钻井液全部赶走了的部分；其深入地层的范围一般约7—8厘米。侵入带是钻井液与地层中流体的混合部分。
石油测井生产安全技术(二)测井设备及主要部位 一、测井绞车 (一)测井绞车的用途 测井、射孔等作业使用的电缆是缠放在绞车滚筒上，滚筒借助于汽车发动机的动力而转动，从而控制电缆在井内按要求的速度上提和下放。 (二）测井绞车的结构 1.汽车底盘。供给绞车动力，装载并运移绞车、电缆及其他配套设备。 2.传动系统。包括动力选择箱、液压泵、液压马达等液压动力传动设备、减速器、传动轴及传动链条。传动系统担负着作业所需动力的传递。 3.绞车。用于测井或射孔时起升或下放电缆、测井仪器及工具。 4.车身和支承底盘。用于支承绞车及其传动系统等，并提供驾驶室、操作室和绞车室。 5.操作装置。包括副排档装置、副油门装置、副离合器、盘绳器及刹车装置。用于作业时对绞车控制或操纵。 6.气路系统。用于设备的控制或操纵。 (三)测井绞车的安全操作 操作绞车就是通过操纵动力和变速系统使电缆滚筒以不同的速度和方向转动，从而使电缆及测井仪器在井中下放或上提，达到完成各项作业的目的。操作绞车只要做到操作措施得当、操作准确并做到井口慢、井底慢、特殊井段慢、遇阻、遇卡慢等，就能做到安全生产。具体说有以下操作要点： 1．测井绞车应摆放在距井25m远的上风头位置，对正井口滑轮，打好掩木。 2．起下电缆时，速度要均匀，不准猛提、猛刹，随时观察电缆运行张力读数，及时判断遇阻、遇卡。在进行井壁取心作业时，拉力增到25kN时，必须立即停车，然后慢速上下活动防止拉断岩心筒的钢丝绳，以免岩心筒落井。 3．仪器(射孔器)放人或起出井口时，应注意听从井口操作手和操作工程师的指挥，防止拉掉或摔坏仪器(射孔器)，甚至发生伤人事故。
4．注意盘齐电缆，同时做好电缆的清洁保养和防锈维护。 5．在斜井、“狗腿子”井等特殊井况下作业时，容易遇阻和遇卡，仪器和电缆下放速度要比直井慢，下放时要保持匀速，不准高速下放。发现遇阻时，不准硬冲，同时应避免仪器在井中长时间停留，要及时上提，防止遇卡。 6．井壁取心上提至套管鞋前，过油管射孔上提到油管喇叭口前，必须放慢速度，等仪器进入套管(油管)后再加快速度，防止卡掉仪器。 7．如果发现遇卡，应立即停车，如果仪器可下放，则慢速下放电缆，然后缓慢上提，这种情况比较容易解卡。若仪器能下放，则可慢速上提电缆，张力不能超过悬重15kN，若仍不能解卡，应采取其他解卡措施。 8．每次作业，应认真填写绞车运行记录。 二、测井电缆 （一）测井电缆的作用 1.输送各种下井仪器和工具。 2.向井下仪器供电和传送控制信号。 3.将井下仪器采集到的信号传送到地面。 （二）测井电缆的性能要求 1.有较强的抗拉强度。 2.有良好的韧性能盘绕在绞车滚筒上。 3.必须具有多股缆芯，导电、绝缘、抗干扰等性能良好，以满足供电和传送电信号的要求。 4.缆芯的绝缘材料必须具有耐高温、耐高压、防腐蚀、防脆化的良好性能。
编号成果类型所属专业（课程）主题发表刊物刊期教案计算机CAD软件应用 精确绘图方法
理论授课教案（首页） 授课日期 班级 课题：精确绘图方法 教学目的要求： 1、掌握如何通过常用的指定坐标法绘制图形，在不输入坐标的情况下快速、精确地绘制图形。 2、掌握各种精确绘图手段，并能够利用它们进行精确绘图。 教学重点、难点： 1、了解坐标系的创建、表示方法，掌握如何设置和使用极轴追踪和极轴捕捉，从而快速绘制倾斜直线。 2、了解对象捕捉和对象捕捉追踪的使用方法。授课方法：任务驱动法、实例演示法、分组讨论法、练习法 教学参考及教具（含多媒体教学设备）： 1、软件教程书籍、软件练习习题图册 2、软件教学课件及示范视频、多媒体电教设备 授课执行情况及分析： 通过讲解和上机操作练习，学生基本能够掌握精确绘图方法。巡回指导过程中，发现同学们的绘图方法和步骤不尽相同， 基本有自己的绘图思路，整体效果良好，课堂气氛活跃，大部分 同学们可以独立完成课堂作业。 板书设计或授课提纲 §4精确绘图方法 AutoCAD的坐标系设置绘图单位与绘图边界 使用栅格、捕捉和正交辅助定位 通过捕捉图形几何点精确定位点 使用追踪精确定位点动态输入与选择预览
教学环节 及时间分配 教学内容教学方法 组织教学 导入新课 5分1.问候，检查学生出勤情况 2.检查学生课前准备情况 与手工画图相比，使用AutoCAD画图的 优势之一在于它为用户提供了众多辅助画 图手段。例如，利用坐标可轻松定位点， 利用捕捉和栅格可控制光标的精确移动， 利用正交和极轴追踪可绘制水平、垂直或 倾斜直线等。这一章我们就来细细讲述各 种精确绘图方法。 §4精确绘图方法 一．AutoCAD 的坐标系 默认情况下，AutoCAD自动将当前坐标 系设置为世界坐标系，即WCS，它包括X轴 和Y轴。如果在三维空间工作，还有一个Z 轴。在用户界面的底部状态栏左侧显示了 当前光标所在位置的三维坐标值，同时在 图形编辑窗口显示了坐标系图标，以指明 当前坐标系的性质。 考勤 登记 提问互 动，引导 学生进入 学习角色
常用测井曲线名称 测井符号英文名称中文名称 Rt true formation resistivity. 地层真电阻率 Rxo flushed zone formation resistivity 冲洗带地层电阻率 Ild deep investigate induction log深探测感应测井 Ilm medium investigate induction log中探测感应测井 Ils shallow investigate induction log 浅探测感应测井 Rd deep investigate double lateral resistivity log 深双侧向电阻率测井Rs shallow investigate double 浅双侧向电阻率测井 lateral resistivity log RMLL micro lateral resistivity log 微侧向电阻率测井 CON induction log 感应测井 AC acoustic 声波时差 DEN density 密度 CN neutron 中子 GR natural gamma ray 自然伽马 SP spontaneous potential 自然电位 CAL borehole diameter 井径 K potassium 钾 TH thorium 钍 U uranium 铀 KTH gamma ray without uranium 无铀伽马 NGR neutron gamma ray 中子伽马 常用测井曲线名称
A1R1 T1R1声波幅度 A1R2 T1R2声波幅度 A2R1 T2R1声波幅度 A2R2 T2R2声波幅度 AAC 声波附加值AAVG 第一扇区平均值 AC 声波时差 AF10 阵列感应电阻率AF20 阵列感应电阻率AF30 阵列感应电阻率AF60 阵列感应电阻率AF90 阵列感应电阻率AFRT 阵列感应电阻率AFRX 阵列感应电阻率 AIMP 声阻抗 AIPD 密度孔隙度 AIPN 中子孔隙度 AMAV 声幅 AMAX 最大声幅 AMIN 最小声幅AMP1 第一扇区的声幅值AMP2 第二扇区的声幅值AMP3 第三扇区的声幅值AMP4 第四扇区的声幅值AMP5 第五扇区的声幅值AMP6 第六扇区的声幅值AMVG 平均声幅 AO10 阵列感应电阻率AO20 阵列感应电阻率AO30 阵列感应电阻率AO60 阵列感应电阻率AO90 阵列感应电阻率 AOFF 截止值AORT 阵列感应电阻率AORX 阵列感应电阻率APLC 补偿中子 AR10 方位电阻率 AR11 方位电阻率 AR12 方位电阻率 ARO1 方位电阻率 ARO2 方位电阻率 ARO3 方位电阻率 ARO4 方位电阻率 ARO5 方位电阻率
测井曲线代码一览表 常用测井曲线名称 测井符号英文名称中文名称 Rt true formation resistivity. 地层真电阻率Rxo flushed zone formation resistivity 冲洗带地层电阻率Ild deep investigate induction log 深探测感应测井Ilm medium investigate induction log 中探测感应测井Ils shallow investigate induction log 浅探测感应测井Rd deep investigate double lateral resistivity log 深双侧向电阻率测井Rs shallow investigate double lateral resistivity log 浅双侧向电阻率测井RMLL micro lateral resistivity log 微侧向电阻率测井CON induction log 感应测井 AC acoustic 声波时差DEN density 密度 CN neutron 中子 GR natural gamma ray 自然伽马 SP spontaneous potential 自然电位 CAL borehole diameter 井径 K potassium 钾 TH thorium 钍 U uranium 铀 KTH gamma ray without uranium 无铀伽马 NGR neutron gamma ray 中子伽马 5700系列的测井项目及曲线名称 Star Imager 微电阻率扫描成像 CBIL 井周声波成像 MAC 多极阵列声波成像 MRIL 核磁共振成像 TBRT 薄层电阻率 DAC 阵列声波 DVRT 数字垂直测井 HDIP 六臂倾角 MPHI 核磁共振有效孔隙度 MBVM 可动流体体积 MBVI 束缚流体体积 MPERM 核磁共振渗透率 Echoes 标准回波数据 T2 Dist T2分布数据 TPOR 总孔隙度 BHTA 声波幅度 BHTT 声波返回时间 Image DIP 图像的倾角 COMP AMP 纵波幅度 Shear AMP 横波幅度 COMP ATTN 纵波衰减 Shear ATTN 横波衰减 RADOUTR 井眼的椭圆度 Dev 井斜
［收稿日期］２００６１０２０ ［作者简介］吴爱平（１９７７），男，２００２年大学毕业，硕士，讲师，现主要从事信号采集与处理方面的教学与研究。基于ＦＰＧＡ的测井深度系统设计 吴爱平，付青青 （长江大学电子信息学院，湖北荆州４３４０２３） ［摘要］针对单片机处理深度信号存在速度慢、外围电路复杂等问题，提出了用现场可编程门阵列实现深 度系统设计的方法，给出了设计的整体硬件框图和ＦＰＧＡ设计的顶层图，并对设计进行了仿真。试验结 果表明，使用该方法设计的深度处理系统结构简单，速度快，可靠性高。 ［关键词］ＦＰＧＡ；测井深度；设计 ［中图分类号］Ｐ６３１畅８８ ［文献标识码］Ａ ［文章编号］１６７３１４０９（２００６）０４０１０２０３ 深度信息是一个非常重要的测井信息，如果深度测量不准，可能会导致测井资料作废，甚至有带来生产事故的危险，所以它在测井过程中起着相当重要的作用。目前进行深度系统设计主要是采用单片 机，用单片机设计的深度系统［１］，外围电路复杂，执行速度慢，但开发风险小。为了克服采用单片机设 计的深度系统的弊端，笔者提出了采用ＦＰＧＡ（现场可编程门阵列）设计深度系统的方法，采用该方法设计的深度系统，硬件电路简单，执行速度快，可靠性高，升级改造容易。 图１ 测井深度系统总体硬件框图１ 硬件设计 系统硬件框图如图１所示，深度信号送至ＦＰＧＡ 处理，同时ＦＰＧＡ通过ＰＣＩ９０５２接受工控机发出的深 度初始设置命令，并回送计算的深度脉冲数给工控 机。ＦＰＧＡ中设计了地址译码、读写控制、数据选 择、方向判断、深度计数、深度中断发生器、数据锁 存７个模块，其中用到了２３个Ｉ／Ｏ（输入／输出端 口）。ＦＰＧＡ选择Ａｌｔｅｒａ公司的ＥＰ１Ｃ３ＴＣ１４４［２］。由于Ｃｙｃｌｏｎｅ系列器件是基于ＳＲＡＭ工艺，掉电后信息不能保存，每次上电时要对芯片进行配置后才可使用， 设计中选用的配置芯片是Ａｌｔｅｒａ公司的ＥＰＣＳ１。通过ＪＴＡＧ（联合测试行动小组）进行在线编程，ＡＳ（主动串行）进行在线配置，配置电路完成ＦＰＧＡ程序的上电配置。 配置寄存器选用９３Ｃ４６，完成ＰＣＩ９０５２上电时的自动配置，确保接口电路按预先设置的方式进行工作。９３Ｃ４６存储了ＰＣＩ９０５２重要的配置信息，如设备号ＤＩＤ、制造商号ＶＩＤ、子设备号ＳＤＩＤ、子制造商号ＳＶＩＤ、中断号、设备类型号、局部空间基地址、局部空间描述符、片选响应以及局部响应控制ＣＮＴＲＬ等信号。配置寄存器的内容非常重要，它直接关系到整个系统能否正常工作，如果配置不正确，即使硬件设计没有一点错误，电路也很难正确地工作［３］。 ２ 软件设计 软件设计选用Ａｌｔｅｒａ公司的ＱｕａｔｒｕｓⅡ５畅０［４］ 软件开发工具。深度信号的处理由ＦＰＧＡ实现，深度信号Ａ和Ｂ从ＦＰＧＡ的全局时钟进入。在ＦＰＧＡ内部做了接口模块、测井方向判决模块、一个３２位的可逆计数器和深度中断发生器（图２）。当进行上提测井时，Ａ信号相位超前Ｂ信号，方向模块输?２０１?长江大学学报（自科版） ２００６年１２月第３卷第４期理工卷ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＵｎｉｖｅｒｓｉｔｙ（ＮａｔＳｃｉＥｄｉｔ） Ｄｅｃ畅２００６，Ｖｏｌ畅３Ｎｏ畅４Ｓｃｉ＆ＥｎｇＶ