CN105701848A - 一种地层界线图层的自动化生成方法 - Google Patents
一种地层界线图层的自动化生成方法 Download PDFInfo
- Publication number
- CN105701848A CN105701848A CN201610023762.2A CN201610023762A CN105701848A CN 105701848 A CN105701848 A CN 105701848A CN 201610023762 A CN201610023762 A CN 201610023762A CN 105701848 A CN105701848 A CN 105701848A
- Authority
- CN
- China
- Prior art keywords
- key element
- layer
- boundary
- stratum
- surface key
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000008021 deposition Effects 0.000 claims abstract description 12
- 239000010410 layer Substances 0.000 claims description 89
- 125000004122 cyclic group Chemical group 0.000 claims description 37
- 239000011435 rock Substances 0.000 claims description 26
- 101000685990 Homo sapiens Specifically androgen-regulated gene protein Proteins 0.000 claims description 13
- 102100023355 Specifically androgen-regulated gene protein Human genes 0.000 claims description 13
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 9
- 238000005755 formation reaction Methods 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000003325 tomography Methods 0.000 claims description 8
- 239000011229 interlayer Substances 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 4
- 230000004069 differentiation Effects 0.000 claims description 2
- 238000000151 deposition Methods 0.000 description 10
- 239000004927 clay Substances 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000019771 cognition Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000000213 tara gum Substances 0.000 description 2
- SFXVRNLGOSAWIV-UHFFFAOYSA-N 4-amino-N-[4-(methoxymethyl)phenyl]-7-(1-methylcyclopropyl)-6-(3-morpholin-4-ylprop-1-ynyl)pyrrolo[2,3-d]pyrimidine-5-carboxamide Chemical class NC=1C2=C(N=CN=1)N(C(=C2C(=O)NC1=CC=C(C=C1)COC)C#CCN1CCOCC1)C1(CC1)C SFXVRNLGOSAWIV-UHFFFAOYSA-N 0.000 description 1
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 1
- 239000001747 Potassium fumarate Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 229910052700 potassium Inorganic materials 0.000 description 1
- 239000011591 potassium Substances 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种地层界线图层的自动化生成方法。该方法包括如下步骤:(1)基于沉积、火山地层界线图层和变质地层单位图层,搜寻出图层中所有邻接地层,构建地层邻接关系图;(2)基于地层邻接关系图提取相邻地层的公共边界,生成地层界线;(3)基于断层图层,识别断层接触关系;(4)基于产状图层,识别基本接触关系,完成地层界线图层的自动化生成。与人工制作的方法相较,本发明方法生成的地层界线图层不仅高效、快捷,且通过降低人工参与程度而尽可能避免人为误差的引入,满足了地质图空间数据库地层界线图层制作的需要。
Description
技术领域
本发明属于地理信息技术与应用领域,具体涉及一种基于沉积和火山地层单位图层、变质地层单位图层、产状图层和断层图层,自动化生成地层界线图层的方法。
背景技术
在地质图空间数据库建设过程中,地层界线图层的制作是一项必不可少的工作。在中国地质调查局所制定的行业标准“数字化地质图图层及属性文件格式[DZ/T0197-1997]”中,明确地对地层界线图层的属性数据做出了规定。
目前,地层界线图层的制作主要通过人工制作来实现。制作方法主要包括两方面工作:一是扫描并矢量化已有的地质图件得到地层界线,并对界线图元进行编辑,完成地层界线的提取;二是依据不同地层接触关系在地质图上的表现,人工识别接触关系并进行属性信息的录入。该制作方法不仅存在工作量大、人工识别效率低,而且易引入人为误差,如接触关系识别错误、属性信息错误、与相关地层单位图层界线不一致等缺点。
在地质图空间数据库建设过程中,着重制作的沉积和火山地层单位图层、变质地层单位图层中,已蕴含了相关地层界线位置和接触关系信息。其中,地层界线的位置在沉积和火山地层单位图层和变质地层单位图层中表现为具有相邻拓扑关系的两地层要素的公共边界。接触关系可以基于地层先后顺序、产状数据、断层数据等进行自动化识别。
发明内容
为此,本发明主要针对人工制作地层界线图层时出现的工作量大、识别效率低,以及易引入人为误差等问题,提出了一种基于沉积和火山地层单位图层、变质地层单位图层、产状图层和断层图层,进行自动化生成地层界线图层的方法。
本发明的技术方案是:
一种地层界线图层的自动化生成方法,包括如下步骤:(1)基于沉积、火山地层界线图层和变质地层单位图层,搜寻出图层中所有邻接地层,构建地层邻接关系图(StratumAdjacencyRelationGraph,简称SARG);(2)基于地层邻接关系图提取相邻地层的公共边界,生成地层界线;(3)基于断层图层,识别断层接触关系;(4)基于产状图层,识别基本接触关系,完成地层界线图层的自动化生成。
所述步骤(1)的具体步骤如下:
步骤11:加载沉积和火山地层单位图层、变质地层单位图层数据,并将相应的地层面要素读入集合S={si|i=1,2,…,n},n为地层面要素的个数;
步骤12:读出每个地层面要素si的顶点集vi={Pj|j=1,2,…,m},其中Pj表示构成地层面要素si的顶点,m为顶点的个数;并对每个地层面要素si的顶点集vi进行如下的预处理:
d)将顶点集vi中的点Pj按顺时针排序;
e)判断每个顶点集vi中是否存在P1=Pm,即首节点等于尾节点,若存在,则跳转下一步
c),对顶点集vi进行处理;否则,跳转步骤13;
f)删除顶点集vi的尾节点Pm,保证首尾不相接,即vi={Pj|j=1,2,…,m-1};
步骤13:对图层中地层面要素的邻接关系进行初判:若地层面要素sa的顶点集va的最小外接矩形Ra和要素sb的顶点集vb的最小外接矩形Rb有交点,则跳转步骤14;否则,跳转步骤16;其中,a为循环变量,a∈[1,n-1],初值为1;b为循环变量,b∈[a+1,n],初值为a+1;
步骤14:对两地层面要素sa和sb的邻接关系进行细判,具体过程如下:
d)计算顶点集va和vb的公共点集T={p1,p2,...,pk},k为公共点的个数,并记录T中点在顶点集va中序号,记为N={n1,n2,...,nk};
e)若公共点集T中点的个数k≥2,则跳转下一步c);否则,跳转步骤16;
f)判断公共点序号集N内是否存在连续的点号nt,nt+1满足∣nt-nt+1∣=1,若存在,即地层面要素s1和s2存在邻接关系,则跳转步骤15;否则,跳转步骤16;
步骤15:构建地层邻接关系图的边d:计算顶点集va与vb的重心点Ca和Cb,点Ca和Cb作为地层邻接关系图的顶点,连接点Ca和Cb作为地层邻接关系图的边d;
步骤16:若循环变量a>n-1,则输出SARG={d1,d2,...dw},w为地层邻接关系图SARG的边的总个数;否则,循环变量b加1,此时若b≤n,则直接跳转步骤13;否则循环变量a加1,循环变量b=a+1,再跳转步骤13。
所述步骤(2)的具体步骤如下:
步骤21:读取SARG={d1,d2,...dw}中第i条边di,i为循环变量,初值为1;
步骤22:求出地层面要素sa和sb的顶点集va={p1,p2,...,pm}和vb={q1,q2,...,qn},m、n分别为地层面要素sa、sb的顶点个数,并对顶点集va和vb进行如下的预处理:
d)将顶点集va={p1,p2,...,pm}和vb={q1,q2,...,qn}中的点均按顺时针排序;
e)判断顶点集va和vb是否同时满足p1=pm与q1=qn,即首节点等于尾节点;若满足,则跳转步骤23;否则,对不符合条件的顶点集va或vb,执行下一步c);
f)对顶点集va或vb继续添加第一个点p1或q1,保证首尾相接;
步骤23:计算顶点集va和vb的公共点集T={t1,t2,...,tj},j为公共点的个数,并记录公共点集T中各点分别在顶点集va和vb中的序号,分别记为Qa和Qb;
步骤24:公共点集T的分段处理,分段处理步骤如下:
d)设k为循环变量,k∈[1,j-1],初值为1;读取序号Qa中第k个点号i'、第k+1个点号j';同样,读取序号Qb中第k个点号ii'、第k+1个点号jj',若同时满足∣i'-j'∣=1且∣ii'-jj'∣=1,则点号i'和点号j'、以及点号ii'和点号jj'所对应的两点属于同一段公共边界;
e)循环变量k加1,若k<j,继续执行上一步a)直至公共点集T分段结束,并依次将属于同一段边界的点构成边e,得到公共边界Ea/b={e1,e2,..,em'},m'表示公共边界的段数;如果m'>2,执行下一步c);否则跳转步骤25;
f)依次两两比较构成边et的最后点坐标tL(x,y)与构成边et+1的第一个点坐标tF(x,y),t∈[1,m'],若存在tL(x,y)=tF(x,y),则将边et与et+1合并,边et与et+1属于同一段公共边界;
步骤25:循环变量i加1,若i≤w,循环执行步骤21-24,得到整个图层的地层界线集合E={Ea/b,Ec/d,...},其中Ea/b,Ec/d分别代表邻接地层sa与sb、邻接地层sc与sd间的地层界线。
所述步骤(3)的具体步骤如下:
步骤31:加载断层图层数据,并以断层图层中的断层线F={f1,f2,...,fm}为中心轴线,以距离参数d生成断层缓冲区FB;
步骤32:提取与断层数据中断层线F={f1,f2,...,fm}相交的地层邻接关系图的边集合D={d1,d2,...,dn},n为地层邻接关系图与断层线F相交的条数;
步骤33:搜寻边集合D中边di所对应的地层面要素Sa与地层面要素Sb间的地层界线Ea/b={e1,e2,...,em},其中i为循环变量,i∈[1,n],初值为1;
步骤34:判断断层接触:若Ea/b={e1,e2,...,em}的某段界线ej在断层缓冲区FB内,则界线ej两侧地层间的接触关系即为断层接触关系;
步骤35:循环变量i加1,若i≤n,则重复执行步骤33、34,直至图层中所有的断层接触提取完毕。
所述步骤(4)的具体步骤如下:
步骤41:基于地层岩性信息,区分侵入体与非侵入体;将每个地层的区分结果记录到沉积和火山地层单位图层、变质地层单位图层的地层面要素的“岩类”属性中;
步骤42:基于岩石地层单位简表,生成序列表X,且较老的地层序号数值较小;
步骤43:读取地层界线集合E中的第i条界线要素Ea/b所对应的地层面要素sa和sb的岩石地层单位Rocknamea和Rocknameb,其中i为循环变量,i∈[1,k],初值为1,k为地层界线的条数;
步骤44:读取地层面要素sa和sb的“岩类”属性;若地层面要素sa和sb的“岩类”均为非侵入体,则跳转步骤45;否则,跳转步骤47;
步骤45:整合接触和不整合接触的判断;读取岩石地层单位Rocknamea和Rocknameb在序列表X上的序号xa、xb,若序号xa、xb相邻或相同,即xa=xb或者∣xa-xb∣=1,则地层面要素sa和sb的接触关系为整合接触;否则,为不整合接触,跳转步骤46,进行下一步判断;
步骤46:平行不整合接触和角度不整合接触判断,具体判断方法如下:
c)加载产状图层数据,并读取地层面要素sa、sb的产状要素:走向Strikea和Strikeb、倾向Dipa和Dipb、以及倾角Anglea和Angleb;若任一产状要素信息无法获取,则地层面要素sa和sb的接触关系仅判断为不整合接触,跳转步骤48;
d)判断Strikea=Strikeb且Dipa=Dipb且∣Anglea-Angleb∣≤δ是否成立,其中0°<δ<10°;若成立,则地层面要素sa和sb的接触关系为平行不整合接触;否则,为角度不整合接触;
步骤47:沉积接触与侵入接触的判断;读取岩石地层单位Rocknamea和Rocknameb在序列表X上的序号xa、xb;当地层面要素sa的“岩类”属性为侵入体时,若xa≥xb,则为侵入接触;否则,则为沉积接触;同样,当地层面要素sb的“岩类”属性为侵入体时,且xa<xb,则为侵入接触;否则,则为沉积接触;
步骤48:记录地层线要素的图元编号和接触关系属性;并循环步骤43-48,直至地层界线集合E中所有界线要素两侧地层间的接触关系均识别结束,完成地层界线图层的制作。
本发明的有益效果:与人工制作的方法相较,本发明方法生成的地层界线图层不仅高效、快捷,且通过降低人工参与程度而尽可能避免人为误差的引入,满足了地质图空间数据库地层界线图层制作的需要。此外,本发明方法又可作为检验地层界线图层质量的有效工具。
附图说明
图1本发明地层界线图层制作流程图;
图2庐山实验数据;
图3庐山的地层邻接关系图;
图4提取的庐山地层界线;
图5断层缓冲区;
图6与断层线相交的地层邻接关系图;
图7断层接触提取结果;
图8庐山地层界线专题图。
具体实施方式
下面结合附图并通过描述一个地层界线图层自动化生成的实例,来进一步说明本发明的效果。本实例选择1:50000的庐山沉积和火山地层单位图层、变质地层单位图层、断层图层和产状图层为实验数据,如图2所示。
具体实施过程如下:
(一)构建SARG
步骤11:加载庐山沉积和火山地层单位图层、变质地层单位图层数据,如图2所示。并将相应的地层面要素读入集合S={si|i=1,2,…,419},本实例中共包含419个地层面要素;
步骤12:读出每个地层面要素si的顶点集vi,并对每个地层面要素si的顶点集进行预处理,以地层面要素s1(大冶组)为例:
a)将顶点集v1中的点按顺时针排序,得到v1={(379043.6582,3277656.1655),(379018.7804,3277657.9852),...,(379043.6582,3277656.1655)},其中p1=(379043.6582,3277656.1655),p18=(379043.6582,3277656.1655),共包含18个顶点;
b)判断顶点集v1中是否存P1=P18。v1中P1=(379043.6582,3277656.1655),P18=(379043.6582,3277656.1655),存在P1=P18(即首节点等于尾节点),则跳转下一步c)对顶点集v1进行处理;
c)删除顶点集v1的尾节点P18,保证首尾不相接。即v1={(379043.6582,3277656.1655),(379018.7804,3277657.9852),...,(379045.2562,3277678.0115)}。
步骤13:对图层中两地层面要素sa(a为循环变量,a∈[1,418],初值为1)和sb(b为循环变量,b∈[2,419],初值为a+1)的邻接关系进行初判,以s1(大冶组)和s2(新港粘土下段)为例进行步骤13-16,即a=1,b=2时:地层面要素s1的顶点集v1的最小外接矩形R1和s2的顶点集v2的最小外接矩形R2有交点,跳转步骤14;
步骤14:对两地层面要素sa和sb的邻接关系进行细判,同样,以s1与s2为例,具体过程如下:
a)计算顶点集v1和顶点集v2的公共点集T={(379043.6582,3277656.1655),(379018.7804,3277657.9852),(378995.2579,3277644.7475)...(379043.6582,3277656.1655)},公共点的个数k=9。并记录T中点在顶点集v1中序号,记为N={1,2,3,...18};
b)公共点集T中点的个数k≥2,则跳转下一步c);
c)公共点集序号N内存在连续的点号1、2,满足∣1-2∣=1,即地层面要素s1和s2存在邻接关系,跳转步骤15。
步骤15:构建地层邻接关系图的边:仍以s1和s2为例,计算v1与v2的重心点C1(378978.798955993,3277704.63529242)和C2(379006.309367991,3277476.48297586),点C1和C2作为地层邻接关系图的顶点,连接C1和C2作为地层邻接关系图的第一条边d1;
步骤16:循环变量a=1<418,则循环变量b加1。b=3≤419,跳转步骤13,判断地层面要素s1和s3是否存在邻接关系。循环执行步骤13-15,直至地层面要素s418和s419的邻接关系判断结束,得到庐山的SARG,输出SARG={d1,d2,...d785},即存在785个地层邻接关系,如图3所示。
(二)基于SARG提取相邻接地层要素间的公共边界,即地层界线的提取
步骤21:读取SARG={d1,d2,...d785}中第i条边di,i为循环变量,初值为1。以d1为例,i=1,d1由地层面要素s1的重心点C1和地层面要素s2的重心点C2连接得到。提取地层面要素s1和地层面要素s2的公共边界,执行步骤22-26;
步骤22:求出地层面要素s1和s2的顶点集v1={(379043.6582,3277656.1655),(379018.7804,3277657.9852),..,(379043.6582,3277656.1655)}共18个顶点和v2={(379043.6582,3277656.1655),(379048.2663,3277622.0509),..,(379043.6582,3277656.1655)}共23个顶点。并对顶点集v1和v2进行如下的预处理:
a)将顶点集v1和v2中的点均按顺时针排序,步骤22中得到的v1和v2已按顺时针方向排序;
b)顶点集v1和v2均满足p1(379043.6582,3277656.1655)=p18(379043.6582,3277656.1655)与q1(379043.6582,3277656.1655)=q23(379043.6582,3277656.1655),即首节点等于尾节点,跳转步骤23。
步骤23:计算顶点集v1和v2的公共点集T={(379043.6582,3277656.1655),(379018.7804,3277657.9852),(378995.2579,3277644.7475)...(379043.6582,3277656.1655)},公共点的个数为9,即j=9。并记录T中各点在顶点集v1和v2中的序号,Q1={1,2,3,4,5,6,7,8,18}和Q2={22,21,20,19,18,17,16,15,0};
步骤24:公共点集T的分段处理,分段处理步骤如下:
a)设k为循环变量,k∈[1,j-1],初值为1。以k=1为例进行分段:循环读取序号Q1中第1个点号1、第2个点号2,同样,读取序号Q2中第1个点号22、第2个点号21,∣1-2∣=1且∣22-21∣=1,则序号Q1中点号1和点号2或序号Q2中点号22和点号21,所对应的两点属于同一段公共边界;
b)循环变量k加1,继续执行k=2、3...、9,进行分段。分段结果,表明在v1中点号1和2、2和3、3和4、4和5、5和6、6和7、7和8所对应的两点均属于同一段公共边界。因而,v1中点号1,2,3,4,5,6,7,8所对应的点,均属于同一段边界。且v1和v2仅有一段公共边界,公共边界E1/2={e1}。即m'=1,不执行c),跳转步骤25。
步骤25:循环变量i加1,i=2,i≤785,继续执行i=2、3...、785,提取地层界线,得到整个图层的地层界线集合E={E1/2,E1/58,...E417/418},共785条(926段)地层界线,并以线图层的方式展示,如图4所示。由于存在地层界线包含多段公共边界的情形,如E16/388={e1,e2},包含2段公共边界,因而地层界线E共包含926段地层界线。
(三)基于断层数据,识别断层接触关系
步骤31:加载庐山断层图层数据(图2)。并以断层图层中的断层线F={f1,f2,...,f51}为中心轴线,距离参数100m(本实例中距离参数d=100m)生成断层缓冲区FB,如图5所示;
步骤32:提取与断层数据中断层线F={f1,f2,...,f51}相交的地层邻接关系图的边集合D={d1,d2,...,d142},本实例中SARG中共有142条边与断层线F相交,如图6所示;
步骤33:搜寻边集合D中边di(i为循环变量,i∈[1,142],初值为1),所对应的地层面要素Sa与地层面要素Sb间的地层界线Ea/b={e1,e2,...,em}。以边d129为例,提取所对应的地层面要素S49与地层面要素S343间的地层界线E49/343=e1;
步骤34:判断断层接触:若Ea/b={e1,e2,...,em}的某段界线ej在断层缓冲区FB内,则ej两侧地层间的接触关系即为断层接触关系。继续以地层界线e1为例,e1在断层缓冲区FB内,则e1两侧地层即茅山组与五通组在e1处的接触关系为断层接触关系;
步骤35:循环变量i加1,若i≤142,则重复执行步骤33、34,直至图层中所有的断层接触提取完毕。断层接触提取的最终结果如图7所示。
(四)基于产状数据,识别基本接触关系,完成地层界线图层的制作。
步骤41:基于地层岩性信息,区分侵入体与非侵入体。本实例中区分结果为s7(脉岩:花岗伟晶岩脉)、s9(脉岩:石英脉)、s393(观音桥片麻岩套:高家岭钾长片麻岩)、s394(星子岩群:归宗寺变粒岩)为侵入体,其余地层面要素均为非侵入体。并将区分结果记录到沉积和火山地层单位图层、变质地层单位图层的地层面要素的“岩类”属性中;
步骤42:基于岩石地层单位简表,生成序列表X,且较老的地层序号数值较小。本实例中基于庐山地层单位简表,生成序列表X,其中最小的是Pt1xq(星子岩群,年代最老),最大的是Qhl(连圩组,年代最新);
步骤43:读取地层界线集合E中的第i条界线要素Ea/b(i为循环变量,i∈[1,k],初值为1,k为地层界线的条数)所对应的地层面要素sa、sb的岩石地层单位Rocknamea、Rocknameb。本实例中分别以地层界线集合E中的第227条界线要素E85/275所对应的地层面要素s85(清水组上段)和s275(新港粘土下段)的岩石地层单位S1q2和Qp3x1,与第732条界线要素E366/394所对应的地层面要素s366(大姑泥砾)、s394(星子岩群:归宗寺变粒岩)的岩石地层单位Qp2d和Pt1xg为例,执行步骤44-48;
步骤44:读取地层面要素s85(清水组上段)和s275(新港粘土下段)的“岩类”属性。地层面要素s85(清水组上段)和s275(新港粘土下段)的“岩类”均为非侵入体,跳转步骤45;而地层面要素s366(大姑泥砾)的“岩类”为非侵入体,s394(星子岩群:归宗寺变粒岩)的“岩类”为侵入体,则跳转步骤47;
步骤45:整合接触和不整合接触的判断。读取岩石地层单位S1q2和Qp3x1在序列表X上的序号分别为29、58,∣29-58∣=29≠1,则地层面要素s85和s275的接触关系为不整合接触,跳转步骤46,进行下一步判断;
步骤46:平行不整合接触和角度不整合接触判断。地层面要素s85(清水组上段)、s275(新港粘土下段)的接触关系进行下一步判断如下:
a)加载庐山产状图层数据(图2),并读取地层面要素s85的走向Strike85=“岩层正常走向”、倾向Dip85=“SE”和倾角Angle85=25°;地层面要素s275的走向Strike275=“岩层正常走向”、倾向Dip275=“SE”和倾角Angle275=20°;
b)Strike85=Strike275且Dip85=Dip275且∣Angle85-Angle275∣=5≤δ(本实例中δ取5°)成立,则地层面要素s85(清水组上段)、s275(新港粘土下段)的接触关系为平行不整合接触。
步骤47:沉积接触与侵入接触的判断。读取岩石地层单位Qp2d(大姑泥砾)和Pt1xg(星子岩群:归宗寺变粒岩)在序列表X上的序号分别为54、1。由于地层面要素s394(星子岩群:归宗寺变粒岩)的“岩类”属性为侵入体,且x366>x394,则为沉积接触;
步骤48:记录地层线要素的图元编号85和275和平行不整合接触关系属性和地层线要素的图元编号366和394和沉积接触关系属性。循环步骤43-48,识别地层界线集合E={E1/2,E1/58,...E417/418}中926段地层界线要素两侧地层间的接触关系。识别结果为整合接触185段,平行不整合接触36段,角度不整合6段,不整合接触663段,断层接触29段,侵入接触4段,沉积接触3段,并以专题图的形式展示,如图8所示,完成地层界线图层的制作。
Claims (5)
1.一种地层界线图层的自动化生成方法,其特征在于,该方法包括如下步骤:
(1)基于沉积、火山地层界线图层和变质地层单位图层,搜寻出图层中所有邻接地层,构建地层邻接关系图;
(2)基于地层邻接关系图提取相邻地层的公共边界,生成地层界线;
(3)基于断层图层,识别断层接触关系;
(4)基于产状图层,识别基本接触关系,完成地层界线图层的自动化生成。
2.根据权利要求1所述的一种地层界线图层的自动化生成方法,其特征在于,所述步骤(1)的具体步骤如下:
步骤11:加载沉积和火山地层单位图层、变质地层单位图层数据,并将相应的地层面要素读入集合S={si|i=1,2,…,n},n为地层面要素的个数;
步骤12:读出每个地层面要素si的顶点集vi={Pj|j=1,2,…,m},其中Pj表示构成地层面要素si的顶点,m为顶点的个数;并对每个地层面要素si的顶点集vi进行如下的预处理:
a)将顶点集vi中的点Pj按顺时针排序;
b)判断每个顶点集vi中是否存在P1=Pm,即首节点等于尾节点,若存在,则跳转下一步c),对顶点集vi进行处理;否则,跳转步骤13;
c)删除顶点集vi的尾节点Pm,保证首尾不相接,即vi={Pj|j=1,2,…,m-1};
步骤13:对图层中地层面要素的邻接关系进行初判:若地层面要素sa的顶点集va的最小外接矩形Ra和要素sb的顶点集vb的最小外接矩形Rb有交点,则跳转步骤14;否则,跳转步骤16;其中,a为循环变量,a∈[1,n-1],初值为1;b为循环变量,b∈[a+1,n],初值为a+1;
步骤14:对两地层面要素sa和sb的邻接关系进行细判,具体过程如下:
a)计算顶点集va和vb的公共点集T={p1,p2,...,pk},k为公共点的个数,并记录T中点在顶点集va中序号,记为N={n1,n2,...,nk};
b)若公共点集T中点的个数k≥2,则跳转下一步c);否则,跳转步骤16;
c)判断公共点序号集N内是否存在连续的点号nt,nt+1满足∣nt-nt+1∣=1,若存在,即地层面要素s1和s2存在邻接关系,则跳转步骤15;否则,跳转步骤16;
步骤15:构建地层邻接关系图的边d:计算顶点集va与vb的重心点Ca和Cb,点Ca和Cb作为地层邻接关系图的顶点,连接点Ca和Cb作为地层邻接关系图的边d;
步骤16:若循环变量a>n-1,则输出SARG={d1,d2,...dw},w为地层邻接关系图SARG的边的总个数;否则,循环变量b加1,此时若b≤n,则直接跳转步骤13;否则循环变量a加1,循环变量b=a+1,再跳转步骤13。
3.根据权利要求2所述的一种地层界线图层的自动化生成方法,其特征在于,所述步骤(2)的具体步骤如下:
步骤21:读取SARG={d1,d2,...dw}中第i条边di,i为循环变量,初值为1;
步骤22:求出地层面要素sa和sb的顶点集va={p1,p2,...,pm}和vb={q1,q2,...,qn},m、n分别为地层面要素sa、sb的顶点个数,并对顶点集va和vb进行如下的预处理:
a)将顶点集va={p1,p2,...,pm}和vb={q1,q2,...,qn}中的点均按顺时针排序;
b)判断顶点集va和vb是否同时满足p1=pm与q1=qn,即首节点等于尾节点;若满足,则跳转步骤23;否则,对不符合条件的顶点集va或vb,执行下一步c);
c)对顶点集va或vb继续添加第一个点p1或q1,保证首尾相接;
步骤23:计算顶点集va和vb的公共点集T={t1,t2,...,tj},j为公共点的个数,并记录公共点集T中各点分别在顶点集va和vb中的序号,分别记为Qa和Qb;
步骤24:公共点集T的分段处理,分段处理步骤如下:
a)设k为循环变量,k∈[1,j-1],初值为1;读取序号Qa中第k个点号i'、第k+1个点号j';同样,读取序号Qb中第k个点号ii'、第k+1个点号jj',若同时满足∣i'-j'∣=1且∣ii'-jj'∣=1,则点号i'和点号j'、以及点号ii'和点号jj'所对应的两点属于同一段公共边界;
b)循环变量k加1,若k<j,继续执行上一步a)直至公共点集T分段结束,并依次将属于同一段边界的点构成边e,得到公共边界Ea/b={e1,e2,..,em'},m'表示公共边界的段数;如果m'>2,执行下一步c);否则跳转步骤25;
c)依次两两比较构成边et的最后点坐标tL(x,y)与构成边et+1的第一个点坐标tF(x,y),t∈[1,m'],若存在tL(x,y)=tF(x,y),则将边et与et+1合并,边et与et+1属于同一段公共边界;
步骤25:循环变量i加1,若i≤w,循环执行步骤21-24,得到整个图层的地层界线集合E={Ea/b,Ec/d,...},其中Ea/b,Ec/d分别代表邻接地层sa与sb、邻接地层sc与sd间的地层界线。
4.根据权利要求3所述的一种地层界线图层的自动化生成方法,其特征在于,所述步骤(3)的具体步骤如下:
步骤31:加载断层图层数据,并以断层图层中的断层线F={f1,f2,...,fm}为中心轴线,以距离参数d生成断层缓冲区FB;
步骤32:提取与断层数据中断层线F={f1,f2,...,fm}相交的地层邻接关系图的边集合D={d1,d2,...,dn},n为地层邻接关系图与断层线F相交的条数;
步骤33:搜寻边集合D中边di所对应的地层面要素Sa与地层面要素Sb间的地层界线Ea/b={e1,e2,...,em},其中i为循环变量,i∈[1,n],初值为1;
步骤34:判断断层接触:若Ea/b={e1,e2,...,em}的某段界线ej在断层缓冲区FB内,则界线ej两侧地层间的接触关系即为断层接触关系;
步骤35:循环变量i加1,若i≤n,则重复执行步骤33、34,直至图层中所有的断层接触提取完毕。
5.根据权利要求4所述的一种地层界线图层的自动化生成方法,其特征在于,所述步骤(4)的具体步骤如下:
步骤41:基于地层岩性信息,区分侵入体与非侵入体;将每个地层的区分结果记录到沉积和火山地层单位图层、变质地层单位图层的地层面要素的“岩类”属性中;
步骤42:基于岩石地层单位简表,生成序列表X,且较老的地层序号数值较小;
步骤43:读取地层界线集合E中的第i条界线要素Ea/b所对应的地层面要素sa和sb的岩石地层单位Rocknamea和Rocknameb,其中i为循环变量,i∈[1,k],初值为1,k为地层界线的条数;
步骤44:读取地层面要素sa和sb的“岩类”属性;若地层面要素sa和sb的“岩类”均为非侵入体,则跳转步骤45;否则,跳转步骤47;
步骤45:整合接触和不整合接触的判断;读取岩石地层单位Rocknamea和Rocknameb在序列表X上的序号xa、xb,若序号xa、xb相邻或相同,即xa=xb或者∣xa-xb∣=1,则地层面要素sa和sb的接触关系为整合接触;否则,为不整合接触,跳转步骤46,进行下一步判断;
步骤46:平行不整合接触和角度不整合接触判断,具体判断方法如下:
a)加载产状图层数据,并读取地层面要素sa、sb的产状要素:走向Strikea和Strikeb、倾向Dipa和Dipb、以及倾角Anglea和Angleb;若任一产状要素信息无法获取,则地层面要素sa和sb的接触关系仅判断为不整合接触,跳转步骤48;
b)判断Strikea=Strikeb且Dipa=Dipb且∣Anglea-Angleb∣≤δ是否成立,其中0°<δ<10°;若成立,则地层面要素sa和sb的接触关系为平行不整合接触;否则,为角度不整合接触;
步骤47:沉积接触与侵入接触的判断;读取岩石地层单位Rocknamea和Rocknameb在序列表X上的序号xa、xb;当地层面要素sa的“岩类”属性为侵入体时,若xa≥xb,则为侵入接触;否则,则为沉积接触;同样,当地层面要素sb的“岩类”属性为侵入体时,且xa<xb,则为侵入接触;否则,则为沉积接触;
步骤48:记录地层线要素的图元编号和接触关系属性;并循环步骤43-48,直至地层界线集合E中所有界线要素两侧地层间的接触关系均识别结束,完成地层界线图层的制作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023762.2A CN105701848B (zh) | 2016-01-14 | 2016-01-14 | 一种地层界线图层的自动化生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023762.2A CN105701848B (zh) | 2016-01-14 | 2016-01-14 | 一种地层界线图层的自动化生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105701848A true CN105701848A (zh) | 2016-06-22 |
CN105701848B CN105701848B (zh) | 2018-07-20 |
Family
ID=56227367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610023762.2A Expired - Fee Related CN105701848B (zh) | 2016-01-14 | 2016-01-14 | 一种地层界线图层的自动化生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105701848B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106384016A (zh) * | 2016-10-08 | 2017-02-08 | 南京师范大学 | 单一倾向断层的自动反演方法 |
CN106897445A (zh) * | 2017-03-01 | 2017-06-27 | 武汉科技大学 | 一种基于ArcGIS的改进边界划定操作的新型方法 |
CN106934357A (zh) * | 2017-02-28 | 2017-07-07 | 南京师范大学 | 一种平行断层的自动识别方法 |
CN107065009A (zh) * | 2017-04-18 | 2017-08-18 | 中国石油大学(华东) | 一种基于地震资料的潜山油气田地下地质图编制方法 |
CN107146283A (zh) * | 2017-05-10 | 2017-09-08 | 南京师范大学 | 一种层状岩质边坡类型的自动划分方法 |
CN111583776A (zh) * | 2020-04-28 | 2020-08-25 | 南京师范大学 | 一种侵入岩体发育时序获取方法 |
CN111681314A (zh) * | 2020-05-29 | 2020-09-18 | 中国地质大学(武汉) | 一种基于虚拟等值线约束的第四系覆盖层三维建模方法 |
CN112231423A (zh) * | 2020-09-22 | 2021-01-15 | 南京师范大学 | 基于地质图的断层构造自动化恢复方法及装置 |
CN113539051A (zh) * | 2021-06-24 | 2021-10-22 | 南京师范大学 | 基于地质图的地层界线逐点岩层产状获取方法及装置 |
CN115544196A (zh) * | 2022-10-12 | 2022-12-30 | 兰州天瑞测绘工程有限公司 | 一种权属界线及走向说明自动生成方法及*** |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW581972B (en) * | 2002-05-06 | 2004-04-01 | Newsoft Technology Corp | Method to determine adjacency relation |
CN101303414B (zh) * | 2008-05-22 | 2011-03-30 | 北京航空航天大学 | 一种基于水平集的地层面及地质体生成方法 |
CN102651143A (zh) * | 2011-02-25 | 2012-08-29 | 中冶沈勘工程技术有限公司 | 一种尖灭地质体三维复杂剖面自动生成方法 |
CN104122602A (zh) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | 一种获取沉积盆地剥蚀厚度的方法及装置 |
-
2016
- 2016-01-14 CN CN201610023762.2A patent/CN105701848B/zh not_active Expired - Fee Related
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106384016B (zh) * | 2016-10-08 | 2019-04-23 | 南京师范大学 | 单一倾向断层的自动反演方法 |
CN106384016A (zh) * | 2016-10-08 | 2017-02-08 | 南京师范大学 | 单一倾向断层的自动反演方法 |
CN106934357A (zh) * | 2017-02-28 | 2017-07-07 | 南京师范大学 | 一种平行断层的自动识别方法 |
CN106934357B (zh) * | 2017-02-28 | 2020-04-03 | 南京师范大学 | 一种平行断层的自动识别方法 |
CN106897445A (zh) * | 2017-03-01 | 2017-06-27 | 武汉科技大学 | 一种基于ArcGIS的改进边界划定操作的新型方法 |
CN107065009A (zh) * | 2017-04-18 | 2017-08-18 | 中国石油大学(华东) | 一种基于地震资料的潜山油气田地下地质图编制方法 |
CN107146283A (zh) * | 2017-05-10 | 2017-09-08 | 南京师范大学 | 一种层状岩质边坡类型的自动划分方法 |
CN107146283B (zh) * | 2017-05-10 | 2020-05-19 | 南京师范大学 | 一种层状岩质边坡类型的自动划分方法 |
CN111583776A (zh) * | 2020-04-28 | 2020-08-25 | 南京师范大学 | 一种侵入岩体发育时序获取方法 |
CN111583776B (zh) * | 2020-04-28 | 2022-02-11 | 南京师范大学 | 一种侵入岩体发育时序获取方法 |
CN111681314A (zh) * | 2020-05-29 | 2020-09-18 | 中国地质大学(武汉) | 一种基于虚拟等值线约束的第四系覆盖层三维建模方法 |
CN112231423A (zh) * | 2020-09-22 | 2021-01-15 | 南京师范大学 | 基于地质图的断层构造自动化恢复方法及装置 |
CN112231423B (zh) * | 2020-09-22 | 2024-02-23 | 南京师范大学 | 基于地质图的断层构造自动化恢复方法及装置 |
CN113539051A (zh) * | 2021-06-24 | 2021-10-22 | 南京师范大学 | 基于地质图的地层界线逐点岩层产状获取方法及装置 |
CN115544196A (zh) * | 2022-10-12 | 2022-12-30 | 兰州天瑞测绘工程有限公司 | 一种权属界线及走向说明自动生成方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN105701848B (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105701848A (zh) | 一种地层界线图层的自动化生成方法 | |
CN107358654B (zh) | 基于多边形变形技术的剖面重构三维表面建模方法及*** | |
CN109903383B (zh) | 一种采煤机在工作面煤层三维模型中精确定位方法 | |
CN111899341A (zh) | 长大带状线性工程三维地质快速建模方法 | |
US20120166160A1 (en) | Block model constructing method for complex geological structures | |
CN116152461B (zh) | 地质建模方法、装置、计算机设备及计算机可读存储介质 | |
CN111311746B (zh) | 一种基于钻孔数据的智能三维地质建模方法 | |
CN110163965B (zh) | 一种单斜岩层构造地貌的实体三维模型构建方法 | |
CN105321204B (zh) | 一种三维地质栅格体模型构建方法 | |
CN104504754B (zh) | 一种油气储层多点统计建模的方法及装置 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN109979011B (zh) | 基于多源异构数据的平原地区三维地质模型建设方法 | |
CN112734926B (zh) | 面向松散层覆盖区域的图切地质剖面自动生成方法 | |
CN109859312A (zh) | 一种基于bim技术的精细化三维地质模型建模方法 | |
Gold et al. | Outside-in: an alternative approach to forest map digitizing | |
CN107766978A (zh) | 不规则井网的智能优化方法 | |
CN108986213B (zh) | 一种基于叠置技术的三维地层建模方法 | |
CN109712239A (zh) | 一种矿床精细三维地质建模方法 | |
CN109410183A (zh) | 基于点云数据的平面提取方法、***、装置及存储介质 | |
CN115113299A (zh) | 智能标志地层定位划分方法、装置、设备及存储介质 | |
CN111899316B (zh) | 一种地层结构的二维地层网格数字化编辑方法及装置 | |
CN112581558B (zh) | 图切地质剖面中侵入构造的模型构建方法及*** | |
CN112562080B (zh) | 一种基于钻孔数据的地质结构降维模型建模方法 | |
CN115294290A (zh) | 一种地质矿产勘查提取数据的方法 | |
CN114399605A (zh) | 一种基于虚拟地形面的地层实体自动建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180720 Termination date: 20210114 |