周报0610--0616（共中心点道集如何抽取）

现有的速度模型

• SEG/EAGE盐和Overthrust模型 是两个3D地质模型，是SEG和EAGE之间多阶段合作的一部分，由美国能源部国家实验室以及来自工业界和学术界的几个机构资助。这些模型由SEG/EAGE 3D建模委员会开发，该委员会由Fred Aminzadeh（南加州大学）和Fabio Rocca（米兰理工大学）共同主持。这两个数据集由三维地下速度模型、模型地表和补充文档组成。有关内容的完整说明，请参阅随附的 CD 小册子的附录（1997）。
  盐模型的合成地震道可用。大部分计算是在桑迪亚和洛斯阿拉莫斯国家实验室完成的。
  3M字节的数据，包含A、B、C三个阶段，是合成数据集
• Marmousi 是一个复杂的2D结构模型及其地震响应由法兰西石油研究院（IFP）设计，用于测试二维偏移算法、比较深度偏移和速度确定模型。 它涉及强烈的水平和垂直速度变化。滑脱背斜下有一个油气成藏（平坦斑）。——1988法国 声学合成数据，基于非洲某盆地的海槽the North Quenguela Trough in the Quanza Basin of Angola
• Marmousi2 即AGL Elastic Marmousi 模型，由休斯顿大学创建，并由其作者 Gary S. Martin、Robert Wiley 和 Kurt Marfurt 于 2019 年经联合地球物理实验室联盟项目主任许可发布。它是原始 Marmousi 模型的弹性扩展。该模型的横向延伸为17 km，深度为3.5 km，包括总共199层地质层，以及450 m深度的延伸水层，以模拟深水环境。这些层源自与原始 Marmousi 模型相同的地层拾取，由法国石油研究所 （IFP） 的 Aline Bourgeois 提供。模拟数据是在休斯顿大学（UH）的超级计算机系统上计算的，该系统由Sun Microsystems提供支持。该数据集包括密度、P波和S波速度的地下模型，以及一系列经过处理的数据集（NMO叠后、叠前和叠后时间和深度偏移部分）。——美国LLNL 偏移量可达 15 km

Versteeg, Roelof (1994). “The Marmousi experience: Velocity model determination on a synthetic complex data set”. The Leading Edge 13 (9): 927–936. doi:10.1190/1.1437051> Martin, G. S., 2004, The Marmousi2 model, elastic synthetic data, and an analysis of imaging and AVO in a structurally complex environment: Master’s thesis. University of HoustonMartin, Gary S.; Marfurt, Kurt J.; Larsen, Shawn (2002). Marmousi‐2: An updated model for the investigation of AVO in structurally complex areas. pp. 1979–1982. doi:10.1190/1.1817083.Martin, Gary S.; Wiley, Robert; Marfurt, Kurt J. (2006). “Marmousi2: An elastic upgrade for Marmousi”. The Leading Edge 25 (2): 156–166. doi:10.1190/1.2172306. ISSN 1070-485X.
模型的宽度从9.2 km扩展到17 km，深度从3 km扩展到3.5 km，以更好地模拟深水环境中的长偏移采集。
最初的 Marmousi 模型靠近 Marmousi2 的中心。我们扩展了原始模型中的地平线以适应扩展模型，并添加了 41 个新地平线，使 Marmousi2 中的地平线总数达到 199 个。
由于对评估 AVO 分析的振幅保持感兴趣，因此增强的 Marmousi2 模型现在包括了几个结构简单但地层复杂的特征，这些特征位于原始复杂结构区域之外。我们还插入了几个储层，以明确模拟Marmousi2模型的原始部分和新部分中的通道和其他碳氢化合物。
Marmousi2 是一个完全弹性的模型，不仅支持压缩波，还支持剪切波、转换波和各种引导波，包括沿海底行进的导波。
如何定义模型：需要为每一层分配一个剪切速度以及原始压缩速度和密度。P波–与Marmousi 模型一样/ V = V0 + kz压实模型/ V0new = V0old − (k∗468)分配了与相邻单元相似的 P 波速度；
S波–应用Greenburg和Castagna（1992）变换；
密度–应用Castagna等人（1993）变换来获得密度。
为什么做这个2模型（3个目标）：提供一个真实的模型来评估叠前深度成像对 AVO 的影响；（没太理解文献中说的真实是怎个意思，为什么说Marmousi、SEAM、Overthrust是真实的，明明也是设计合成的？——前面看一篇去噪的论文，里面的ground truth是说”干净“标签的意思，也就是无噪声的）
1.希望合成数据量能够为研究人员、软件供应商、服务提供商和石油公司提供一种评估、校准和比较速度分析、地震成像、AVO和阻抗反演新发展的方法
2.使用 Marmousi2 作为更复杂模型的基础，这些模型将包括衰减、各向异性、复杂地层、更复杂的碳氢化合物流体取代公式和更具创新性的采集技术的影响
3.最终目标是提供一个足够大和复杂的模型，供研究人员评估算法和计算机架构的改进，这将为真实地震勘测的真正三维数值建模奠定基础。
为什么做一些改变
1.有点极端的值被放入模型中–以确保这些单元具有非常强的反射率响应。
2.引入了几个沙子单元–不会发生 P 波阻抗变化，横波表现出较小的横波阻抗变化。因为更强（更高的声阻抗）同相轴往往是大部分模式转换的来源。所以波阻抗小意味着波的类型不易发生转换，比如P波转换成S波。

• OpenFWI 是一个大规模多结构基准数据集的集合，以促进对FWI的多样化、严谨和可重复的研究。特别是，OpenFWI由来自多个震源的12 个数据集（总2.1 TB）组成。它涵盖了地球物理学的不同领域（界面、断层、CO2 储层等），涵盖了不同的地质地下结构（平坦、曲线等），并包含各种数量的数据样本 （2K - 67K）。它还包括一个用于 3D FWI 的数据集。——美国
• 𝔼${}^{\mathbf{F}\mathbf{W}\mathbf{I}}$ 是一个开源的弹性FWI数据集。 𝔼${}^{\mathbf{F}\mathbf{W}\mathbf{I}}$包含八个数据集（无断层的-平坦AB、弯曲AB，有断层的-平坦AB、弯曲AB），包括三种深度学习方法的基准。随𝔼${}^{\mathbf{F}\mathbf{W}\mathbf{I}}$发布的数据集提供了不同的 P 波和 S 波速度，专门解决了多参数反演中遇到的耦合问题。最初的基准测试在某些数据集上显示出有希望的结果，而其他数据集可能需要进一步研究。此外，还进行了耦合测试，以深入了解多参数反演问题的网络设计。——美国LANL

https://imageevent.aapg.org/portals/26/3896358.pdf

国内石油公司
中海油东方地球物理公司（BGP）——中石油旗下全资物探专业化子公司
SEG勘探地球物理学家学会EAGE欧洲地球科学家和工程师协会

方位角（Azimuth）是什么宽/窄方位，是为采集地震数据指示一个方向，由震源和检波器的方位角确定，想象震源是个手电筒，它和检波器决定着要照明的方向。方位角分析（by玫瑰图）反映了地震波传播方向的偏好
宽方位地震数据的加载是各向异性分析必需的

图片源自毛宁波教授的PPT，详见https://www.bilibili.com/video/BV1qK4y1V77Z/?vd_source=d6b1de7f052664abab680fc242ef9bc1
SEG-Y道头信息里没有直接提供方位角信息
一般来讲，CMP未做偏移，CDP是已经做好偏移的做好所有常规处理的
各向异性往往发生在大偏移距上
走时随着方位角的变化而变化，这就是速度的各向异性

地震数据标定，是为了建立地震数据的反射波（反射时间）和钻井资料（深度）的联系，帮助我们来了解地震反射波所对应的地质层位、地质时代、岩性组合、含油气层位等特征。

观测系统相关内容纵测线与非纵测线都是获得震源和检波点下方的地质信息。纵测线：震源和检波点在同一条直线上，不论多少个。比如二维地震勘探非纵测线：所有的震源和检波点不在同一直线上。比如三维地震勘探偏移距≠炮检距
偏移距是震源与第一个检波点之间的距离，要小于最浅的目的层
炮检距就是震源与检波点之间的距离

角道集
角道集指的是地震勘探中组织数据的一种方式。角指的是角点，通常称为“共中心点”。换句话说角道集就是共中心点道集。
共中心点道集
是从共炮点道集抽取的（是同一个震源，但是换位置激发多次，实现对地下同一点的多次覆盖）——共中心点多次覆盖观测系统

不论是共反射点道集、等炮检距道集（不同的激发位置，不同的接受位置）还是共接收点道集，同一个震源，但是换位置激发多次，实现对地下同一点的多次覆盖。
比如，针对每次都是一点激发的2D纵测线观测系统，抽取共反射点道集从综合平面图示法看：
> 注意共中心点和共反射点有区别：虽然都是激发点和检波点之间的点，但是共中心点指地表/地上的位置而共反射点指地下的位置。因此它们仅在反射界面水平时是重合的。（正是因为地下不水平时CMP和CDP不是一个，所以CMP的地震剖面都是经过常规处理后的数据，属于叠后数据。这就是“用CMP数据做反演属于叠后反演”说法的由来。）

想法＋论文
波形嵌入：把波形用自编码器（DNN+残差块）做了重建处理，去掉不相关噪声再做地层拾取的任务不会受横向地质扭曲的影响。效果比直接用波形相关性要好

Waveform embedding: Automatic horizon picking with unsupervised deep learning

下一步计划：阅读Mapping full seismic waveforms to vertical velocity profiles by deep learning，并复现代码