CN106021651B - 一种复杂地质条件下直流电流的有限元计算方法 - Google Patents
一种复杂地质条件下直流电流的有限元计算方法 Download PDFInfo
- Publication number
- CN106021651B CN106021651B CN201610301433.XA CN201610301433A CN106021651B CN 106021651 B CN106021651 B CN 106021651B CN 201610301433 A CN201610301433 A CN 201610301433A CN 106021651 B CN106021651 B CN 106021651B
- Authority
- CN
- China
- Prior art keywords
- model
- current
- geological
- result
- finite element
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种复杂地质条件下直流电流的有限元计算方法,该方法根据真实地质环境在ANSYS中建立模型进行计算,该方法包括下列步骤:建立原始模型,对原始模型进行直流电流的有限元计算;加入一种地质模型,进行直流电流的有限元计算,并将计算得到的结果与原始模型中得到的结果进行比较;改变加入地质模型的参数,进行计算,并将得出的结果与未改变参数的结果进行比较;加入两种或两种以上的地质模型,进行计算,将得到的结果与之前的结果进行比较;总结比较结果,得出不同地质条件对直流电流的有限元计算的影响规律;根据规律与当前地质条件进行直流电流的有限元计算。与现有技术相比,本发明具有实用性高以及计算结果准确等优点。
Description
技术领域
本发明涉及高压直流输电单极大地方式运行引起的直流偏磁的计算与治理,尤其是涉及一种复杂地质条件下直流电流的有限元计算方法。
背景技术
随着我国电力系统不断发展,超特高压直流输电线路逐渐成为跨地区电力输送的重要走廊。直流输电工程通常采用双极系统接线方式,正常运行时,直流电流不流经大地(或海水)。但在系统调试、检修或故障时,转为单极大地回线运行方式,大地(或海水)相当于直流输电线路的一根导线,其电流场会使地电位发生变化。对中心点接地的交流系统而言,处于不同直流电位的变电站经输电线路、变压器绕组构成直流回路,使变压器产生直流偏磁现象,严重影响供电质量和供电安全。
直流偏磁的计算和治理建立在对接地极散流特性准确把握的基础之上,而接地极的散流特性受其附近地质条件的影响很大,需要考虑土壤的水平分层以及远离接地极时土壤电阻率的变化情况。在综合考虑土壤的水平分层和垂直分层时,接地电流场是一个复杂的物理系统,某一区域地质条件的变化会引起整体的电流场发生显著改变,进而对直流偏磁的预测结果产生影响。而目前对接地极散流特性的研究主要集中在计算方法和水平分层的影响方面,对地质结构和地质条件的影响规律研究还比较缺乏。因而提出利用基于矩量法的有限元计算对接地极的散流特性进行仿真研究,进而通过加入和改变不同的地质条件来系统地分析总结山体、盆地、河流、湖泊等对散流特性的影响。
发明内容
本发明的目的是针对上述问题提供一种复杂地质条件下直流电流的有限元计算方法。
本发明的目的可以通过以下技术方案来实现:
1.一种复杂地质条件下直流电流的有限元计算方法,该方法根据真实地质环境在ANSYS中建立模型进行计算,其特征在于,该方法包括下列步骤:
1)建立原始模型,对原始模型进行直流电流的有限元计算;
2)对步骤1)中得到的原始模型加入一种地质模型,进行直流电流的有限元计算,并将计算得到的结果与步骤1)中得到的结果进行比较;
3)改变步骤2)中加入地质模型的参数,进行直流电流的有限元计算,并将得出的结果与步骤2)中得到的结果进行比较;
4)对步骤1)得到的原始模型加入两种或两种以上的地质模型,进行直流电流的有限元计算,并将得到的结果与步骤1)和步骤2)中得到的结果进行比较;
5)总结步骤2)、3)与步骤4)中的比较结果,得出不同地质条件对直流电流的有限元计算的影响规律;
6)根据步骤5)中得出的规律与当前地质条件,对当前地质条件进行直流电流的有限元计算。
所述直流电流的有限元计算包括下列步骤:
11)根据建立的模型选取模型单元,添加单元材料属性;
12)根据步骤11)中选取的模型单元,对建立的模型进行剖分;
13)对步骤12)中进行过剖分的模型施加电流载荷以及边界条件;
14)选取求解算法对步骤13)中施加过电流载荷以及边界条件的模型进行求解,得出整个模型的电位分布以及电流密度分布。
所述模型单元为十节点四面体SOLID232。
所述单元材料属性为电阻率属性。
所述求解算法为JCG算法。
所述原始模型包括土壤、海洋、表层基岩以及底层基岩。
所述地质模型包括山体模型、河流模型、湖泊模型以及盆地模型。
所述地质模型的参数包括大小、电阻率以及方向。
与现有技术相比,本发明通过加入和改变不同的地质条件来系统的分析总结山体、盆地、河流和湖泊等对散流特性的影响,全面的考虑了不同的地质结构与地质条件,使得对于直流偏磁的计算更加准确,同时也具有较高的实用价值;同时总结了地质结构和地质条件对于直流偏磁计算的影响规律,便于在遇到相似或有规律的地质情况时避免重复计算,节省了计算时间与计算量,节省了成本。
附图说明
图1为原始模型俯视图;
图2为原始模型侧视图;
图3为山体模型图;
图4为湖泊与盆地模型图;
图5为河流模型图;
图6为加入山体模型后的电位等位线图;
图7为加入湖泊模型后的电位等位线图;
图8为加入河流模型后的电位等位线图;
图9为加入盆地模型后的电位等位线图;
图10为复杂地质条件下直流电流的有限元计算方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图10所示为复杂地质条件下直流电流的有限元计算方法流程图,该方法根据真实地质环境在ANSYS中建立模型进行计算,该方法包括下列步骤:
(1)建立原始模型,对原始模型进行直流电流的有限元计算;
(2)对(1)中得到的原始模型加入一种地质模型,对新模型进行直流电流的有限元计算,并将计算得到的结果与(1)中得到的结果进行比较;
(3)改变(2)中加入地质模型的参数,对新模型进行直流电流的有限元计算,并将得出的结果与(2)中得到的结果进行比较;
(4)对(1)得到的原始模型加入两种或两种以上的地质模型,将得到的结果与(1)和(2)中得到的结果进行比较;
(5)总结(2)、(3)与(4)中的比较结果,得出规律;
(6)根据(5)中得出的规律与当前地质条件,对当前地质条件进行直流电流的有限元计算。
在上述步骤中,涉及到的具体技术方案如下:
(a)建立模型
根据真实环境建立一个模型,如图1及图2所示,此模型中包括了土壤、海洋、表层基岩及底层基岩。
在此模型基础上,可添加山体、河流、湖泊以及盆地模型等,以接近实际的地理环境。山体模型如图3所示,山体底部连接着底层基岩,顶部突处于土壤;湖泊与盆地模型如图4所示,为镶嵌在土壤中的圆柱形,盆地的模型与湖泊类似,盆地中的介质为空气;河流模型如图5所示,河流也是镶嵌在土壤中。这些模型可以单独添加,也可以一起添加在上面的模型中。
(b)选取单元类型并添加材料属性
为了控制运算量,减少自由度,选取十节点四面体SOLID232作为模型单元。ANSYS中稳态电流场分析只需要添加电阻率属性,因此添加单元材料属性如表1所示。
表1单元材料属性
(c)模型分析
为了得到平滑的电流分布曲线,需对模型进行较为细密的剖分。综合考虑计算结果精度和硬件配置,将模型剖分成百万数量级进行计算。Solid 232是四面体型单元,因此网格划分只能采用自由剖分。在ANSYS中打开智能网格划分选项,利用lesize命令设定各边剖分精度,使单元和节点数能保证计算结果的准确性。
(d)施加载荷及边界条件
给模型施加电流载荷,电流大小为3000A,模型四周施加电压为0的边界条件。
(e)模型求解
由于模型自由度达到百万级,而计算平台的内存和磁盘大小有一定限制,因此有必要分析各种算法对硬件的要求,以便选择合适的算法。ANSYS提供的求解算法主要分为3类:一是内存占用大,磁盘占用小;二是内存占用小,磁盘占用大;三是内存和磁盘占用较为均衡。由于本模型计算问题为电流场单场领域,计算平台内存120G,磁盘空间200G,综合所有因素,选择JCG算法进行计算。
(f)计算结果
通过以上的建模计算,得出整个模型的电位分布与电流密度分布,通过比较得出各种地质条件对散流的影响。电位等位线图如图6~图9所示。
Claims (7)
1.一种复杂地质条件下直流电流的有限元计算方法,该方法根据真实地质环境在ANSYS中建立模型进行计算,其特征在于,该方法包括下列步骤:
1)建立原始模型,对原始模型进行直流电流的有限元计算,所述原始模型包括土壤、海洋、表层基岩以及底层基岩;
2)对步骤1)中得到的原始模型加入一种地质模型,进行直流电流的有限元计算,并将计算得到的结果与步骤1)中得到的结果进行比较;
3)改变步骤2)中加入地质模型的参数,进行直流电流的有限元计算,并将得出的结果与步骤2)中得到的结果进行比较;
4)对步骤1)得到的原始模型加入两种以上的地质模型,进行直流电流的有限元计算,并将得到的结果与步骤1)和步骤2)中得到的结果进行比较;
5)总结步骤2)、3)与步骤4)中的比较结果,得出不同地质条件对直流电流的有限元计算的影响规律;
6)根据步骤5)中得出的规律与当前地质条件,对当前地质条件进行直流电流的有限元计算。
2.根据权利要求1所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述直流电流的有限元计算包括下列步骤:
11)根据建立的模型选取模型单元,添加单元材料属性;
12)根据步骤11)中选取的模型单元,对建立的模型进行剖分;
13)对步骤12)中进行过剖分的模型施加电流载荷以及边界条件;
14)选取求解算法对步骤13)中施加过电流载荷以及边界条件的模型进行求解,得出整个模型的电位分布以及电流密度分布。
3.根据权利要求2所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述模型单元为十节点四面体SOLID232。
4.根据权利要求2所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述单元材料属性为电阻率属性。
5.根据权利要求2所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述求解算法为JCG算法。
6.根据权利要求1所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述地质模型包括山体模型、河流模型、湖泊模型以及盆地模型。
7.根据权利要求1所述的复杂地质条件下直流电流的有限元计算方法,其特征在于,所述地质模型的参数包括大小、电阻率以及方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610301433.XA CN106021651B (zh) | 2016-05-09 | 2016-05-09 | 一种复杂地质条件下直流电流的有限元计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610301433.XA CN106021651B (zh) | 2016-05-09 | 2016-05-09 | 一种复杂地质条件下直流电流的有限元计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021651A CN106021651A (zh) | 2016-10-12 |
CN106021651B true CN106021651B (zh) | 2019-03-22 |
Family
ID=57098887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610301433.XA Active CN106021651B (zh) | 2016-05-09 | 2016-05-09 | 一种复杂地质条件下直流电流的有限元计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021651B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109858074B (zh) * | 2018-12-13 | 2023-06-13 | 云南电网有限责任公司电力科学研究院 | 一种基于有限体积法的变压器油流涌动的仿真模拟方法 |
CN112668212B (zh) * | 2020-09-02 | 2023-03-24 | 国网内蒙古东部电力有限公司检修分公司 | 基于有限元的不同土壤模型下接地极溢流特性分析方法 |
CN114706012A (zh) * | 2022-03-17 | 2022-07-05 | 云南电网有限责任公司电力科学研究院 | 一种变压器接地电流的计算方法、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593523A (zh) * | 2013-11-12 | 2014-02-19 | 国网上海市电力公司 | 基于有限元理论的多直流落点条件下直流偏磁抑制方法 |
CN104090169A (zh) * | 2014-06-25 | 2014-10-08 | 河海大学 | 一种测量坡面分层土壤结构中地网接地电阻的方法 |
CN104408272A (zh) * | 2014-12-19 | 2015-03-11 | 国网上海市电力公司 | 一种山区环境下的直流偏磁影响评价系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8190368B2 (en) * | 2008-04-07 | 2012-05-29 | Baker Hughes Incorporated | Method of finite-element discretization in heterogeneous and highly conductive grid cells |
-
2016
- 2016-05-09 CN CN201610301433.XA patent/CN106021651B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593523A (zh) * | 2013-11-12 | 2014-02-19 | 国网上海市电力公司 | 基于有限元理论的多直流落点条件下直流偏磁抑制方法 |
CN104090169A (zh) * | 2014-06-25 | 2014-10-08 | 河海大学 | 一种测量坡面分层土壤结构中地网接地电阻的方法 |
CN104408272A (zh) * | 2014-12-19 | 2015-03-11 | 国网上海市电力公司 | 一种山区环境下的直流偏磁影响评价系统及方法 |
Non-Patent Citations (3)
Title |
---|
基于有限元分析的直流偏磁电流计算模型研究;傅晨钊;《华东电力》;20121031;第40卷(第10期);第1739-1742页 |
复杂地形三维直流电阻率有限元数值模拟;李勇;《地球物理学进展》;20090630;第24卷(第3期);第1039-1046页 |
大电流及复杂地质结构下接地系统散流规律及参数计算;曹晓斌;《万方学术论文库》;20111130;第1-95页 |
Also Published As
Publication number | Publication date |
---|---|
CN106021651A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Akpınar et al. | Wave energy potential along the south-east coasts of the Black Sea | |
CN106021651B (zh) | 一种复杂地质条件下直流电流的有限元计算方法 | |
CN103593523A (zh) | 基于有限元理论的多直流落点条件下直流偏磁抑制方法 | |
Sanchez et al. | Floating vs. bottom-fixed turbines for tidal stream energy: A comparative impact assessment | |
Overbye et al. | The use of geographic data views to help with wide-area electric grid situational awareness | |
Guo et al. | Impact of the EHV power system on geomagnetically induced currents in the UHV power system | |
CN104573234B (zh) | 一种基于ArcGIS的流域库容密度空间化方法 | |
CN103699960A (zh) | 一种中压配电网运行态势可视化方法 | |
Sabiha et al. | Assessment of grounding grid for enhancing wind turbine service sustainability | |
CN109975596A (zh) | 一种单极大地运行方式下地中电流分布研究方法 | |
Xia et al. | Comparison of morphodynamic models for the Lower Yellow River 1 | |
Zhang et al. | Hydrologic information extraction based on arc hydro tool and DEM | |
Syed et al. | Simultaneous optimal placement of DGs and fixed capacitor banks in radial distribution systems using BSA optimization | |
CN102201669A (zh) | 一种±660kV阀厅直流金具选型方法 | |
CN109697521A (zh) | 基于任意断面构建洪水预报方案的方法 | |
CN104408272B (zh) | 一种山区环境下的直流偏磁影响评价系统及方法 | |
Bujack et al. | Visualizing electrical power systems as flow fields | |
CN103647270A (zh) | ±800kV特高压直流线路接地的设置方法及系统 | |
Kara et al. | Optimum grounding grid design by using genetic algorithms | |
Ma et al. | Computerized analysis of grounding plates in multilayer soils | |
CN108763790A (zh) | 一种基于扩展临界阻尼调整法的电力系统电磁暂态仿真方法 | |
CN114879788A (zh) | 一种确定变电站直流偏磁电流分布的方法及系统 | |
Fang et al. | Feasibility study on using jacket structure as natural grounding electrode of offshore wind turbines | |
Liu et al. | Research on potential distributions of soil surface caused by ground electrode of HVDC system in heterogeneous soil structure | |
CN107045488A (zh) | 一种基于Bezier曲面四边形参数方程的曲面边界元计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |