西北某铁路隧道采空区地质选线勘察.pdf

·地质勘察·收稿日期：２０１６０７２５基金项目：中铁第五勘察设计院集团有限公司科研开发计划项目（Ｙ２０１２００５）西北某铁路隧道采空区地质选线勘察孙鹏（中铁第五勘察设计院集团有限公司北京１０２６００）摘要在铁路地质选线工作中，采空区选线一直以来是一个重要研究内容。在西北某煤矿专用线铁路采空区地质选线工作中，我们运用资料收集、地质调绘、物探、钻探等综合勘察方法，查明了选线工作中控制性工点Ｄ１Ｋ３３＋１６０隧道的采空区分布范围、特征等，并对穿越采空区的控制性因素安全煤柱宽度和安全顶板厚度进行了分析确定，指导了采空区线路选线和优化工作。关键词采空区铁路隧道地质选线中图分类号Ｕ４５２．１＋３文献标识码Ａ文章编号１００９４５３９（２０１６）１１００２９０５ＧｅｏｌｏｇｉｃａｌＬｉｎｅＳｅｌｅｃｔｉｏｎｏｆａＲａｉｌｗａｙＴｕｎｎｅｌｉｎＧｏａｆｏｆＮｏｒｔｈｗｅｓｔＲｅｇｉｏｎＳｕｎＰｅｎｇ（ＣｈｉｎａＲａｉｌｗａｙＦｉｆｔｈＳｕｒｖｅｙａｎｄＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅＧｒｏｕｐＣｏ．Ｌｔｄ．，Ｂｅｉｊｉｎｇ１０２６００，Ｃｈｉｎａ）ＡｂｓｔｒａｃｔＩｔｈａｓｂｅｅｎａｎｉｍｐｏｒｔａｎｔｒｅｓｅａｒｃｈｃｏｎｔｅｎｔ⁃⁃ｆｏｒｔｈｅｌｏｃａｔｉｏｎｏｆｍｉｎｅｄｏｕｔａｒｅａｉｎｔｈｅｒａｉｌｗａｙｇｅｏｌｏｇｉｃａｌｌｉｎｅｓｅｌｅｃｔｉｏｎ．Ｉｎｔｈｅｇｅｏｌｏｇｉｃａｌｌｉｎｅｓｅｌｅｃｔｉｏｎｏｆａｓｐｅｃｉａｌｒａｉｌｗａｙｌｉｎｅｆｏｒｃｏａｌｍｉｎｅｉｎｇｏａｆｏｆｎｏｒｔｈｗｅｓｔｒｅｇｉｏｎ，ｓｏｍｅｃｏｍｐｒｅｈｅｎｓｉｖｅｓｕｒｖｅｙｍｅｔｈｏｄｓｓｕｃｈａｓｄａｔａｃｏｌｌｅｃｔｉｏｎ，ｇｅｏｌｏｇｉｃａｌｓｕｒｖｅｙｉｎｇａｎｄｃｏｍｐｉｌｉｎｇ，ｇｅｏｐｈｙｓｉｃａｌ⁃ｐｒｏｓｐｅｃｔｉｎｇａｎｄｄｒｉｌｌｉｎｇｗｅｒｅａｐｐｌｉｅｄｔｏｉｄｅｎｔｉｆｙｔｈｅｇｏａｆｄｉｓｔｒｉｂｕｔｉｏｎｒａｎｇｅ，ｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆＤ１Ｋ３３＋１６０ｔｕｎｎｅｌｉｎｔｈｅｃｏｎｔｒｏｌｌｅｄｅｎｇｉｎｅｅｒｉｎｇ．⁃Ｉｎａｄｄｉｔｉｏｎ，ｔｈｅｋｅｙｆａｃｔｏｒｓｏｆｓａｆｅｃｏａｌｐｉｌｌａｒｗｉｄｔｈａｎｄｔｈｉｃｋｎｅｓｓｏｆ⁃ｔｈｅｒｏｏｆｉｎｍｉｎｅｄｏｕｔａｒｅａｗｅｒｅａｎａｌｙｚｅｄａｎｄｄｅｔｅｒｍｉｎｅｄ．Ｉｔｃｏｕｌｄｇｕｉｄｅｔｈｅｌｉｎｅｓｅｌｅｃｔｉｏｎａｎｄｏｐｔｉｍｉｚａｔｉｏｎｉｎｇｏａｆ．Ｋｅｙｗｏｒｄｓｇｏａｆ；ｒａｉｌｗａｙｔｕｎｎｅｌ；ｇｅｏｌｏｇｉｃａｌｌｉｎｅｓｅｌｅｃｔｉｏｎ在铁路地质选线中，采空区地质选线是一个重要研究内容，特别在煤矿铁路专用线勘察设计时，这个问题尤为突出。在对我国西北地区某运煤专用线穿越煤矿采空区段落进行勘察设计时，我们首先应运用综合勘察手段查明采空区的分布情况［１］，再评价采空区的稳定性及对隧道工点的影响，提出立体选线的理念，确定了线路以隧道形式通过采空区的方案。１工程概况该铁路运煤专用线位于内蒙古自治区西南部鄂尔多斯市境内，为当地铁路煤炭物流通道中重要的一段，设计线路全长１２８．１０２ｋｍ。线路所经地区主要为低中山区，地层、岩性多变，水文地质条件复杂，发育有采空区、滑坡、风沙、湿陷性黄土的地质问题；其中线路ＤＫ０～ＤＫ４０段煤矿分布众多，沿线开采煤矿形成的采空区是地质选线勘察的重点、难点［２］。在设计过程中，经过充分研究当地煤矿、采空区分布情况，以及深入的经济效益比选后，初步确定了绕矿方案和穿矿方案，绕矿方案线路较长，且不符合当地的规划要求，但穿矿方案中，Ｄ１Ｋ３３＋１６０隧道穿越煤矿之间的保安煤柱，距最近煤矿采空边界仅３８ｍ。勘察期间我们对Ｄ１Ｋ３３＋１６０隧道进行了专项勘察、采空区稳定性评价、铁路隧道穿越采空区控制因素分析等地质选线工作［３］，为线路方案的设计优化提供了依据。９２铁道建筑技术ＲＡＩＬＷＡＹＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ２０１６（１１）万方数据·地质勘察·２隧道采空区专项勘察为完成Ｄ１Ｋ３３＋１６０隧道重点工程穿越采空区地质选线的勘察设计需求，确保铁路运行安全，我们于２０１０年４月至６月、２０１０年８月至１０月、２０１１年４月至５月开展了煤矿采空区地质选线的专题勘察。勘察过程中采用资料收集、工程地质调绘、物探、岩土测试、水文地质、工程测量等相结合的综合勘察方法［４］，共完成∶１２０００工程地质及水文地质测绘工作３ｋｍ２；钻孔总数３６个，总进尺２４９４．２ｍ，取岩样１０组，水质分析样品２组；高密度电法测线２５３０点／１２６２５ｍ，电测深９７点，地震波ＣＴ１５对。专项综合勘察工作流程如图１所示。图１Ｄ１Ｋ３３＋１６０隧道综合勘察工作流程２．１资料收集及应用在内蒙古国土资源厅、鄂尔多斯国土局、线路沿线煤矿等单位收集了大量相关资料，包括∶１２０万区域地质、水文地质调查报告、区域地质及区域水文地质图∶、１５万区域地质图及报告。我们在勘察工作准备期间及勘察过程中充分收集了测区煤矿采空区井上井下对照图、煤层等高线图及瓦斯地质资料［５］，对本隧道的采空区地质勘察工作有重要参考价值。２．２工程地质调绘地质组对隧道及采空区通过地段进行了１／２０００、１／５００平面图的工程地质调绘。结合访问，工程地质调绘工作初步查明了隧址区的工程地质条件和水文地质条件。对隧址区的含煤地层、煤矿采空区等不良地质体范围进行了调查，并进入了矿井下核实采空区分布范围。其中，采空区处断面图由测量人员和地质人员协同绘制，实测岩层分界线，逐段填绘工程地质、水文地质界线［６］。２．３钻探勘察中投入多台—ＸＹ１８０型—、ＸＹ１５０型—、ＸＹ３００型钻机承担Ｄ１Ｋ３３＋１６０隧道采空区的钻探工作。钻探按《铁路工程地质钻探规程》要求进行。２．４物探本次专项勘察采用了综合物探方法，包括电测深、高密度电法、瞬变电磁法、孔间ＣＴ法等［７］，通过上述方法先初步确定采空区异常部位，再对部位进行钻探工作，运用孔间ＣＴ法定位采空区具体分布位置；期间均采用了两种或以上的物探方法互相验证、互相补充，提高物探解译精度。本次物探工作的野外数据采集、资料整理和解释严格按照《铁路工程物理勘探规程》的要求进行。２．５其他方法专项勘察期间，我们进行了含煤地层有害气体测定和评价、岩土测试、水文地质调查、水文地质实验、工程测量等工作［８］。专项地质各项工作内容、精度和资料整编均符合相关技术规范、规程的要求，质量良好，满足勘察技术要求和施工图设计需要，达到了预期目的。３隧道与煤矿采空区的关系本次隧道地质选线工作，根据查明的隧道与采空区的位置关系，进行采空区安全煤柱宽度分析、计算下穿采空区安全顶板厚度，确定隧道穿越煤矿采空区的控制要素，指导线路方案的选择。Ｄ１Ｋ３３＋１６０隧道与沿线煤矿、采空区边界位置关系如图２所示。图２Ｄ１Ｋ３３＋１６０隧道与沿线煤矿、采空区边界位置示意图３．１Ａ煤矿采空区Ａ煤矿与线路关系如图３所示，线路经过范围于２００５年之前采用巷道方式开采，后因某种原因停采，所采煤层为３＃煤层，煤层底板标高１３６１～１３６２ｍ，０３铁道建筑技术ＲＡＩＬＷＡＹＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ２０１６（１１）万方数据·地质勘察·煤层平均层厚约４～５ｍ。经测绘和调查，Ａ煤矿采空区沉陷影响主要表现为地面裂缝，基本分布在Ｄ１Ｋ３２＋３８０～Ｄ１Ｋ３２＋４６０，主要裂缝有２组。根据勘察成果：采空区分布在Ｄ１Ｋ３２＋３４０～Ｄ１Ｋ３２＋７１０段，采空标高底板为１３６１～１３６２ｍ，开采厚度约３ｍ。图３Ａ煤矿采空区边界、地面裂缝与线路位置示意图３．２Ｂ煤矿采空区Ｂ煤矿与线路关系如图４所示，该矿所采煤层为５＃煤层，采煤方式为综采。经测绘和调查，Ｂ煤矿采空区在地表的主要破坏形式为地面裂缝，基本分布在Ｄ１Ｋ３３＋９２０～Ｄ１Ｋ３４＋０８０右侧７０～１１００ｍ范围内，主要裂缝共５组。根据综合勘察成果：采空区分布在Ｄ１Ｋ３３＋９００～Ｄ１Ｋ３４＋７７０段右侧５０ｍ以外，该处采空区分布在１３１０～１３２０ｍ，高度约２～４ｍ。图４Ｂ煤矿采空区边界、地面裂缝与线路位置示意图３．３Ｃ煤矿采空区Ｃ煤矿与线路关系如图５所示，该煤矿位于线路左侧４３～２６５ｍ范围，该矿所采煤层为２＃、４＃煤层；经测绘和调查，其采空区在地表的主要破坏形式为地面裂缝，基本分布在Ｄ１Ｋ３３＋８００～Ｄ１Ｋ３４＋３５０，主要裂缝共３组。根据综合勘察成果沿线路方向在Ｄ１Ｋ３２＋８００～Ｄ１Ｋ３４＋３５０段左侧４３～２６５ｍ外为采空区范围，其中Ｄ１Ｋ３２＋８００～Ｄ１Ｋ３３＋７５０段采空底板标高约１３６２～１３７４ｍ，厚度约３ｍ，Ｄ１Ｋ３３＋７５０～Ｄ１Ｋ３４＋３５０段采空底板标高为１３１０～１３３０ｍ，厚约３ｍ，采空区地面均出现变形裂缝。图５Ｃ煤矿采空区边界、地面裂缝与线路位置示意图４铁路隧道穿越采空区的控制因素分析本次隧道地质选线工作，根据查明的隧道与采空区的位置关系，进行采空区安全煤柱宽度分析、计算下穿采空区安全顶板厚度，确定隧道穿越煤矿采空区的控制要素，指导线路方案的选择。４．１保安煤柱计算通过调查、物探、钻探等，煤矿采空区基本同采矿权登记面积相同。根据调查、煤矿矿产详查报告等资料，该区域煤层倾角较小，约３～７°，勘测期间实测岩层倾向为Ｎ６０°Ｗ／６°Ｓ，根据《铁路工程地质手册》中大面积采空区的地表移动变形预测一章中，关于移动盆地边界预测的走向边界角的公式计算：δ０＝δ－１５°（其中δ为基岩移动角）按照手册要求，侏罗系砂岩夹泥岩地层走向移动角取７２°。根据本次勘察查明的采空区分布情况及隧道范围内的工程地质条件等［９］，隧道Ｄ１Ｋ３２＋２４０～Ｄ１Ｋ３３＋１１０段Ａ煤矿采空区位于隧道上方，对隧道顶板稳定性有一定影响，控制该段线路方案选择；隧道Ｄ１Ｋ３４＋３００～３５０段左侧Ｃ煤矿采空区距离隧道最近处近４４～５０ｍ，采空标高１３１０ｍ，对隧道工程有影响；隧道右侧Ｂ煤矿巷道距隧道水平距离约８３ｍ，对隧道影响较小。Ｄ１Ｋ３３＋１６０隧道沿线各煤矿主要煤层保安煤柱宽度分析如下所述。４．１．１Ａ煤矿（Ｄ１Ｋ３２＋２４０～Ｄ１Ｋ３３＋１１０）段采空区保安煤柱３＃煤层：煤层底板高程为１３６１～１３７８ｍ，煤层厚度约２ｍ，该层煤已开采，采空区位于隧道洞顶上方，线路所经分为采空区底板高程１３６１～１３６３ｍ，对隧１３铁道建筑技术ＲＡＩＬＷＡＹＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ２０１６（１１）万方数据·地质勘察·道顶板围岩稳定性有一定影响，应注意控制采空区底板与隧道顶板之间的距离。计算示意图见图６。图６ＤＫ３２＋４５０处Ａ煤矿采空区移动盆地计算示意图４．１．２Ｂ煤矿、Ｃ煤矿（Ｄ１Ｋ３３＋１１０～Ｄ２Ｋ３４＋７７０）段保安煤柱３＃煤层：煤层底板高程为１３８０～１３９３ｍ，煤层厚约３ｍ，该层煤层已经开采，采空边界距离线路最近处位于隧道路肩上方３５ｍ左侧４４ｍ，对隧道影响较小。５＃煤层：煤层底板高程为１３１０～１３３５．５１ｍ，煤层厚约３．１ｍ，线路所经范围煤层底板标高为１３１０ｍ，该层煤已经开采，采空边界距离线路最近处为隧道设计路肩面以下３５ｍ左侧４４ｍ，保安煤柱范围左侧５１ｍ，右侧５５ｍ。图７隧道下穿采空区安全顶板厚度分析模型４．２隧道下穿采空区安全顶板厚度确定方法采空区底板的破坏深度是确定下穿隧道合理安全顶板厚度的重要因素，采空区底板破坏与采空区垮落后支撑压力、采空区开挖后底板应力分布有关［１０］。煤矿开采后采空区底板围岩破坏与下穿采空区隧道围岩破坏相互影响，围岩破坏区之间几何关系如图７所示。由图７所示围岩破坏区几何关系分析，隧道下穿采空区安全顶板厚度Ｈ可按下式确定：Ｈ≥Ｄ＋ｓ＋（Ｒ－ｒ）式中，Ｄ为采空区底板破坏深度；ｓ为原岩保护层厚度；Ｒ为隧道中心至拱顶破裂区范围；ｒ为隧道半径。由式中可知，隧道下穿采空区安全顶板厚度与采空区底板破坏区深度、隧道围岩松动与破坏区范围、原岩破坏层厚度相关，这三个参数是确定安全顶板厚度的基本依据。４．２．１采空区底板破坏深度采空区下底板围岩破坏深度可通过弹性理论进行分析，采空区底板最大破坏深度可按下式求解：Ｄ＝Ｌ·ｃｏｓφ２·ｃｏｓ（π／４＋φ／２）·ｅ（π４＋φ２）·ｔａｎφ式中，Ｌ为煤体区服区宽度，按经验取Ｌ＝０．０１５Ｈｃ，Ｈｃ为煤层底板埋深；φ为底板岩层内摩擦角。４．２．２原岩保护层厚度原岩保护层是采空区底板破坏区与隧道拱部破坏区之间的预留安全岩层［１１］。顶水采煤相关研究给出了确定该层厚度的经验公式：ｓ＝１．５ｈ１／ｆ＋ｃ式中，ｓ为原岩保护层厚度；ｈ１为坑道宽度；ｆ为岩石坚固系数；ｃ为岩层强风化带厚度。４．２．３采空区下隧道围岩应力的弹塑性分析采空区下隧道开挖应力重分布后，围岩稳定处于弹性状态时成为弹性应力场。当隧道距离地面埋深较大、围岩应力较大或围岩的强度偏小时，隧道围岩的强度条件小于或等于围岩的应力分布时，围岩会从弹性应力状态转换为塑性应力状态。径向应力：σｒｅ＝ｐ（１－Ｒ２ｒ２）＋（ｐ１＋ｃｃｏｔφ）（Ｒａ）２ｓｉｎφ１－ｓｉｎφ－ｃｃｏｔφéëêêùûúúＲ２ｒ２切向应力：σｔｅ＝ｐ（１＋Ｒ２ｒ２）－（ｐ１＋ｃｃｏｔφ）（Ｒａ）２ｓｉｎφ１－ｓｉｎφ－ｃｃｏｔφéëêêùûúúＲ２ｒ２塑性区的应力：径向应力：σｒｐ＝（ｐ１＋ｃｃｏｔφ）（ｒ／ａ）２ｓｉｎφ１－ｓｉｎφ－ｃｃｏｔφ切向应力：σｔｐ＝（ｐ１＋ｃｃｏｔφ）１＋ｓｉｎφ１－ｓｉｎφ（ｒａ）２ｓｉｎφ１－ｓｉｎφ－ｃｃｏｔφ塑性区半径：Ｒ＝ａ（ｐ＋ｃｃｏｔφ）（１－ｓｉｎφ）ｐ１＋ｃｃｏｔφéëêêùûúú２ｓｉｎφ１－ｓｉｎφ式中，ｐ为原岩应力；ｐ１为支护反力；ｒ为所求应力处２３铁道建筑技术ＲＡＩＬＷＡＹＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ２０１６（１１）万方数据·地质勘察·的半径；Ｒ为塑性区半径；φ为内摩擦角；ｃ为黏聚力。４．２．４隧道拱顶围岩破裂区隧道拱顶围岩破裂区在工程实践中，常采用普氏压力拱法进行预估。隧道拱顶上部岩体塌落的高度，可按下式计算：Ｈｔ＝［Ｂ２＋ｈｔａｎ（４５°－φ２）］／ｆ式中，Ｈｔ为隧道塌落拱高度；Ｂ为隧道宽度；ｈ为隧道高度；φ为内摩擦角；ｆ为岩石坚固系数。采用前文介绍的基于多因素的理论分析方法，确定Ｄ１Ｋ３３＋１６０隧道下穿采空区安全顶板厚度。计算时，煤层底板埋深３０ｍ，坑道宽度３ｍ，采空区底板岩层内摩擦角３２°，岩石坚固系数１．０，隧道宽度为１１．８ｍ，隧道净高１０．５ｍ。按多因素理论分析方法求得的采空区底板岩层最大破坏深度为１．１４ｍ，隧道塌落拱高度１１．７２ｍ，原岩保护层厚度２．６ｍ，得出安全顶板最小厚度１５．４６ｍ。４．３立体选线理念隧道与其下方采空区边界底板的水平距离称为隧道上跨采空区保安煤柱；隧道顶部距上方采空区底板岩层厚度表示了下穿隧道与采空区的竖向位置关系，又称为隧道下穿采空区安全顶板厚度。Ｄ１Ｋ３３＋１６０隧道穿越矿区段选线时，通过线路优化规避Ｂ煤矿、Ｃ煤矿采空区的影响，降低了部分隧道穿越采空区的工程风险，但在部分段落须下穿Ａ煤矿采空区，通过降低铁路标高得以增大下穿采空区时安全顶板厚度，施工风险降低［１２］。因此，隧道穿越采空区设计方案应从投资、施工风险、工期等方面进行比选，而保安煤柱宽度、安全顶板厚度的确定是整个设计的关键，成为控制隧道投资、施工风险、工期的决定性因素。５结束语（１）Ｄ１Ｋ３３＋１６０隧道运用资料收集分析、地质调绘、物探、钻探等综合勘察手段查明了采空区的分布范围、特征等，为线路方案选线设计提供了有力依据。（２）Ｄ１Ｋ３３＋１６０隧道在地质选线过程中，分析线路平面绕避安全煤柱设计宽度、下穿采空区安全顶板厚度等控制性因素，采用立体选线的思路，为线路穿越煤矿采空区段落提供了准确、合理的指导。（３）设计专业通过采空区地质选线勘察的成果，对线路方案和隧道防护措施进行了优化。现Ｄ１Ｋ３３＋１６０隧道已经建成，并通车两年，未出现工程地质病害问题。参考文献［１］李鹏．隧道穿越采空区施工技术研究［Ｊ］．桥隧施工与机械，２０１３（１０）：８６．［２］孙鹏．某在建铁路工程滑坡成因分析与工程治理浅析［Ｊ］．科技创新与应用，２０１２（３３）：２０２．［３］刘士春．永固煤矿采场底板破坏深度及防突水技术研究［Ｊ］．煤炭技术，２００８（１）：７８．［４］ＴＢ１００１２－２００７铁路工程地质勘察规范［Ｓ］．［５］曹祥．山西中南部铁路通道吕梁山越岭段方案研究［Ｊ］．铁路标准设计，２０１４（２）：８．［６］杜宇本．成兰铁路主要地质灾害与地质选线［Ｊ］．铁道工程学报，２０１２（８）：１１．［７］周革．采空区铁路工程地质选线［Ｊ］．铁路勘察，２０１２（３８）：４１．［８］孙宝忠．包西铁路毛乌素沙地大面积煤矿采空区地质选线研究［Ｊ］．铁路标准设计，２０１４（６）：１．［９］白雪源．高频大地电磁法在云雾山隧道勘察中的应用［Ｊ］．铁道建筑技术，２０１３（８）：３５．［１０］苗庆库．弹性波ＣＴ在岩溶勘察中的应用［Ｊ］．铁道建筑技术，２００８（Ｓ１）：５３５．［１１］覃羡安．常家山隧道穿越煤系地层及采空区段的施工过程数值模拟研究［Ｄ］．西安：西安科技大学，２００５．［１２］铁道部第一勘测设计院．铁路工程地质手册［Ｍ］．北京：中国铁道出版社，１９９９：５７７－５７８．������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������（上接第４页）［１２］孙宝忠，李光伟，吴为禄．三维影像遥感技术在滇藏线预可行性研究中的应用［Ｊ］．铁道工程学报，２００２（３）：７－９．［１３］冯光胜．三维遥感技术在十宜铁路地质选线中的应用［Ｊ］．地理空间信息，２０１０，８（２）：９６－９９．［１４］高山．三维遥感可视化技术在高盖山隧道工程地质调查中的应用研究［Ｊ］．铁道勘察，２００７，３３（６）：５０－５４．［１５］高文峰，甘俊，王长进．ＤＩｎＳＡＲ沉降监测技术在铁路勘察设计中的应用［Ｊ］．遥感信息，２０１５，３０（５）：８３－８７．［１６］谭衢霖，魏庆朝，杨松林，等．铁路工程地质遥感雷达图像应用分析［Ｊ］．铁道工程学报，２００９（１）：１７－２２．３３铁道建筑技术ＲＡＩＬＷＡＹＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ２０１６（１１）万方数据