基于大地电磁测深数据的大地电位分布计算.docxVIP

基于大地电磁测深数据的大地电位分布计算 第 PAGE \* Arabic \* MERGEFORMAT 4 页第 PAGE \* Arabic \* MERGEFORMAT 4 页 基于大地电磁测深数据的直流接地极地表电位分布及电网偏磁电流计算 摘要：直流接地极入地电流会导致接地极周围地表电位变化，从而引起电网中性点接地变压器偏磁电流。地表电位分布与接地极入地电流、大地电阻率和接地极型式等很多因素有关。其中，大地电阻率是重要因素之一。为了提高地表电位分布计算精度，本文提出基于大地电磁测深数据的大地电阻率模型建模方法；在此基础上，建立±800kV天中直流工程送端电网哈密接地极(E94.2°,N42.8°)为中心约4°×4°范围内的大地电阻率模型，计算研究范围内地表电位分布及哈密电网偏磁电流水平，验证了大地电阻率建模方法的可行性。 关键词：大地电磁；地表电位；偏磁电流；有限元 0引言 直流输电工程单极大地回线运行方式下在地面形成地表电位[1,2]，从而会引起电网变压器直流偏磁问题[3-5]，威胁电网安全[6,7]。 确定地表电位的分布要考虑大地电性结构的实际分布，尤其是深层大地电阻率对地表电位的影响[8-10]。特高压直流输电工程单极运行时接地极入地电流达到几千安，上海庙-山东特高压直流输电工程额定入地电流将达到6250A，其影响范围的空间尺度达到一百多公里甚至几百公里[11,12]，可能包含大地断层、海陆交界和不同地质结构交界面的大地电阻率横向变化，应考虑大地电阻率的空间三维分布特性。文献[13-15]采用大地电阻率典型值考虑了山川、湖泊和海洋对大地地位的影响。文献[16,17]基于深部探测中大地电磁测深数据建立华北地区三维电导率模型并运用有限元计算磁暴感应地电场。现有文献未发现有采用大范围大地电磁测深数据建立大地电阻率模型计算直流接地极所引起的地面电位。 本文采用我国往年所获得的大地电磁测深数据构建天中直流哈密接地极周围大地电阻率模型，应用有限元法求解地表电位分布。根据地表电位的计算结果，基于全节点电压模型[18]计算天中直流工程送端电网偏磁电流，并与实测值进行比较，验证大地电阻率模型和地表电位及偏磁电流计算的有效性，为进一步制定直流输电工程送/受端电网防范接地极侵害的技术方案和接地极址选址提供依据。 1地面电位计算基本原理 直流流入大地后，大地表面的电位分布是一个电流场的计算，流入土壤中的电流与地表电位满足恒定电流场方程： 场方程 （不含电流源） （包含电流源） 在地面 A1 处满足第二类的齐次 Norman 边界条件： 在无穷远参考点处A2 处满足第一类 Dirichlet 边界条件： 根据大地电性构造求解接地极附近地表电位分布问题可以用电阻率不同的各个区域内的 Laplace 方程或 Poisson 方程进行描述，各区域内满足一定的边界条件，求解方程即可获得地表电位分布。有限元法将偏微分方程表征的连续微分函数所在的封闭场域划分为有限个小区域，每个小区域用选定的近似函数来逼近，以获得整个该场域中函数的近似数值。由于有限元方法能够自动满足媒质交界面条件，因此有限元法是一种特别适合于处理多层复杂媒质问题的数值算法。 2我国大地电磁测深数据及大地电阻率建模 大地电阻率数据获取的最理想的方法是按照接地极影响区域进行大规模测量。但是由于接地极影响区域广阔，以及财力物力的限制，无法对影响区域进行完备的电阻率测量。根据已有的大地电磁数据建立接地极周围电阻率模型是一种较为经济的方法。 2.1 我国大地电磁测深数据 2.1.1中国地球物理图集 《中国地球物理图集》[19]电阻率图是根据原地质矿产部系统的大地电磁测深（MT）野外实地调查结果，以及收集由国家地震局、各地质院校、中国科学院等单位绘制中国地学大断面时所做大地电磁测量总共约900个测点整理汇编而成。具体大地电磁剖面如下图所示。 2.1.2中国地学大断面 本文用于我国大地电阻率建模的大地电磁测深数据来源于国家地震局地质断面编委会出版的地学断面的大地电磁测深结果[36,37]。1987年，中国岩石圈委员会决定编制11条中国地学断面(CGT)，并成立了中国地学大断面协调组，由地质矿产部、国家地震局、中国科学院和石油工业部(现为石油天然气总公司)代表组成。 2.1.3深部探测（SinoProbe）数据 深部探测技术与实验研究专项(SinoProbe，2008-2012)是我国历史上实施的规模最大的地球深部探测计划。专项开展了全国4°×4°、华北和青藏高原1°×1°的大地电磁阵列观测。“观测网”的范围和各测点如下图所示。采用多种大地电磁反演方法，得到全国实验区岩石圈不同深度电阻率数据体及三维电磁参数模型。 2.2大地电阻率建模 本文以《中国地球物理图集