石油地球物理勘探  2021, Vol. 56 Issue (1): 1-8  DOI: 10.13810/j.cnki.issn.1000-7210.2021.01.001
0
文章快速检索     高级检索

引用本文 

常德宽, 雍学善, 王一惠, 杨午阳, 李海山, 张广智. 基于深度卷积神经网络的地震数据断层识别方法. 石油地球物理勘探, 2021, 56(1): 1-8. DOI: 10.13810/j.cnki.issn.1000-7210.2021.01.001.
CHANG Dekuan, YONG Xueshan, WANG Yihui, YANG Wuyang, LI Hai-shan, ZHANG Guangzhi. Seismic fault interpretation based on deep convolutional neural networks. Oil Geophysical Prospecting, 2021, 56(1): 1-8. DOI: 10.13810/j.cnki.issn.1000-7210.2021.01.001.

本项研究受国家自然科学基金项目“页岩气储层有效应力分布规律的精细地震预测方法研究”(41674130)和中国石油天然气集团有限公司科学研究与技术开发项目“深层及非常规物探新方法新技术”(2019A-3312)联合资助

作者简介

常德宽  博士研究生, 工程师, 1990年生; 2012、2015年分别获中国石油大学(华东)勘查技术与工程专业学士学位和地质资源与地质工程专业硕士学位, 2019年开始攻读该校能源与环保专业博士学位; 现为SEG和EAGE会员, 就职于中国石油勘探开发研究院, 主要从事地震资料处理、解释等领域的科研工作, 重点研究人工智能在地球物理勘探领域的应用

常德宽, 甘肃省兰州市雁儿湾路535号中国石油勘探开发研究院西北分院, 730020。Email:changdk@petrochina.com.cn

文章历史

本文于2020年2月29日收到,最终修改稿于同年9月15日收到
基于深度卷积神经网络的地震数据断层识别方法
常德宽①② , 雍学善 , 王一惠 , 杨午阳 , 李海山 , 张广智     
1 中国石油大学(华东), 山东青岛 266555;
2 中国石油勘探开发研究院西北分院, 甘肃兰州 730020
摘要:对地震数据进行断层解释一直是油气勘探开发过程中的一项重点工作。传统的断层解释主要是以人机交互方式进行的,效率低,并且人为因素可能增大断层解释结果的不确定性;而常规的断层识别方法则通常需要设置多个控制参数,导致断层识别的结果严重依赖参数设置的准确性。为此,提出一种基于卷积深度神经网络的地震数据断层识别方法,该方法利用ResNet深度残差网络可有效训练深层卷积神经网络和U-Net架构可表征多尺度、多层次特征信息的优势,将ResNet和U-Net架构联合,构建了用于地震数据断层识别的网络架构(SeisFault-Net)。其中,U-Net架构由编码和解码两个子网络组成,使SeisFault-Net以端到端的方式进行模型训练;同时,利用残差神经网络克服深层网络梯度弥散的问题,有效提高SeisFault-Net的训练效率。训练后的SeisFault-Net无需设置任何参数即可对地震数据进行断层识别,避免了常规方法中人为设置参数的经验误差和不确定性。数据实验表明,提出的SeisFault-Net方法可准确地识别断层位置,且识别的断层垂向连续性好,断层轮廓清晰。与相干算法相比,SeisFault-Net方法识别的断层细节更丰富,断层解释更准确;同时,可有效提高地震断层识别的效率。
关键词断层识别    深度学习    深度残差网络    U-Net架构    地震数据解释    
Seismic fault interpretation based on deep convolutional neural networks
CHANG Dekuan①② , YONG Xueshan , WANG Yihui , YANG Wuyang , LI Hai-shan , ZHANG Guangzhi     
1 China University of Petroleum (East China), Qingdao, Shandong 266555, China;
2 Research Institute of Petroleum Exploration & Development-Northwest, PetroChina, Lanzhou, Gansu 730020, China
Abstract: Seismic fault interpretation has always been a key task in the process of oil and gas exploration and development. Conventional fault interpretation is mainly based on human-computer interaction, which is of low efficiency and causes the results with many uncertainties. In addition, conventional methods for fault interpretation usua-lly set multiple parameters, whose controls accuracy of the predicted faults. This paper proposes a method using seismic data based on convolutional deep neural networks. Taking the advantages of ResNet for effectively training deep convolutional neural network and U-Net architecture for characterizing multi-scale and multi-layer characteristic information, this method combines deep residual neural network and U-Net architecture to construct a network architecture (SeisFault-Net) for fault interpretation based on seismic data. The U-Net architecture consists of an encoding sub-network and a decoding sub-network. They enable the SeisFault-Net to train models in an end-to-end manner. The residual neural network can suppress the gradient dispersion of deep network, and effectively improve the training efficiency of the SeisFault-Net. After trained, the SeisFault-Net can perform fault interpretation based on seismic data without setting any parameters. This avoids the empirical error and uncertainties caused by parameters artificially set in conventional methods. Applications to raw data have proved that the SeisFault-Net me-thod can effectively and accurately detect fault loca-tions, and the faults have good vertical continuity and clear outlines. The detailed information of faults interpretated by the SeisFault-Net method is more abundant and accurate than the coherent algorithm. And the calculating efficiency of the SeisFault-Net method is very high in seismic fault interpretation.
Keywords: fault interpreation    deep learning    resi-dual neural network    U-Net    seismic data interpretation    
0 引言

地震数据解释过程中一项重要的任务是断层解释,因为断层可以成为油气运移的主要通道,也可以侧向封堵而形成断块油气藏。因此,准确识别断层对油气的勘探开发具有重要意义。断层在地震剖面上表现为同相轴低连续或高度不连续性。传统的断层解释多由人机交互方式完成,效率低,人为因素可能导致结果不确定性大,加大了油气勘探开发的成本和风险。

为了提高断层解释的效率,研究人员提出了很多方法。Bahorich等[1]提出互相关方法,对每道地震数据及其邻域的地震道进行归一化互相关以识别断层。然而,该方法不适用于低信噪比的地震数据。Marfurt等[2]和Hale[3]将相似性算法用于检测地震数据中的不连续性。Marfurt等[4]提出特征结构相干算法,用于构造具有特征分解矩阵的协方差矩阵和特征分解矩阵。这些算法对断层识别的精度均很高,但对于低信噪比条件下的地震数据则需要更多的计算时间。Pedersen等[5]首次将蚂蚁追踪方法用于提取小断层,且合并形成较大的断层面。Aurnhammer等[6]提出遗传算法用于自动化地震图像中的断层识别。印兴耀等[7]提出使用离心窗倾角扫描的曲率属性识别小断层。陈雷等[8]利用相似性传播聚类与主成分分析进行断层识别。赵凤全等[9]将构造导向滤波用于断裂识别,取得了较好的应用效果。Wu等[10]提出3D图像处理方法,自动计算与断层相邻的图像样本的断层面和倾斜滑动矢量。上述方法可用于断层识别,但基于经验的参数设置可影响断层识别结果的准确性。

近年来,深度学习在多个领域取得了显著的成就,特别是在分类、特征提取等方面具有突出的优势。在地震勘探领域,研究人员探索了深度学习方法在地震数据处理和解释中的应用。Srisutthiyakorn[11]将深度学习方法用于预测渗透率。Lewis等[12]应用深度学习方法,利用FWI反演模型的地震图像建立地质模型。Huang等[13]应用深度学习方法,利用地震属性分析地质特征。Serfaty等[14]通过主成分分析和局部角度域深度学习分离地震波场。Chang等[15-16]将生成对抗网络用于地震数据插值。Wang等[17]将深度残余网络用于地震数据插值。王钰清等[18]利用卷积神经网络(Convolutional Neural Network,CNN)压制地震随机噪声,通过数据增广的方式生成标签数据,以提高去噪模型的泛化能力。张玉玺等[19]利用深度学习进行多属性盐丘自动识别。李海山等[20]利用深度残差网络压制地震数据噪声。Xiong等[21]和Wu等[22]应用深度学习方法进行断层智能化识别,取得了初步效果。Xiong等[21]使用相干体作为断层标签数据,训练深度学习模型,预测结果与应用相干体技术预测断层结果相似,精度并未明显提高。Wu等[22]使用正演合成数据制作标签数据集,用于深度学习模型训练,在构造简单的实际工区数据断层识别中取得了较好效果,但在复杂构造工区数据的断层识别中存在明显假象。Chang等[23]提出利用U-Net网络进行断层识别,利用合成数据和实际数据一起训练深度学习模型,试图解决智能断层识别网络的泛化能力弱的问题。

为了提高地震数据断层识别的精度和效率,本文结合U-Net架构和深度残差网络(ResNet),构建智能地震数据断层识别网络(SeisFault-Net)。U-Net架构适用于数据特征检测和提取,ResNet网络结构有助于训练深度深层网络,因此利用两种网络架构的优势可有效提高深层网络训练收敛速度和断层识别精度。同时,SeisFault-Net利用GPU进行数据处理,因此可大幅提高断层检测效率。

1 方法原理 1.1 ResNet结构

ResNet是一种用于训练更深层网络的框架[20],通过在卷积神经网络块(图 1a)上添加一个跳跃连接(Skip Connection),将输入向量和卷积层的输出融合在一起构建残差学习块(图 1b)。利用图 1a所示的卷积神经网络块构建深层网络时,层数的增加会导致深层存在网络梯度消失或者爆炸现象,使深度网络模型难以有效训练,出现模型退化问题。He等[24]提出的残差神经网络块(图 1b)用于解决深层网络模型退化问题,该网络结构基于恒等映射的思想,即每堆叠一层网络,堆叠前的输入与堆叠后的输出相同,使新增网络层只需学习一个残差函数即可。

图 1 卷积神经网络块(a)和残差神经网络块(b) x为输入地震数据;F(x)为残差函数;ReLU为激活函数

深度残差网络由堆叠的多个残差学习块组成,每个残差学习块的残差函数为

$ \begin{array}{c} I\left(\boldsymbol{x}_{l}\right)=\boldsymbol{x}_{l}+F\left(\boldsymbol{x}_{l}, \theta_{l}\right) \end{array} $ (1)
$ \boldsymbol{x}_{l+1}=f\left[I\left(\boldsymbol{x}_{l}\right)\right] $ (2)

式中:$ I\left(\boldsymbol{x}_{l}\right)$为同等映射函数;xlxl+1表示第l个残差学习块的输入和输出,lNN为网络堆叠层数;F是残差函数;θl表示第l个残差学习块的卷积核;激活函数f使用的是整流线性单元函数ReLU [25]

同等映射函数允许将较深层中的梯度值快速传递到较浅层,这有效解决了网络训练过程中由于链式求导法则导致的梯度爆炸或梯度消失的问题,从而提高网络训练过程的稳定性和效率。

1.2 SeisFault-Net结构

U-Net架构最初由Ronneberger等[26]提出并用于生物医学图像的目标分割。U-Net架构具有多尺度、多层次的特征,能有效识别、定位目标的边界。U-Net由编码子网络和解码子网络构成。编码子网络对数据进行编码学习,池化层和上采样层将数据流放大、缩小到不同尺度,实现对数据降维或升维计算。解码子网络学习在不同尺度下的数据流中的特性信息,使U-Net架构可以进行端到端训练,实现数据的语义分割。在编码子网络和解码子网络之间利用跳跃连接将多层信息进行融合,可提高目标特征识别的精度。

鉴于断层识别需要准确刻画边界和位置,本文将U-Net网络和ResNet结合构建用于地震数据断层识别的SeisFault-Net。SeisFault-Net将ResNet训练深层网络的优势与U-Net架构的特征分割相结合,实现对地震数据断层的高精度识别。

SeisFault-Net如图 2所示,设计为对称结构,由三类层构成,包括ResNet块、池化层(Pooling)和上采样层(Up-Sampling)。其中,ResNet块由3个卷积层和1个跳跃连接构成,第一、第二、第三层卷积核分别为1×1、3×3、1×1,共11个块。利用两个1×1的卷积层代替3×3的卷积层降低整个网络模型的学习参数量,在不影响残差学习块的性能的同时,对原始3×3卷积核网络的冗余参数进行压缩。池化层为步长2×2的最大池化,共4层,用于地震数据降维,提高后续的卷积层对地震数据的感受野。上采样层是卷积核为2×2、步长为1×1的反卷积层,共4层。使用反卷积层可以从映射层中学习重建的特征,使上采样的数据能够更准确地重建断层特征。

图 2 SeisFault-Net网络 (a)地震数据;(b) SeisFault-Net网络结构;(c)标签数据

SeisFault-Net可表示为

$ \hat{y}_{i}=\text { SeisFault-Net }\left(\boldsymbol{x}_{i} ; \boldsymbol{W}, b\right) $ (3)

式中:$ \hat{y}_{i}$表示SeisFault-Net预测的第i个断层数据;W为卷积核;b为偏置。

对SeisFault-Net网络模型进行优化的损失函数为二分类的交叉熵函数,即

$ \begin{aligned} \ln \left(y_{i} \mid \hat{y}_{i}\right)=&-\left[y_{i} \cdot \ln \left(\hat{y}_{i}\right)+\right.\\ &\left.\left(1-y_{i}\right) \cdot \ln \left(1-\hat{y}_{i}\right)\right] \end{aligned} $ (4)

式中yi为第i个断层标签数据。获得损失函数后,利用梯度下降求取$ \hat{y}_{i}$对(W, b)的导数,通过反向传播算法更新(WW, bb)(ΔW、Δb表示每一次训练更新的量),使损失函数最小化,达到模型训练的目的。

2 数据试算 2.1 训练数据集

为了提高SeisFault-Net的泛化能力,本文共选取7个工区的地震数据(图 3a),将地震数据分割为256×256的数据块。其中,5个工区的数据用于制作训练数据集和验证数据集,共有800个地震数据块;2个工区的地震数据用于测试数据集,共有150个地震数据块。实际地震数据集中包含线性的和非线性的断层类型,标注断层可获得标签数据集(图 3b)。

图 3 训练数据集 (a)地震数据;(b)标签数据
2.2 SeisFault-Net训练

首先,利用训练地震数据集训练SeisFault-Net,使用动态学习率,即

$ \left\{ {\begin{array}{*{20}{l}} {{I_k} = {I_0}}&{1 \le k \le 20}\\ {{I_k} = {I_0}{r^{\frac{g}{{{s_{k - 1}}}}}}}&{k > 20} \end{array}} \right. $ (5)

其中

$ s_{k-1}=20 \operatorname{int}\left(\frac{k-1}{20}\right) $

式中:I0=0.01为初始学习率;r为衰减系数;s为衰减步长;g为全局步长;k为训练次数;int(·)表示取整。随着训练次数的增加,学习率降低,这有助于提高模型收敛速度,降低模型训练时间。图 4a为训练实际数据集的损失函数曲线。由图可见,经过100次训练后,损失函数值基本稳定。图 4b为预测断层准确率曲线,验证集的断层预测准确率在80次训练后基本稳定在0.92附近。

图 4 SeisFault-Net训练曲线 (a)损失函数曲线;(b)断层预测准确率
2.3 SeisFault-Net模型验证

利用验证数据集和测试数据集测试训练后Seis-Fault-Net网络模型的断层预测性能。图 5为验证集数据断层识别结果。由图可见,经过训练后的SeisFault-Net在验证数据集上断层预测准确,与标签数据一致。

图 5 SeisFault-Net验证数据集断层识别结果 (a)地震数据;(b)断层标签数据;(c)SeisFault-Net预测断层

测试数据集断层识别结果如图 6所示。训练后的SeisFault-Net在测试数据上断层预测效果较好,准确刻画了断层位置和边界。这说明SeisFault-Net具备应用于实际断层识别的泛化能力。

图 6 SeisFault-Net测试数据集断层识别结果 (a)地震数据;(b)SeisFault-Net预测断层

利用测试数据进行断层识别准确率和效率测试。图 7为SeisFault-Net方法和相干体技术断层识别效果对比,可以看到SeisFault-Net方法明显优于相干算法,且计算效率为相干算法的18倍。

图 7 SeisFault-Net与相干算法断层识别效果对比 (a)地震数据;(b)相干算法;(c)SeisFault-Net方法
3 应用效果 3.1 泛化性能

使用实际数据体测试SeisFault-Net断层识别能力。图 8为DJ工区地震剖面识别断层的投影显示,通过训练后的SeisFault-Net无需设置任何参数即可进行线性和非线性断层的识别。由图 8可见,SeisFault-Net识别的断层位置与地震同相轴不连续处基本重合,且断层垂向连续性好;黄色虚线框内识别的断层清晰,位置准确;蓝色框内由于地震资料分辨率较低,断层的垂向连续性受到影响,但SeisFault-Net也可以精确识别断层的位置和方向,这表明经过训练后的SeisFault-Net的泛化能力较强。

图 8 实际数据体预测断层的投影显示 黑色投影线为SeisFault-Net方法识别的断层
3.2 断层识别

利用DQ工区三维地震数据体测试SeisFault-Net,并且与相干算法进行比较。三维地震数据体大小为1.02GB,使用机器型号为Nvidia P100 GPU。SeisFault-Net识别断层耗时5.53s,相干算法识别断层耗时107.00s,即SeisFault-Net计算效率是相干算法的19.35倍。

图 9为DQ工区三维数据体的等时切片。由图可见,SeisFault-Net识别的断层位置准确,横向展布清晰,断层连续性好(如图中黄色箭头处);与相干算法识别的断层(图 9b)相比,SeisFault-Net识别的断层(图 9c)细节特征更丰富,断层展布清晰(如图中黄色方框内)。图 10为地震剖面对比结果,可见相干算法识别断层连续性差,存在水平干扰现象,而SeisFault-Net识别的断层位置准则,不存在水平干扰。

图 9 实际数据等时切片 (a)地震数据等时切片;(b)相干算法识别断层切片;(c)SeisFault-Net识别断层切片

图 10 地震剖面不同方法断层识别结果对比 (a)地震数据;(b)相干算法;(c)SeisFault-Net方法
4 结论

本文提出基于深度残差网络与U-Net架构联合识别地震资料断层的方法,能有效识别实际数据中的断层。SeisFault-Net网络利用了U-Net架构,通过编码子网络和解码子网络在高维空间有效表征地震数据中的断层特征,利用跳跃连接将多层的特征信息进行融合,提高了断层识别准确性。同时,SeisFault-Net利用ResNet层,提高了网络训练效率。相较于相干算法,本文方法能更有效、准确识别断层位置,识别的断层的垂向连续性更好。SeisFault-Net用于地震断层识别的计算效率高,能极大缩短地震数据的断层解释时间,降低断层解释结果的人为不确定性。

参考文献
[1]
Bahorich M, Farmer S. 3-D seismic discontinuity for faults and stratigraphic features:The coherence cube[J]. The Leading Edge, 1995, 14(10): 1053-1058. DOI:10.1190/1.1437077
[2]
Marfurt K J, Kirlin R L, Farmer S L, et al. 3-D seismic attributes using a semblance-based coherency algorithm[J]. Geophysics, 1998, 63(4): 1150-1165. DOI:10.1190/1.1444415
[3]
Hale D. Methods to compute fault images, extract fault surfaces, and estimate fault throws from 3D seismic images[J]. Geophysics, 2013, 78(2): O33-O43. DOI:10.1190/geo2012-0331.1
[4]
Marfurt K J, Sudhaker V, Gersztenkorn A, et al. Coherency calculations in the presence of structural dip[J]. Geophysics, 1999, 64(1): 104-111. DOI:10.1190/1.1444508
[5]
Pedersen S I, Randen T, Sonneland L, et al. Automatic fault extraction using artificial ants[C]. SEG Technical Program Expanded Abstracts, 2002, 21: 512-515.
[6]
Aurnhammer M, Tonnies K D. A genetic algorithm for automated horizon correlation across faults in seismic images[J]. IEEE Transactions on Evolutionary Computation, 2005, 9(2): 201-210. DOI:10.1109/TEVC.2004.841307
[7]
印兴耀, 高京华, 宗兆云. 基于离心窗倾角扫描的曲率属性提取[J]. 地球物理学报, 2014, 57(10): 3411-3421.
YIN Xingyao, GAO Jinghua, ZONG Zhaoyun. Curvature attribute based on dip scan with eccentric window[J]. Chinese Journal of Geophysics, 2014, 57(10): 3411-3421. DOI:10.6038/cjg20141027
[8]
陈雷, 肖创柏, 禹晶, 等. 基于相似性传播聚类与主成分分析的断层识别方法[J]. 石油地球物理勘探, 2017, 52(4): 826-833.
CHEN Lei, XIAO Chuangbai, YU Jing, et al. Fault recognition based on affinity propagation clustering and principal component analysis[J]. Oil Geophysical Prospecting, 2017, 52(4): 826-833.
[9]
赵凤全, 崔德育, 康婷婷, 等. 构造导向滤波技术在断裂识别中的应用[J]. 石油地球物理勘探, 2018, 53(增刊1): 214-218.
ZHAO Fengquan, CUI Deyu, KANG Tingting, et al. Fault identification with structure-oriented filtering[J]. Oil Geophysical Prospecting, 2018, 53(S1): 214-218.
[10]
Wu X M, Hale D. 3D seismic image processing for faults[J]. Geophysics, 2016, 81(2): IM1-IM11. DOI:10.1190/geo2015-0380.1
[11]
Srisutthiyakorn N.Deep-learning methods for predicting permeability from 2D/3D binary-segmented images[C]. SEG Technical Program Expanded Abstracts, 2016, 35: 3042-3046.
[12]
Lewis W, Vigh D.Deep learning prior models from seismic images for full-waveform inversion[C]. SEG Technical Program Expanded Abstracts, 2017, 36: 1512-1517.
[13]
Huang L, Dong X, Clee T E. A scalable deep learning platform for identifying geologic features from seismic attributes[J]. The Leading Edge, 2017, 36(3): 249-256. DOI:10.1190/tle36030249.1
[14]
Serfaty Y, Itan L, Chase D, et al.Wavefield separation via principle component analysis and deep learning in the local angle domain[C]. SEG Technical Program Expanded Abstracts, 2017, 36: 991-995.
[15]
Chang D K, Yang W Y, Yong X S, et al. Generative adversarial networks for seismic data interpolation[C]. SEG 2018 Workshop, SEG Maximizing Asset Value Through Artificial Intelligence and Machine Learning, Beijing, China, 2018.
[16]
Chang D K, Yang W Y, Yong X S, et al. Seismic data interpolation using dual-domain conditional generative adversarial networks[J]. IEEE Geoscience and Remote Sensing Letters, 2020, 99: 1-5.
[17]
Wang B F, Zhang N, Lu W K, et al. Deep-learning-based seismic data interpolation:A preliminary result[J]. Geophysics, 2019, 84(1): V11-V20. DOI:10.1190/geo2017-0495.1
[18]
王钰清, 陆文凯, 刘金林, 等. 基于数据增广和CNN的地震随机噪声压制[J]. 地球物理学报, 2019, 62(1): 421-433.
WANG Yuqing, LU Wenkai, LIU Jinlin, et al. Random seismic noise attenuation based on data augmentation and CNN[J]. Chinese Journal of Geophysics, 2019, 62(1): 421-433.
[19]
张玉玺, 刘洋, 张浩然, 等. 基于深度学习的多属性盐丘自动识别方法[J]. 石油地球物理勘探, 2020, 55(3): 475-483.
ZHANG Yuxi, LIU Yang, ZHANG Haoran, et al. Multi-attribute automatic interpretation of salt domes based on deep learning[J]. Oil Geophysical Prospecting, 2020, 55(3): 475-483.
[20]
李海山, 陈德武, 吴杰, 等. 叠前随机噪声深度残差网络压制方法[J]. 石油地球物理勘探, 2020, 55(3): 493-503.
LI Haishan, CHEN Dewu, WU Jie, et al. Pre-stack random noise suppression with deep residual network[J]. Oil Geophysical Prospecting, 2020, 55(3): 493-503.
[21]
Xiong W, Ji X, Ma Y, et al. Seismic fault detection with convolutional neural network[J]. Geophysics, 2018, 83(5): O97-O103. DOI:10.1190/geo2017-0666.1
[22]
Wu X M, Shi Y Z, Fomel S, et al.Convolutional neural networks for fault interpretation in seismic images[C]. SEG Technical Program Expanded Abstracts, 2018, 37: 1946-1950.
[23]
Chang D, Yong X, Yang W, et al.U_Net & residual neural networks for seismic fault interpretation[C]. Extended Abstracts of 81st EAGE Conference and Exhibition, 2019, 1-5.
[24]
He K M, Zhang X, Ren S, et al.Identity mappings in deep residual networks[C]. European Conference on Computer Vision, 2016, 630-645.
[25]
Nair V, Hinton G E.Rectified linear units improve restricted Boltzmann machines[C]. International Conference on Machine Learning, 2010, 807-814.
[26]
Ronneberger O, Fischer P, Brox T. U-Net: convolutional networks for biomedical image segmentation[C]. Medical Image Computing and Computer-assisted Intervention-MICCAI 2015, Springer, 2015, 234-241.