Paleoscan断层数据输出格式和加载入其他软件的格式整理方式2018.3.pptx

1、,1,2,3,4,PaleoScan软件可以输出5种断层格式，但实际上前两个基本内容差不多，所以我把其分为四类,如何把PaleoScan输出的断层变为landmark可以加载的断层格式，在最后面介绍,第一类格式，是把PaleoScan软件中内插的断面按照你选的地震测网在inline方向上每5道一抽稀，同时保存的文件名与在软件中是一致的。这种格式输出的断层棍 id 编号入下图所示，一个剖面上的断点的id是一样的，这和landmark 1 3或1223这样的编号方式不一样。这一类格式最大的问题是：比如F1断层输出后断层棍在inline1、6、11、16，F2断层的断层棍可能会在inline 3、8

2、、13、18，F3断层的断层棍可能会在inline 4、9、14、19上,1,2,3,4,id,第二类格式，是把PaleoScan软件中我们解释的断层棍输出，但是同时保存的文件名是按照序号重新命名，与软件中的名字无关联。这种格式输出的断层棍 id 编号入下图所示，但是需要注意的是其输出是将一条比如有3个解释节点的断层，变成有分开的两段组成的如图2所示后面我们会介绍如何在这种格式基础上编辑为landmark可以识别的格式如果需要paleoscan解释的断层可以在其他软件可以用，那就要和常规软件一样按照一定的统一间隔解释断层,id,示意,实际中是重叠显示，这里只是示意,第三类格式，是PaleoSc

3、an中内插的断面按照所选的测网按照xyz格式1X1的输出，但是没有id号，文件名也是重新命名，与软件中不一致如果按照这种格式输出断层，就不存在几X几的问题，可以较好的发挥paleoscan的优势，我们可以把输出结果当层位加载入其他软件中。我公司目前准备是在这种输出格式基础上，写一个小程序进行抽稀，并变为landmark、geoeast等软件的格式,第四种格式 是gocad软件格式，是建模软件用的，petrel可以加载，landmark等好像不行,后面将提供1、将第二种整理为geoeast、landmark格式为例2、将第一类格式整理为landmark格式为例3、将第三类格式整理为landmar

4、k格式为例,需要有以下知识：理解landmark等断层格式内容，会用ultraedit软件的列对齐等功能，最好理解excel的数据计算功能推荐看用第三类进行格式转换！,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,增加一列并填上断层名，保存为txt文件,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,所有断层都这么处理,所有断层都保存在一个文件夹内，并新建一个文件000.txt,打开键入，保存,改后缀为1.bat,双击执行,所有文件合并为fff.dat,用excel打开,删掉copy *.* fff.dat,断层格式整理，将第二类格式整

5、理为geoeast、landmark格式为例，,这一步是要合并所有单条断层到一个文件当中,依次手动键入,=IF(AND(B2=B3,C2=C3,D2=D3),0,1),=IF(AND(F2=F1,E2=E1),3,0),=F2+G2,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,加序列号，并把全部数据按照数据重新粘贴，不保留公式,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,排序,删除该列为0的行,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,排序,断层格式整理，将第二类格式整理为geoeast、land

6、mark格式为例，,删除没用的,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,删掉,4代表一条断层棍在一个剖面上的结束,手动键入第一个剖面的id号,=IF(F2=4,IF(B3=B2,G2+1,1),G2),断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,这里的F2、b3、b2、g2等参数要根据实际数据调整,完成111格式的初步整理，后续按照软件输出的断层格式，用UltraEdit软件，把tab替换为空格后，按照参考格式去对齐，即可加载,另存为文本文件，注意删除尾行,断层格式整理，将第二类格式整理为geoeast、landmark格式为

7、例，,删掉,手动键入，第一个剖面的id号,=IF(F2=1,IF(F1=4,1,2),3),123格式的初步整理即完成,然后再想办法按照landmark断层格式用UltraEdit对齐数据，增加一些必要的内容即可加载,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,删除,断层格式整理，将第一类格式整理为landmark格式为例，,第一类格式，先把所有断层合并为一个文件，方法前面有,=IF(E2=E1,IF(E2E3,3,2),1),手动键入1,最后一行手动改为3,断层格式整理，将第一类格式整理为landmark格式为例，,剩下的就是用ultraedit按照landm

8、ark对其格式，增加内容即可,断层格式整理，将第三类格式整理为geoeast、landmark格式，以inline方向为例,大致步骤1，根据测网数据建立xy和线道的数学关系2，将xyz格式断面数据转化为线道z格式的断面数据（保留原始xy）3，根据想抽稀的测网对断面进行抽稀4，抽稀后数据加id号5，根据其他软件断层格式进行列对其并添加相应内容,增加一列并填上断层名，保存为txt文件,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,所有断层都这么处理,所有断层都保存在一个文件夹内，并新建一个文件000.txt,打开键入，保存,改后缀为1.bat,双击执行,所有文件合并为

9、fff.dat,用excel打开,删掉copy *.* fff.dat,断层格式整理，将第二类格式整理为geoeast、landmark格式为例，,这一步是要合并所有单条断层到一个文件当中,读取地震测网4点坐标，以及工区内至少任意一点的线道xy（推荐有多个），以及线道与xy坐标之间的夹角，尽量用道头读取软件在sgy道头里面读取，不是在平面上图上去图取新打开一个excel表填写进去,=RADIANS(C10)把度换算为弧度,=C2*COS(1.54219786124659)-D2*SIN(1.54219786124659),=C2*SIN(1.54219786124659)+D2*COS(1.5

10、4219786124659),坐标转换,这个值是转换后的弧度值,用线和X1做交会图，二者一条直线说明是正确的，横坐标是X1，纵坐标是线，然后进行回归，得到公式，如果不对可以试着用线和Y1做交汇分别得到线道号和x1y1的数学关系,公式的小数点后位数一定要够，否则误差很大,公式外加一个四舍五入函数即可,=C1*COS(1.54219786124659)-D1*SIN(1.54219786124659),=C1*SIN(1.54219786124659)+D1*COS(1.54219786124659),x,x1,y1,y,线,道,将前面整理的公式带入断层文件当中,x,x1,y1,y,线,道,Z,只

11、保留数据清除公式,假设，从inline100开始，inline每4道一个断层棍如果数据点太密加载后断层棍显示会呈阶梯状，所有xline方向也要抽稀，该间隔普遍较大，这里以16为例从道号10开始，,=IF(INT(G1-100)/4)=(G1-100)/4),G1,0),=IF(INT(H1-10)/16)=(H1-10)/16),H1,0),只保留值数值，然后把计算结果是0的行全部删掉,x,x1,y1,y,线,道,Z,=IF(AND(G2=G1,H2=H1)0,1),手动键入1,先仅保存数值，再把这一列是0的数据行也删掉,线,道,删除完后，先按线号排序，再按照断层名排序,x,x1,y1,y,线,道,Z,保留线号和道号后面可能有用,我这没有保留线道号,手动输入1,=IF(G2=G1,IF(G2G1,H1+1,H1),1),在这个基础上,id,这是geoest格式断层id,线,线,整理完成后，可以加载入paleoscan，三维显示检查下再按照对应软件格式进行增加内容对其，如果加载后有问题可以尝试把线道号重新算为xy坐标再去加载,