判断古主应力方向的方法及电子设备与流程

1.本发明属于基础地质研究领域，具体涉及一种利用地层缩短量定量判 断古主应力方向的方法及电子设备。


背景技术：

2.古应力场的分布基本控制了古地质历史时期的构造面貌，断裂的走向 及性质、褶皱的轴迹也主要受控于古应力的方向，在断裂的多期活动与改 造过程中，古应力场的演化也决定了断裂的最终性质和形态。现有的判断 构造应力方向的方法主要有地貌学及测量学方法、野外观测法、浅部应力 测量法、有限元法模拟、显微构造分析、震源机制解法等方法。
3.①
地貌测量法，主要通过直接观测和卫星、航测等技术对构造应力产 生的构造地貌进行测量分析，对区域水系、山系、山脊沟谷、地表及地下 裂隙等构造行迹以及地形变进行测量、统计，并结合多种数学方法进行分 析，以期得到应力场特征(王强,虎雄林，朱荣欢,等.构造应力场研究综述[j]. 云南大学学报(自然科学版),2014,36(s1)：122-129)。
②
野外观测法，通 过对野外所观测到的岩石中已存在的构造变形特征来反推变形作用发生时 的构造应力状态,其中节理、构造线理、岩墙、褶皱、擦痕、断层等是最常 用的构造变形，通常需对所得到的构造变形数据进行预处理，如利用共轭 节理判断古应力场时，需先将其所在地层产状恢复至水平(目前多采用相 关构造节理恢复软件，如johannnes duyster 1999年所编写的stereonett2.4)； 或者某些构造变形数据可直接指示古构造应力方向，如基性岩墙群的产状 及时代或褶皱两翼的平均产状，将处理后的数利用赤平投影可求取最大主 应力方向(万天丰.1988.古构造应力场[m].北京:地质出版社:4
–
10；关 宝文,郑建京,赵恒,等.博格达山地区二叠纪以来构造应力场解析及地质 意义[j].大地构造与成矿学,2017,41(1):42-49；吴志远,杨德芳，马丽红,等. 大宁-吉县地区古构造应力场恢复[j].华北科技学院学报,2018,15(2):43-48)。
ꢀ③
显微构造分析法，在应力作用下，构成岩石的矿物等也会发生变形或表 现出定向排列的现象，也可见到有微观尺度的破裂现象，这些现象统称为 微观构造。可通过对微观构造如矿物的定向构造、裂隙、带状构造、显微 褶皱的统计，绘制极密图，进而判断主应力方向(孟献真,王邵祥.甘肃文县 阳山金矿带显微构造研究[j].安徽地质,2013,23(1)：40-44)。以上几种方法主 要基于安德森模式，是目前常用的应力方向判别方法，但是鉴于所应用的 数据均为现今构造面貌或形态，受限于所能观测到的构造，仅能判断多期 构造叠加下的最终变形的应力方向和地表构造应力方向。
④
浅部应力测量 法，包括平面变形法、原地应力测量方法、恢复法、孔壁压裂法(钻孔崩落)、 岩性换算法(波速、磁差、电阻率、同位素)、x射线法等等(秦向辉,陈群 策,郝跃进,等.安徽金寨沙坪沟矿区地应力测量与现今应力场研究[j].地球 学报,2016,37(2):223-231；)。此类方法主要是通过传感器等仪器、设备直接 测得地表或近地表(井下)的应力绝对值，或通过测得的应变值转换为应 力值。该种方法不仅能够得到主应力的方向，还可得到应力的绝对值，但 是该方法多应用于局部应力方向的判断，且受限于井数据或者矿区的分布。
ꢀ⑤
有
限元模拟法，多采用各类模拟软件如anysys等，该方法流程是先建 立目标区地质模型，选取力学参数并将模型划分为有限个单元(二维模型 划分为任意多边形，三维模型划分为多面体)并划分单元网格密度，而后 添加边界条件，在分析计算后，将结果可视化。有限元模拟法受限于技术 手段和方法，所得到的结果受人为因素干扰较大，模拟准确性受到地质模 型精度、边界条件、网格密度等多条件限制。(胡秋媛，李理.鲁西地区晚 中生代-古近纪伸展构造的应力场数值模拟[j].石油实验地质,2015，37(2)： 259-266；郑淳方,候贵廷,詹彦，等.库车坳陷新生代构造应力场恢复[j]. 2016,35(1)：130-139)。
⑥
震源机制解法，震源机制解是对地震做出合理性 解释的一种结果，用于判断震后断面走向、节面、倾角等参数，利用震源 机制解求取构造应力场时，首先将记录到的各台站地震数据汇集处理成可 读取的格式求取断层面解、震源深度和矩震级，以震源机制解为基础数据， 将研究区域划分为网格，使用阻尼线性逆推法计算区域构造应力场(崔子 健,陈章立,王勤采,等.南北地震带区域构造应力场反演[j].地震学报, 2019，41(2)：219-229；吴宝峰,郝永梅，徐建权，等.中俄交界4.8级 地震综合震源机制解及构造应力场[j].地震地磁观测与研究,2013, 34(1/2):31-35)。因震源机制解法基于天然地震，在应力场的判别时受限于 天然地震的发生频率和区域，其应用起来区域局限性较大。
[0004]
以上应力场大小或方向的求取及判别均有其各自的优越性，在实际地 质问题的解决中得到了应用，然而这些方法多基于现今构造形态给出主应 力特征，在处理多期、多方向应力场叠加的构造时具有一定局限性，且仅 能定性说明应力场的大致方向，存在较大不准确性，尚未有方法进行定量、 相对准确判断各期古应力方向。
[0005]
因此，特别需要一种方法精确的判断应力场方向角度，准确反映不同 期次构造应力方向及构造活动强度。


技术实现要素：

[0006]
本发明的目的是提出一种精确的判断应力场方向角度，准确反映不同 期次构造应力方向及构造活动强度。
[0007]
本发明提供一种利用地层缩短量定量判断古主应力方向的方法，包 括：选取研究工区的边界，根据所述边界确定所述研究工区的中心点；在 所述研究工区的三维地震数据中，以所述中心点为圆心，在预设角度范围 内按预设间隔角度提取通过所述中心点的多个剖面；划分所述研究工区地 质的构造层，获得多个构造子层，并确定每个所述构造子层的标志层；计 算每个所述构造子层的标志层在每个剖面上的净缩短量；分别针对每个构 造子层，基于所述构造子层的标志层在每个剖面上的净缩短量，绘制所述 构造子层对应的缩短量玫瑰花图，基于所述构造子层对应的缩短量玫瑰花 图，确定所述构造子层对应的主应力方向。
[0008]
可选的，所述划分所述研究工区地质的构造层，获得多个构造子层包 括：根据所述研究工区的区域构造演化史、地层接触关系及纵向构造形态 差异，划分所述研究工区地质的构造层，获得多个构造子层。
[0009]
可选的，采用下述公式计算构造子层的标志层在剖面上的净缩短量：
[0010]
l
nrm
＝lnrs
m-lnrs
m-1
[0011]
其中，l
nrm
为第m个构造子层的标志层在第n个剖面上的净缩短量， n＝1,2,3
……
n，m＝2，3
……
m，n表示剖面的数量，m表示构造子层的数 量，lnrsm为第m个构造子层的标志层在第n个剖面上的总缩短量，lnrs
m-1
为第m-1个构造子层的标志层在第n个剖面上的总缩短量，第m-1个构造 子层位于第m个构造子层的上面。
[0012]
可选的，采用下述步骤计算构造子层的标志层在剖面上的总缩短量： 根据剖面及层位追踪、断层解释结果，测量所述构造子层的标志层在所述 剖面上的现今剖面长度；在所述构造子层的所述剖面上，将所述构造子层 的标志层拉直，获得所述构造子层的标志层在所述剖面上的原始层长；基 于所述构造子层在所述剖面上的所述现今剖面长度和原始层长，获得所述 构造子层的标志层在所述剖面上的总缩短量。
[0013]
可选的，采用下述公式计算构造子层的标志层在剖面上的总缩短量：
[0014]
l
nr＝
l
ny-l
nj
[0015]
其中，lnr为所述构造子层的标志层在第n个剖面上的总缩短量，l
ny
为所述构造子层的标志层在第n个剖面上的原始层长，l
nj
为所述构造子层 的标志层在第n个剖面上的现今剖面长度。
[0016]
可选的，所述基于所述构造子层的标志层在每个剖面上的净缩短量， 绘制所述构造子层对应的缩短量玫瑰花图包括：以所述中心点为圆心，以 剖面的方位角为角度坐标，以净缩短量为数量大小坐标绘制坐标图，将所 述构造子层的标志层在每个剖面上的净缩短量标注在所述坐标图中，将全 部标注的点依次连接起来，获得所述构造子层对应的缩短量玫瑰花图。
[0017]
可选的，所述基于所述构造子层对应的缩短量玫瑰花图，确定所述构 造子层对应的主应力方向包括：基于所述构造子层对应的缩短量玫瑰花图， 将最大净缩短量对应的剖面的方位角作为所述构造子层对应的主应力方 向。
[0018]
可选的，将所述构造子层内保留齐全、层位稳定的层作为所述构造子 层的标志层。
[0019]
可选的，所述预设角度范围为0-180
°
。
[0020]
本发明还提供一种电子设备，所述电子设备包括：存储器，存储有可 执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以 实现上述利用地层缩短量定量判断古主应力方向的方法。
[0021]
本发明的有益效果在于：本发明的利用地层缩短量定量判断古主应力 方向的方法用地层缩短量定量判断古主应力方向的方法通过计算各构造子 层各向地层定量缩短量，绘制地层缩短量玫瑰花图，可判断各构造子层最 大缩短量方向，进而较准确判断古应力场方向，较为准确反映不同期次构 造应力方向及构造活动强度，为构造恢复、断裂成因及演化、古石油地质 条件研究提供较精确依据。
[0022]
本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图 和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随 后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本 发明的特定原理。
附图说明
[0023]
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的 上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性 实施方式中，相同的参考标
号通常代表相同部件。
[0024]
图1示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的流程图。
[0025]
图2示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的研究区地震三维工区范围及中心点。
[0026]
图3示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的某研究区的剖面选取情况。
[0027]
图4示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的某研究区第1构造子层的缩短量玫瑰花图。
[0028]
图5示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的某研究区第2构造子层的缩短量玫瑰花图。
具体实施方式
[0029]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明 的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这 里阐述的实施方式所限制。
[0030]
本发明提供一种利用地层缩短量定量判断古主应力方向的方法，包 括：选取研究工区的边界，根据边界确定研究工区的中心点；在研究工区 的三维地震数据中，以中心点为圆心，在预设角度范围内按预设间隔角度 提取通过中心点的多个剖面；划分研究工区地质的构造层，获得多个构造 子层，并确定每个构造子层的标志层；计算每个构造子层的标志层在每个 剖面上的净缩短量；分别针对每个构造子层，基于构造子层的标志层在每 个剖面上的净缩短量，绘制构造子层对应的缩短量玫瑰花图，基于构造子 层对应的缩短量玫瑰花图，确定构造子层对应的主应力方向。
[0031]
具体的，古应力场控制断裂及构造地貌的形成与演化，并进一步控制 了古成藏期油气藏的石油地质条件，鉴于现近地质面貌多为多期应力场下 多期构造运动改造的结果，厘定单一期次构造运动特征对探究古成藏期石 油地质条件尤为重要。
[0032]
古构造应力场是控制各期古构造变形的直接因素，根据安德森模式， 在主应力方向上，地层应具有最大缩短(伸展)量。以安德森模式为理论 依据，认为在考虑区域应力方向的判断时，可忽略同一构造层中水平方向 岩石的不均一性，沿最大主应力方向，地层应具有最大缩短量，通过剥离 不同构造期次的构造活动。在假设盆地演化过程中层长守恒的前提下，通 过求取各构造层或各层在不同方向上的缩短量，进而对比分析缩短量，将 缩短量最大的方向定义为最大主应力方向，进而可以恢复研究区最大古主 应力方向，并可以各期地层缩短量为据，判断各期研究区内古应力相对大 小。首先，在研究区三维地震资料解释基础之上，确定研究区边界及中心 点，以中心点为圆心，在0-180
°
之间，按20
°
间隔(可根据需求精度加大 或减小间隔)，选取剖面；而后对所选取的剖面由新至老主要界面进行非塑 性地层平衡恢复，在恢复过程中，应剥除下覆地层在后期应力场作用下的 缩短量，在得到各主要界面缩短量后，最终可绘制各层地层缩短量玫瑰花 图，判断各古地质历史时期古应力场方向。
[0033]
研究区边界的选取应不包含盆地边界造山带，中心点的选取应尽量保 证距离各
边界尽量相等。
[0034]
以所选取的工区中心点为圆心，在0-180
°
范围内按等角度间隔(如 20
°
)选取过中心点剖面，将剖面编号为l1，l2，
……
ln,在选取剖面时， 剖面为直线，且保持各剖面长度尽量等长。
[0035]
根据示例性的实施方式，利用地层缩短量定量判断古主应力方向的方 法通过计算各构造子层各向地层定量缩短量，绘制地层缩短量玫瑰花图， 可判断各构造子层最大缩短量方向，进而较准确判断古应力场方向，较为 准确反映不同期次构造应力方向及构造活动强度，为构造恢复、断裂成因 及演化、古石油地质条件研究提供较精确依据。
[0036]
作为可选方案，划分研究工区地质的构造层，获得多个构造子层包括： 根据研究工区的区域构造演化史、地层接触关系及纵向构造形态差异，划 分研究工区地质的构造层，获得多个构造子层。
[0037]
具体的，根据研究区区域构造演化史、地层接触关系、纵向构造形态 差异等，划分构造层，将不同期次构造变形层进行划分，并选取各构造层 内保留较全、层位稳定(同相轴连续性好、反射强)的层作为该构造层的 标志层。
[0038]
作为可选方案，采用下述公式计算构造子层的标志层在剖面上的净缩 短量：
[0039]
l
nrm
＝lnrs
m-lnrs
m-1
[0040]
其中，l
nrm
为第m个构造子层的标志层在第n个剖面上的净缩短量， n＝1,2,3
……
n，m＝2，3
……
m，n表示剖面的数量，m表示构造子层的数 量，lnrsm为第m个构造子层的标志层在第n个剖面上的总缩短量，lnrs
m-1
为第m-1个构造子层的标志层在第n个剖面上的总缩短量，第m-1个构造 子层位于第m个构造子层的上面。
[0041]
具体的，针对一个构造子层，计算该构造子层在一个剖面中的净缩短 量，例如l
nrm
为第m个构造子层在第n个剖面上的净缩短量，先测量第 m个构造子层在第n个剖面上的原始长度，再获得第m个构造子层在第n 个剖面上的现今剖面长度，通过原始长度和现今剖面长度获得第m个构造 子层在第n个剖面上的总缩短量lnsrm，再剔除第m个构造子层上面相邻 的第m-1个构造子层在第n个剖面上的总缩短量，将这两个总缩短量做差， 即获得该构造子层在该剖面上的净缩短量lnrs
m-1
，即获得了第m个构造子 层的标志层在第n个剖面上的净缩短量。采用同样的方法，分别计算每个 构造子层上在每个剖面上的净缩短量。
[0042]
针对第1个构造子层，分别只需计算第1个构造子层在各剖面中的总 缩短量，将其总缩短量作为其净缩短量lnr1(n为剖面编号，n＝1,2,3
……
)。
[0043]
作为可选方案，采用下述步骤计算构造子层的标志层在剖面上的总缩 短量：根据剖面及层位追踪、断层解释结果，测量构造子层的标志层在剖 面上的现今剖面长度；在构造子层的剖面上，将构造子层的标志层拉直， 获得构造子层的标志层在剖面上的原始层长；基于构造子层在剖面上的现 今剖面长度和原始层长，获得构造子层的标志层在剖面上的总缩短量。
[0044]
具体的，针对一个构造子层的一个剖面，根据所选剖面及层位追踪、 断层解释结果，测量、计算该剖面的现今剖面长度，再测算现今一套构造 子层在该剖面中的原始层长，获得了该剖面的原始层长和现今剖面长度， 将该剖面的原始层长减去现今剖面长度，即获得该构造子层的该剖面的总 缩短量。采用同样的方法，获得每个构造子层在每个剖面上的总缩短量。
[0045]
作为可选方案，采用下述公式计算构造子层的标志层在剖面上的总缩 短量：
[0046]
l
nr＝
l
ny-l
nj
[0047]
其中，lnr为构造子层的标志层在第n个剖面上的总缩短量，l
ny
为构 造子层的标志层在第n个剖面上的原始层长，l
nj
为构造子层的标志层在第 n个剖面上的现今剖面长度。
[0048]
具体的，针对一个构造子层的一个剖面，将该剖面的原始层长减去现 今剖面长度，即获得总缩短量。采用同样的方法，获得每个构造子层在每 个剖面上的总缩短量。
[0049]
作为可选方案，基于构造子层的标志层在每个剖面上的净缩短量，绘 制构造子层对应的缩短量玫瑰花图包括：以中心点为圆心，以剖面的方位 角为角度坐标，以净缩短量为数量大小坐标绘制坐标图，将构造子层的标 志层在每个剖面上的净缩短量标注在坐标图中，将全部标注的点依次连接 起来，获得构造子层对应的缩短量玫瑰花图。
[0050]
具体的，针对一个构造子层，以中心点为圆心，以剖面方位角θ为角 度坐标，以净缩短量为数量大小坐标绘制坐标图，将该构造子层在各剖面 中计算的净缩短量标注在坐标图中，将全部标注的点依次连接起来，获得 该构造子层对应的缩短量玫瑰花图。用同样的方法，可以获得每个构造子 层对应的缩短量玫瑰花图。
[0051]
作为可选方案，基于构造子层对应的缩短量玫瑰花图，确定构造子层 对应的主应力方向包括：基于构造子层对应的缩短量玫瑰花图，将最大净 缩短量对应的剖面的方位角作为构造子层对应的主应力方向。
[0052]
具体的，分别针对每个构造子层，根据其对应的缩短量玫瑰花图，将 最大净缩短量对应的剖面的方位角作为该构造子层对应的主应力方向。通 过玫瑰花图结果可得到各期主应力方向，并可根据不同构造层的玫瑰花图 对比结果判断各期相对应力大小。
[0053]
作为可选方案，将构造子层内保留齐全、层位稳定的层作为构造子层 的标志层。
[0054]
具体的，选取各构造子层内保留较全、层位稳定(同相轴连续性好、 反射强)的层作为该构造子层的标志层。
[0055]
作为可选方案，预设角度范围为0-180
°
。
[0056]
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行 指令；处理器，处理器运行存储器中的可执行指令，以实现上述利用地层 缩短量定量判断古主应力方向的方法。
[0057]
实施例一
[0058]
图1示出了根据本发明的一个实施例的一种利用地层缩短量定量判断 古主应力方向的方法的流程图。图2示出了根据本发明的一个实施例的一 种利用地层缩短量定量判断古主应力方向的方法的研究区地震三维工区范 围及中心点。图3示出了根据本发明的一个实施例的一种利用地层缩短量 定量判断古主应力方向的方法的某研究区的剖面选取情况。图4示出了根 据本发明的一个实施例的一种利用地层缩短量定量判断古主应力方向的方 法的某研究区第1构造子层的缩短量玫瑰花图。图5示出了根据本发明的 一个实施例的一种利用地层缩短量定量判断古主应力方向的方法的某研究 区第2构造子层的缩短量玫瑰花图。
[0059]
结合图1、图2、图3、图4和图5所示，该利用地层缩短量定量判断 古主应力方向的方法，包括：
[0060]
步骤1：选取研究工区的边界，根据边界确定研究工区的中心点；
[0061]
步骤2：在研究工区的三维地震数据中，以中心点为圆心，在预设角 度范围内按预设间隔角度提取通过中心点的多个剖面；
[0062]
其中，预设角度范围为0-180
°
。
[0063]
步骤3：划分研究工区地质的构造层，获得多个构造子层，并确定每 个构造子层的标志层；
[0064]
其中，划分研究工区地质的构造层，获得多个构造子层包括：根据研 究工区的区域构造演化史、地层接触关系及纵向构造形态差异，划分研究 工区地质的构造层，获得多个构造子层。
[0065]
其中，将构造子层内保留齐全、层位稳定的层作为构造子层的标志层。
[0066]
步骤4：计算每个构造子层的标志层在每个剖面上的净缩短量；
[0067]
其中，采用下述公式计算构造子层的标志层在剖面上的净缩短量：
[0068]
l
nrm
＝lnrs
m-lnrs
m-1
[0069]
其中，l
nrm
为第m个构造子层的标志层在第n个剖面上的净缩短量， n＝1,2,3
……
n，m＝2，3
……
m，n表示剖面的数量，m表示构造子层的数 量，lnrsm为第m个构造子层的标志层在第n个剖面上的总缩短量，lnrs
m-1
为第m-1个构造子层的标志层在第n个剖面上的总缩短量，第m-1个构造 子层位于第m个构造子层的上面。
[0070]
其中，采用下述步骤计算构造子层的标志层在剖面上的总缩短量：根 据剖面及层位追踪、断层解释结果，测量构造子层的标志层在剖面上的现 今剖面长度；在构造子层的剖面上，将构造子层的标志层拉直，获得构造 子层的标志层在剖面上的原始层长；基于构造子层在剖面上的现今剖面长 度和原始层长，获得构造子层的标志层在剖面上的总缩短量。
[0071]
其中，采用下述公式计算构造子层的标志层在剖面上的总缩短量：
[0072]
l
nr＝
l
ny-l
nj
[0073]
其中，lnr为构造子层的标志层在第n个剖面上的总缩短量，l
ny
为构 造子层的标志层在第n个剖面上的原始层长，l
nj
为构造子层的标志层在第 n个剖面上的现今剖面长度。
[0074]
步骤5：分别针对每个构造子层，基于构造子层的标志层在每个剖面 上的净缩短量，绘制构造子层对应的缩短量玫瑰花图，基于构造子层对应 的缩短量玫瑰花图，确定构造子层对应的主应力方向。
[0075]
其中，基于构造子层的标志层在每个剖面上的净缩短量，绘制构造子 层对应的缩短量玫瑰花图包括：以中心点为圆心，以剖面的方位角为角度 坐标，以净缩短量为数量大小坐标绘制坐标图，将构造子层的标志层在每 个剖面上的净缩短量标注在坐标图中，将全部标注的点依次连接起来，获 得构造子层对应的缩短量玫瑰花图。
[0076]
其中，基于构造子层对应的缩短量玫瑰花图，确定构造子层对应的主 应力方向包括：基于构造子层对应的缩短量玫瑰花图，将最大净缩短量对 应的剖面的方位角作为构造子层对应的主应力方向。
[0077]
以塔里木盆地内部某研究区为例，应用本发明进行分析，对该地区各 主要构造活动期应力方向进行求取，取得良好的应用效果。
[0078]
以塔里木盆地内某一三维地震工区为研究范围(已远离盆地边界)，以 地震工区边界为研究边界，确定工区内地震采集数据覆盖范围大小，以数 据所覆盖范围为准，选取
中心点，如图2所示。
[0079]
以所选取的中心点为圆心，在0-180
°
之间，每间隔20
°
选择直线剖面， 并将剖面依次编号为l1，l2,
……
,l9，各剖面长度尽量保证等长，如图3 所示。
[0080]
根据盆地区域构造背景及地层接触关系等，以t74界面之上的不整合 面为界，将各剖面中地层划分为两个主要构造子层，分别以t46界面和t81 界面作为两构造子层的标志层。
[0081]
分别计算各剖面的现今剖面长度、构造子层1标志层原始层长、构造 子层1缩短量l1r1,
……
lnr1，构造子层2原始层长、构造子层2缩短量l1r2见附表1。
[0082]
按照所选剖面角度及缩短量大小，分构造层绘制缩短量玫瑰花图，可 明显判断出各期应力场方向，构造子层1基本上可代表印支-燕山期挤压应 力来源，期方向为北西120
°
左右和北东80
°
左右。如附图5所示，构造 子层2基本代表海西晚期构造运动应力来源，其主应力方向为北西140
°
左 右；此外在该方法应用的过程中，需要注意的是，在最终结果中所反映的 构造层应力方向似呈现出非单一方向的特征，即在另外的方向上也呈现出 具有较大缩短量特征，该种现象与构造层的划分有关，在具体应用时，构 造层和剖面角度间隔的划分越细，所得到的应力方向准确性越高。
[0083]
附表1不同剖面构造层1、构造层2缩短量统计
[0084][0085][0086]
实施例二
[0087]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可 执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述利用 地层缩短量定量判断古主应力方向的方法。
[0088]
根据本公开实施例的电子设备包括存储器和处理器。
[0089]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包 括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计 算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存 储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache) 等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0090]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指 令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件 以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储 器中存储的该计算机可读指令。
[0091]
本领域技术人员应能理解，为了解决如何获取良好用户体验效果的技 术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些 公知的结构也应包含在本公开的保护范围之内。
[0092]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此 不再赘述。
[0093]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽 性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范 围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更 都是显而易见的。