收藏本站
《同济大学》 2004年
收藏 | 手机打开
二维码
手机客户端打开本文

地震反射走时层析理论与应用研究

成谷  
【摘要】: 速度建模是地震勘探中的一项重要技术。其中偏移速度分析和层析成像是两种很重要的手段。从理论上看,速度建模的精度应以层析成像为高,但层析成像在实践过程中受诸多因素的制约,其难点不在于理论方法上,而在于具体实现的技巧上。 本文对两种速度建模方法进行了系统的研究,并重点研究了利用反射地震资料进行走时层析成像速度建模的方法。考虑到层析成像在地震中应用的很多理论和实际上的困难,更由于反射层析中速度和反射深度的耦合性,本文确立了用偏移速度分析的结果为层析成像提供初始速度和反射深度的思想,使得层析成像能够具有一个好的起点。 本文对偏移速度分析的发展现状进行了简单的介绍,详细阐述了层析成像中存在的问题和国内外的发展现状,并针对这两种速度建模方法进行了详细的比较。本文对偏移速度分析的要求不高,偏移速度分析的结果只需给出常速模型即可,本文更对偏移速度分析提供的初始反射深度感兴趣。鉴于速度分析的迭代型和交互性,本文选取了速度较快的局部区域的Kirchhoff偏移方法和沿层进行偏移速度分析的策略。 关于地震层析成像,本文主要介绍了几个方面的内容:层析方程组的数值解法、正则化和归一化。层析方程组的求解,经过了矩阵直接求逆和迭代解两个阶段。矩阵直接求逆的解法,只对小型的矩阵适用。迭代解法主要分为ART类方法和投影类方法。ART类方法是一种行处理方法,是针对早期计算机能力的限制而采取的方法。投影类主要有Lanczos方法、共轭梯度法(CG)、最小二乘共轭梯度法(LSCG)以及最小二乘QR分解法(LSQR)。是目前最经常使用的一类方法。在地震层析中利用正则化对欠定分量和零空间分量进行约束是非常重要的,正则化对欠定分量和零空间分量的约束主要体现为最简单解的思想:包括最平坦解、最光滑解和最小长度解思想,最平坦解思想的出发点是使解的分量间只发生平缓的变化,最光滑解的思想是解的分量间的变化是光滑的,最小长度解思想是使解矢量相对于某个先验矢量发生最小的变化。本文将正则化的具体作用方式分为加法型和乘法型两种。加法型是指将正则化方程组补在原方程组的下面的处理方式,乘法型是指将正则化方程组与原方程组相乘的处理方式。在对正则化的作用范围和作用方式进行介绍后,本文又对正则化的矩阵表述形式、压缩存储进行了阐述,并针对简单模型对导数型正则化的效果进行了分析。归一化是地震层析成像中另一个非常重要的问题。反演计算中,数据d和模型参数m的物理量纲不同,而在多参数类型同时反演时,模型参数和数据中的元本身就带有不同的物理量纲。此时Frechet微商矩阵的列向量的模相差极大,导致问题的性状很坏。因此在进行数据反演计算之前,应该把所有的模型参数和数据进行归一化化为无量纲的量。 针对反射层析中速度和深度的耦合现象,本文综合多人的研究结果得出结论如下:速度深度的耦合性在本质上是由短偏移距引起的,大偏移距的存在能够缓解这种耦合。速度深度的耦合程度受几方面因素的影响:(1)偏移距与深度的比值,偏移距与深度的比值越大,速度深度的耦合性越弱;(2)拾取误差项,拾取误差越小,速度深度的耦合性越弱;(3)界面倾角,倾角越小,速度深度的耦合性越弱:(4)扰动波长与特征波长的比值,当扰动波长大于特征波长时慢度扰动占主要地位,当扰动波长小于特征波长时,深度扰动占主要地位:当扰动波长与特征波长相等时,速度和反射深度发生耦合。特征波长取决于异常体的厚度和异常体与反射面的距离。此外本文还分析了走时误差对反射层析中利用走时扰动分辨薄层内速度扰动的稳定性的影响,并得出结论:当层厚较薄时,由于走时信息中误差的影响,利用走时扰动分辨薄层内的速度扰动是非常不稳定的。层内速度越高、偏移距越小,不稳定性越强。本文针对反射层析的横向和垂向分辨能力也进行了分析,分析证明由于反射地震勘探中以近垂直方向射线为主,因而反射层析的横向分辨能力要优于垂向分辨能力。 本文针对矩形网格参数化的形式建立了层析反演体系,在正演模拟中选用了Langan射线追踪方法,这种方法主要有三个优点:(1)网格内速度的三参数(对二维情况)表述方式在本质上照顾到了射线理论对速度平滑性的要求,方法本身有一种内在的平滑效应;(2)Langan法射线追踪具有完美的解析表达形式,且能够保证精度和效率;(3)在模型矢量选为慢度矢量差,数据矢量选为走时差的情况下,射线在网格内的长度就是反映数据变化对参数变化灵敏度的Frechet微商。考虑到均一矩形网格参数化导致的巨大的矩阵维数给计算机存储量和计算量造成的压力,在对矩形网格参数化层析反演中,本文采用了行索引的压缩存储方法,并以压缩存储的方式补入了一阶导数型正则化和二阶导数型正则化矩阵,最后针对补入正则化矩阵的压缩矩阵采用了LSQR方法迭代求解,大大节省了存储量、计算量和计算时间。 尽管在目前的地震层析成像反演中,采用最多的是矩形网格参数化的方式,但矩形网格参数化方式在很多方面不容乐观:在矩形网格参数化形式下,对地下介质的精细描述必然导致巨大的矩阵维数(各种正则化约束加入后尤甚),方程组的大型稀疏特性使得在求解过程中耗费了大量的计算机存储空间和计算时间。尽管本文针对方程组的稀疏特性设计了一种行索引的压缩存储方法以及针对压缩存储矩阵的迭代解法,大大节省了存储量、计算量和计算时间,但考虑到层析速度建模的迭代性和交互性,存储量和计算时间仍然是矩形网格参数化层析反演应用的一个限制因素。除此之外,矩形网格的均一剖分形式灵活性差,不论在构造简单区域和构造复杂区域都不能区别对待,增大了未知参数的个数,不仅在正演模拟时增加了计算量和计算时间,而且在层析反演时使问题的性态变得更坏,增加了层析反演的难度。矩形网格参数化在对界面的描述上也同样存在着缺陷,即速度模型和界面模型的不一致,由于矩形网格参数化只能用锯齿状的分布来表征界面上下的速度变化,在对界面的处理上增加了难度和引入了误差。针对矩形网格参数化方式在剖分上的不灵活给正反演方法带来的诸多不利因素,考虑到三角网格剖分的灵活性和本质上的正则化效应,本文进一步建立了可变三角网格参数化的层析反演体系,得到了很大的改善。三角网格参数化层析反演相对于矩形网格参数化层析反演无论在计算时间、计算量,还是在存储量上都占有非常大的优势,经验证两种参数化方式得到的层析结果具有可比性。
【学位授予单位】:同济大学
【学位级别】:博士
【学位授予年份】:2004
【分类号】:P631.4

【引证文献】
中国期刊全文数据库 前3条
1 孟宪海;金颖;李吉刚;谭文磊;杨钦;;基于三角域快速行进法的地震波走时计算[J];软件;2011年11期
2 秦宁;李振春;杨晓东;;基于角道集的井约束层析速度反演[J];石油地球物理勘探;2011年05期
3 成谷;张宝金;;基于参数模型两种剖分方式的射线追踪比较[J];中山大学学报(自然科学版);2008年05期
中国重要会议论文全文数据库 前1条
1 黄光南;刘洋;;层析成像数据和射线不均匀覆盖问题的解决方法[A];中国地球物理学会第二十七届年会论文集[C];2011年
中国博士学位论文全文数据库 前1条
1 李军;复杂地表模型菲涅尔体层析反演静校正方法研究[D];成都理工大学;2011年
中国硕士学位论文全文数据库 前3条
1 孟庆利;基于大炮初至的层析静校正算法研究与实现[D];中国石油大学;2011年
2 李国磊;基于大炮初至数据的层析成像算法研究[D];中国石油大学;2010年
3 宋立芳;金属矿地震勘探成像方法技术应用研究[D];长安大学;2012年
【共引文献】
中国期刊全文数据库 前10条
1 郜延红,周云轩,刘万崧;地球物理位场可视化建模初步探讨[J];长春科技大学学报;2000年02期
2 孙建国;何洋;;基于波前构建的射线追踪:一种Java实现[J];吉林大学学报(地球科学版);2007年04期
3 韩复兴;孙建国;杨昊;;基于二维三次卷积插值算法的波前构建射线追踪[J];吉林大学学报(地球科学版);2008年02期
4 张营;冯晅;刘财;赵旭;兰慧田;张瑾;卢勇旭;杨庆节;;转换点位置对转换波速度分析精度的影响[J];吉林大学学报(地球科学版);2010年S1期
5 孔雪;李振春;叶月明;黄建平;徐秀刚;;多尺度地震资料联合成像[J];吉林大学学报(地球科学版);2011年S1期
6 曹小林,洪学海,曹俊兴;面波波形反演中的模拟退火法[J];成都理工学院学报;2000年03期
7 王志,贺振华,薛彬,黄德济,王熙明,周亚同;广角反射研究及高斯射线正演模拟分析[J];成都理工学院学报;2002年01期
8 吴建平,冯锐;井间电磁波CT中的非块体分段反演算法[J];CT理论与应用研究;1993年03期
9 张兆田,邱佩璋;一种优化的图象重建外插算法[J];CT理论与应用研究;1993年04期
10 闻则刚,何璟琳;用于势场反演的特殊解法[J];CT理论与应用研究;2001年02期
中国重要会议论文全文数据库 前10条
1 刘福平;杨长春;;孔隙度敏感系数的数值计算[A];中国科学院地质与地球物理研究所二○○三学术论文汇编·第四卷(油气资源)[C];2003年
2 耿小兵;张仲山;;核钻地弹爆炸地下爆震破坏范围确定方法探讨[A];全国第六届核仪器及其应用学术会议论文集[C];2007年
3 龙士国;马天朗;;稀土超磁致伸缩换能器的研制及其在桥梁检测中的应用[A];第六届中国功能材料及其应用学术会议论文集(3)[C];2007年
4 余易豪;舒玉峰;丁国辉;李铭铭;;震波CT探测技术及其应用[A];矿山地质灾害成灾机理与防治技术研究与应用[C];2009年
5 李满树;;同时反演复杂地区速度与界面的速度建模方法[A];油气地球物理实用新技术——中国石化石油勘探开发研究院南京石油物探研究所2004年学术交流会论文集[C];2004年
6 龙士国;马天朗;;超磁震源声波CT仪的研制及其在工程无损评价中的应用[A];中国仪器仪表学会第九届青年学术会议论文集[C];2007年
7 罗水余;申秋学;;井间地震射线层析成像的分辨率讨论[A];勘探地球物理2005学术交流会论文集[C];2005年
8 夏建军;宋强弓;尹业才;宋锡熊;;复杂地表条件下的地震波场外推法模拟分析[A];中国地球物理第二十一届年会论文集[C];2005年
9 李满树;;一种利用多种地震波走时建立表层速度模型的方法[A];中国地球物理学会第22届年会论文集[C];2006年
10 徐兆涛;;单程波动方程叠前深度偏移并行计算[A];中国地球物理学会第二十三届年会论文集[C];2007年
中国博士学位论文全文数据库 前10条
1 刘俊峰;海拉尔盆地复杂储层地震识别方法研究[D];中国地质大学(北京);2010年
2 鲁晶津;地球电磁三维数值模拟的多重网格方法及其应用研究[D];中国科学技术大学;2010年
3 曹令敏;地震层析成像方法的研究及其在金属矿复杂模型模拟中的应用[D];中国地质大学(北京);2011年
4 王连山;廊固凹陷大兴砾岩体天然气藏评价预测研究[D];中国地质大学(北京);2011年
5 汪自军;基于卫星临边辐射的大气痕量气体含量反演研究[D];吉林大学;2011年
6 张向宇;工业炉温度场可视化与辐射特性参数解耦重建研究[D];华中科技大学;2011年
7 张风雪;有限频体波走时层析成像及其在华北地区的应用[D];中国地震局地球物理研究所;2011年
8 彭远新;能源消费与经济发展的多尺度分析和反演[D];南京师范大学;2011年
9 韩翀;海相碳酸盐岩地震储层预测技术及应用[D];成都理工大学;2011年
10 朱四新;琼东南盆地深水区低幅构造速度研究与PSDM应用[D];成都理工大学;2011年
中国硕士学位论文全文数据库 前10条
1 田禹;煤岩双相介质弹性波场数值模拟[D];山东科技大学;2010年
2 李要;稀土超磁致伸缩换能器的动态特性及其应用[D];湘潭大学;2010年
3 宋小超;南海均衡重力异常研究[D];中国地质大学(北京);2011年
4 郭伟;隧道掘进超前地质勘探的研究与仿真[D];沈阳工业大学;2011年
5 赵旭;转换波多次覆盖道集抽取研究[D];吉林大学;2011年
6 杨文;利用背景噪声方法研究云南强震前后的波速变化[D];中国地震局地震预测研究所;2011年
7 刘乾;茶汗通古盆地侏罗纪地层分布特征研究[D];新疆大学;2011年
8 冯红刚;隐患空区群下采场结构参数优化研究[D];中南大学;2011年
9 王冉;复杂地电模型的2.5-D直流电阻率边界单元数值模拟[D];中南大学;2011年
10 杨婷;华北地区地壳上地幔S波速度结构研究[D];中国地震局地球物理研究所;2011年
【同被引文献】
中国期刊全文数据库 前10条
1 李红;;复杂地区静校正技术应用研究[J];安徽地质;2007年01期
2 彭直兴,张芬,周熙襄,钟本善;地震波旅行时非线性/线性联合插值法[J];成都理工大学学报(自然科学版);2005年03期
3 王金峰;罗省贤;李录明;;改进的神经网络并行算法及其在地震初至拾取中的应用[J];成都理工大学学报(自然科学版);2007年03期
4 徐明才,高景华,荣立新,刘建勋;从金属矿地震方法的试验效果探讨其应用前景[J];中国地质;2004年01期
5 皮金云;段云卿;康平;;综合折射静校正方法在柴达木盆地的应用[J];地球科学(中国地质大学学报);2007年04期
6 朱丹;波前快速推进法三维走时计算技术[J];大庆石油地质与开发;2001年05期
7 高磊;潘树林;周熙襄;史艳华;;利用转换波初至进行转换波剩余静校正[J];大庆石油地质与开发;2008年01期
8 高星;地震层析成像研究的回顾与展望[J];地球物理学进展;2000年04期
9 熊熊,滕吉文,许厚泽;地幔热柱动力学研究的新进展[J];地球物理学进展;2001年01期
10 傅容珊,郑勇,常筱华,黄建华;地震层析成像板块构造及地幔演化动力学[J];地球物理学进展;2001年04期
中国重要会议论文全文数据库 前2条
1 何晓冬;夏忠谋;张国保;;模型法静校正在西部地区的应用[A];中国地球物理·2009[C];2009年
2 张正峰;韩立国;;变换坐标系下的转换波波动方程延拓静校正[A];中国地球物理·2009[C];2009年
中国博士学位论文全文数据库 前8条
1 巩向博;金属矿地震高精度成像与数据处理方法研究[D];吉林大学;2011年
2 王克斌;复杂地表条件下初至折射波静校正方法研究[D];成都理工大学;2004年
3 邹强;山地静校正若干问题研究[D];成都理工大学;2004年
4 于富文;中国西部地区静校正方法对比研究[D];中国地质大学(北京);2006年
5 马德堂;弹性波场数值模拟及井间地震初至波旅行时层析成像[D];长安大学;2006年
6 高磊;地震初至波自动拾取和折射波表层调查方法研究[D];成都理工大学;2008年
7 刘素芹;基于应用网格环境的复杂地表波动方程基准面静校正研究[D];中国石油大学;2008年
8 王兆国;东北地区地震综合信息联合成像与震源机制反演[D];吉林大学;2010年
中国硕士学位论文全文数据库 前10条
1 朱龙生;多方位角逆VSP层析成像[D];长安大学;2003年
2 赵德亨;折射波广义互换法在城市地质调查中的应用[D];吉林大学;2005年
3 彭直兴;地震波初至旅行时射线追踪正演方法研究[D];成都理工大学;2005年
4 刘洪雷;复杂地区的折射波静校正应用研究[D];中国地质大学(北京);2006年
5 秦亚玲;沙漠地区地震资料处理技术研究[D];中国地质大学(北京);2006年
6 尹纪超;基于网格的球面波波前射线追踪方法研究[D];成都理工大学;2006年
7 张丽;复杂地区静校正技术研究[D];长安大学;2006年
8 姜国博;初至波表层模型层析正反演及层析静校正研究[D];中国地质大学(北京);2007年
9 喻维秋;初至波表层模型层析反演及其应用[D];成都理工大学;2007年
10 田连玉;基于模拟退火和遗传算法的剩余静校正方法研究[D];中国石油大学;2007年
【二级引证文献】
中国期刊全文数据库 前5条
1 秦宁;李振春;杨晓东;贾正虹;岳玉波;;叠前多级优化联合偏移速度建模[J];地球物理学进展;2013年01期
2 秦宁;李振春;杨晓东;张凯;王俊;;共散射点道集与角道集串级优化叠前偏移速度分析[J];吉林大学学报(地球科学版);2013年02期
3 刘润泽;于师建;田清伟;陆二男;;基于波前最小走时单元的三角网射线追踪全局算法[J];地球物理学进展;2013年02期
4 秦宁;李振春;杨晓东;周卿;陈富强;桑运云;;自动拾取的成像空间域走时层析速度反演[J];石油地球物理勘探;2012年03期
5 秦宁;李振春;杨晓东;岳玉波;张凯;张晴;;高斯波束角道集的旅行时层析速度分析[J];石油地球物理勘探;2013年03期
中国重要会议论文全文数据库 前1条
1 秦宁;李振春;;基于角度域共成像点道集的走时层析速度反演[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
中国硕士学位论文全文数据库 前2条
1 孟庆利;基于大炮初至的层析静校正算法研究与实现[D];中国石油大学;2011年
2 高杨;三维高密点地震数据预处理方法研究[D];南京理工大学;2013年
【相似文献】
中国期刊全文数据库 前10条
1 张东;童海杰;杨艳;秦前清;徐凌;;基于波前扩展的线性走时插值射线追踪算法[J];武汉大学学报(理学版);2010年06期
2 王云安;艾东海;;基于电磁波CT的岩溶地质异常探测[J];人民黄河;2009年12期
3 黄靓,黄政宇;线性插值射线追踪的改进方法[J];湘潭大学自然科学学报;2002年04期
4 赵连锋,朱介寿,曹俊兴;并行化交错网格法地震层析成像[J];石油物探;2003年01期
5 杨慧珠,聂建新,谭桂华;反射波地震层析成像算法研究[J];辽宁工程技术大学学报;2004年02期
6 王祯祥,高尔根,徐果明,刘同庆;一种层状介质条件下射线的全路径迭代追踪法[J];煤田地质与勘探;2002年04期
7 ;Geophysics,No.1,1992论文文摘[J];勘探地球物理进展;1992年05期
8 许永忠,柴振友,宋端智;煤矿地应力异常层析成像方法的研究[J];中国煤田地质;2000年01期
9 周兵,朱介寿,S.A.Greenhalgh,C.Sinadinovski;一种新的地震射线层析成像计算方法[J];石油物探;1994年01期
10 林伯香,孙晶梅,刘清林;层析成像低速带速度反演和静校正方法[J];石油物探;2002年02期
中国重要会议论文全文数据库 前10条
1 周连庆;赵翠萍;陈章立;;利用地方震层析成像研究龙滩水库库区三维V_P和V_P/V_S层析成像[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
2 孙晶梅;李博;林伯香;徐颖;;复杂近地表速度分析及校正[A];油气地球物理实用新技术——中国石化石油勘探开发研究院南京石油物探研究所2005年学术交流会论文集[C];2005年
3 周连庆;赵翠萍;郑现;陈章立;;三峡库区三维Vp,Vp/Vs,Qp和Qs精细结构层析成像[A];中国地球物理学会第二十七届年会论文集[C];2011年
4 李玉喜;刘四新;;跨孔层析成像四种算法的研究[A];中国地球物理·2009[C];2009年
5 曾吕明;邢达;谷怀民;杨迪武;杨思华;向良忠;;基于多元线性阵列探测器的快速微波致热声层析成像[A];第十次中国生物物理学术大会论文摘要集[C];2006年
6 裴顺平;许忠淮;汪素云;;中国大陆Sn波速度层析成像[A];中国地球物理学会年刊2002——中国地球物理学会第十八届年会论文集[C];2002年
7 高丽丽;何琳;徐果明;朱良保;;层析成像中反演参数选取的影响[A];中国地球物理学会第二十届年会论文集[C];2004年
8 房立华;吴建平;王未来;明跃红;;华北地区噪声层析成像研究[A];中国地震学会成立三十年学术研讨会论文摘要集[C];2009年
9 张宝龙;徐义贤;陈伟;蒋礼;;基于噪声的面波层析成像与观测实验[A];中国地球物理学会第二十七届年会论文集[C];2011年
10 肖化;胡广莉;何惠玲;保宗悌;;基于图象分割的人工神经网络在两相流层析成像中应用研究[A];中国图象图形科学技术新进展——第九届全国图象图形科技大会论文集[C];1998年
中国重要报纸全文数据库 前10条
1 华为技术有限公司 卢鹏;射线追踪技术为网络规划导航[N];通信产业报;2006年
2 黄寰;层析成像技术揭示地球奥秘[N];中国矿业报;2004年
3 本报记者 孟凡君;资源勘探的得力助手[N];中国矿业报;2008年
4 王五平;声波CT为混凝土大坝做体检[N];中国水利报;2005年
5 本报特约记者 曹菲;杨文采与他的应用地球物理反演理论[N];地质勘查导报;2006年
6 记者 刘如春 通讯员 刘文祥;地震CT技术洞悉岩溶变化[N];地质勘查导报;2008年
7 刘洪宇;英特尔Larrabee渲染架构[N];中国计算机报;2008年
8 石岩;中天山地幔精细结构研究项目[N];科技日报;2007年
9 强声;美国研制能透视救灾的无线网络系统[N];人民公安报·消防周刊;2010年
10 ;探索地球科学服务人类社会[N];科技日报;2003年
中国博士学位论文全文数据库 前10条
1 成谷;地震反射走时层析理论与应用研究[D];同济大学;2004年
2 王浩全;超声阵列层析成像检测方法研究[D];中北大学;2011年
3 马辉;基于MapReduce的分布式地震射线追踪方法研究[D];中国地质大学(北京);2012年
4 袁志亮;井间声波电磁波层析成像技术应用研究与软件研发[D];中国地质大学(北京);2007年
5 郝奇;VTI介质速度和各向异性参数建模研究[D];吉林大学;2010年
6 黄靓;混凝土超声波层析成像的理论方法和试验研究[D];湖南大学;2008年
7 喻振华;工程井间地震波正演模拟及层析成像技术研究[D];中南大学;2008年
8 白志明;深地震测深剖面的层析成像研究及其应用[D];中国地震局地球物理研究所;2002年
9 王振宇;土木工程的层析成像与广义反演研究[D];浙江大学;2003年
10 王兆国;东北地区地震综合信息联合成像与震源机制反演[D];吉林大学;2010年
中国硕士学位论文全文数据库 前10条
1 李国磊;基于大炮初至数据的层析成像算法研究[D];中国石油大学;2010年
2 郭彦刚;井间地震层析成像技术及应用研究[D];成都理工大学;2012年
3 周小仙;利用地震反射和折射联合层析成像探测深部构造方法研究[D];中国地质大学(北京);2010年
4 张赛民;跨孔地震层析成像研究[D];中南大学;2003年
5 燕利芳;反射波层析成像的数值模拟及在西安地裂缝中的应用研究[D];长安大学;2011年
6 李志辉;CT技术在路基病害勘查治理中的研究及应用[D];西南交通大学;2005年
7 金溪;浅层地震初至波层析成像方法技术研究[D];长安大学;2010年
8 罗文歆;井间地震层析成像研究及其在隐伏岩溶勘察中的应用[D];中南大学;2011年
9 黄靓;混凝土超声CT的数值模拟与试验研究[D];湖南大学;2003年
10 戚秀真;混凝土超声无损检测层析成像技术研究[D];长安大学;2006年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026