CN112765773B - Method for determining public electric field exposure dose near extra-high voltage direct current transmission line - Google Patents

Method for determining public electric field exposure dose near extra-high voltage direct current transmission line Download PDF

Info

Publication number
CN112765773B
CN112765773B CN202011564338.1A CN202011564338A CN112765773B CN 112765773 B CN112765773 B CN 112765773B CN 202011564338 A CN202011564338 A CN 202011564338A CN 112765773 B CN112765773 B CN 112765773B
Authority
CN
China
Prior art keywords
human body
direct current
current transmission
transmission line
high voltage
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
Application number
CN202011564338.1A
Other languages
Chinese (zh)
Other versions
CN112765773A (en
Inventor
翟国庆
项君莉
董丽
郭婧怡
赵菲菲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202011564338.1A priority Critical patent/CN112765773B/en
Publication of CN112765773A publication Critical patent/CN112765773A/en
Application granted granted Critical
Publication of CN112765773B publication Critical patent/CN112765773B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0857Dosimetry, i.e. measuring the time integral of radiation intensity; Level warning devices for personal safety use
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The invention relates to a method for determining public electric field exposure dose near an extra-high voltage direct current transmission line, and belongs to the technical field of environmental exposure dose evaluation. According to simulation calculation results, respectively fitting and establishing a relation model between the human body exposure field intensity E in the outdoor, indoor and automobile, the voltage grade U of the extra-high voltage direct current transmission line, the ground-to-ground height H of the wire, the loop number g and the horizontal distance x between the center of the line and the human body, determining the direct current transmission line nearest to the human body and the position state (outdoor, indoor or in automobile) of the human body according to the mobile phone positioning signal, and calculating the public electric field exposure dose. The public electric field exposure dose caused by the extra-high voltage direct current transmission line can be determined without adopting an electric field dosimeter, and the method is more convenient and faster and has higher universality.

Description

Method for determining public electric field exposure dose near extra-high voltage direct current transmission line
Technical Field
The invention relates to the technical field of environmental exposure dose evaluation, in particular to a method for determining public electric field exposure dose near an extra-high voltage direct current transmission line.
Background
The development of the ultra-high voltage transmission project can effectively promote the interconnection and intercommunication of energy sources. Compared with the ultra-high voltage alternating current transmission, the ultra-high voltage direct current transmission has certain advantages in the aspects of economy and technology, compared with the ultra-high voltage alternating current transmission line with the same voltage class, the capacity of the direct current transmission line is higher, and when the transmission distance is larger than 1000km, the economic benefit of the direct current transmission is better. At present, the voltage class of extra-high voltage direct current transmission reaches 1100kV, and the maximum ground synthetic field intensity under a transmission line can reach 30 kV/m. When a human body moves in an area with high electric field intensity, the skin feels slight stabbing pain, and some sensitive people even feel obvious electric shock. The population density of China is high, the direct-current line mileage spanning farmlands, highways, mountains and even direct neighboring residential houses is larger and larger along with the construction of an extra-high voltage direct-current transmission main line network, and the possible electric field exposure strength and dosage of the public near the line are increased. Health risks caused by long-term exposure to electric fields near ultra-high voltage direct current transmission lines are of great public concern.
The human health risks of electric fields near ultra-high voltage transmission lines depend on the frequency, strength, and exposure time of the electric field. The electric field intensity near the dc line can be measured by using a field intensity measuring instrument, such as a high voltage dc electric field detecting device disclosed in patent application publication No. CN110794225A, and a device and method for detecting the electric field intensity under the extra-high voltage dc transmission line disclosed in patent application publication No. CN 111579889A. However, the above patent applications all measure the field strength of the human body at the non-measured point, and the distribution of the electromagnetic field is closely related to the boundary conditions, so the presence of the human body can cause the distortion of the electric field, and the electric fields measured at the measured points are completely different if there is a human body, that is, the difference between the electric field strength under the power transmission line measured by the field strength measuring instrument and the actual exposed field strength of the human body is large. It is noted in the literature that the electric field exposure dose near the transmission line can be measured using an electric field dosimeter. The electric field dosimeter can be worn on various parts of a human body in a head-wearing manner, an arm-wearing manner, a pocket type manner and the like. However, electric field dosimeters are costly and inconvenient to wear on a daily basis, and are generally not purchased and worn by the public for determining electric field exposure doses in daily life.
Disclosure of Invention
The invention aims to provide a method for determining public electric field exposure dose near an extra-high voltage direct current transmission line, which can determine the public electric field exposure dose caused by the extra-high voltage direct current transmission line without adopting an electric field dosimeter, and is more convenient and more universal.
In order to achieve the above object, the present invention provides a method for determining public electric field exposure dose in the vicinity of an extra-high voltage direct current transmission line, comprising the steps of:
1) constructing an extra-high voltage direct current transmission line-building or vehicle-earth-human body electromagnetic simulation model, respectively establishing a relation model between the human body exposed field intensity E and the voltage grade U of the extra-high voltage direct current transmission line, the ground clearance H of a lead, the number g of loops and the horizontal distance x between the center of the line and the human body in the outdoor, indoor and automobile through simulation calculation, wherein the relation model is respectively E out =f out (U,H,g,x)、E in-b =f in-b (U, H, g, x) and E in-c =f in-c (U,H,g,x);E out Exposing field intensity for outdoor human body; e in-b For the indoor personBody exposure field strength; e in-c Exposing the field intensity for human bodies in the automobile;
2) acquiring voltage grades and loop numbers of all domestic ultrahigh voltage direct current transmission lines, positions of all iron towers on a line and stringing height information of the iron towers, and establishing an ultrahigh voltage direct current transmission engineering information base;
3) according to the position of an iron tower of the extra-high voltage direct current transmission line, numbering the iron tower from south to north by using 1,2,3, n-1 and n;
4) sequentially acquiring human body positions through a mobile phone positioning system according to equal sampling time intervals delta T, wherein the delta T is T i+1 -t i Where i is 1,2,3, …, m, i is the human body position sampling point number, t i The moment when the position of the human body is acquired for the ith time;
5) at t i At the moment, a dichotomy is adopted to search the number j (j belongs to [1, n ]) of the iron tower closest to the position of the human body]) According to the line where the iron tower j is located, two iron tower numbers (marked as j 'and j') adjacent to the iron tower j on the line are further found on the basis;
6) establishing a three-dimensional coordinate system, respectively determining the coordinates of a human body and iron towers with the numbers of j, j 'and j', simulating a direct current line between two base iron towers by using a straight line segment connecting height points of overhead lines on adjacent iron towers, and calculating and determining a direct current transmission line closest to the human body and the closest horizontal distance x between the direct current transmission line and the human body;
7) searching the extra-high voltage direct current transmission project information base to obtain t i The voltage grade, the number of loops and the ground clearance of the conducting wire of the ultra-high voltage direct current transmission line which is closest to the human body at any moment;
8) the signal intensity of a GPS or Beidou satellite and the positioning position of a human body are obtained through a mobile phone positioning system, the moving speed and the acceleration of the human body are determined according to the signal intensity and the positioning position of the human body, and t is judged i The state of the human body position (outdoor, indoor or in the automobile) at the moment, and stores t i These information of the time of day;
9) according to the position state (outdoor, indoor or in the automobile) of the human body determined in the step 8), substituting the horizontal distance from the actual human body to the center of the nearest extra-high voltage direct current line, the ground clearance of the point on the line closest to the human body, the voltage grade of the line and the number of loops into the corresponding positions established in the step 1)An electric field exposure intensity calculation model calculates to obtain a sampling time t i Human body exposure field intensity E i
10) According to D i =E i Δ T calculation to obtain T i To t i+1 Electric field exposure dose D between moments i (ii) a According to
Figure BDA0002861439010000041
Calculating to obtain the human body electric field exposure dose D within the time T T And m is the number of sampling points in the time T.
According to the technical scheme, according to simulation calculation results, a relation model between outdoor, indoor and automobile human body exposure field intensity E and the voltage level U of the extra-high voltage direct current transmission line, the ground clearance H of a conducting wire, the loop number g and the horizontal distance x between the center of the line and a human body is respectively fitted and established, the direct current transmission line closest to the human body and the position state (outdoor, indoor or automobile) of the human body are determined according to mobile phone positioning signals, and public electric field exposure dose is calculated.
The method comprises the following steps of 1), establishing a relation model between human body exposed field intensity and the direct current transmission line, namely, establishing a geometric model after appropriately simplifying actual ultrahigh voltage direct current transmission lines, the earth, buildings, automobiles and human bodies by using professional software, setting model parameters such as voltage levels and loop numbers of the direct current transmission lines, relative dielectric constants and conductivities of the earth, the buildings, the automobiles and the human bodies, and establishing an ultrahigh voltage direct current transmission line-building (or automobile) -earth-human body electromagnetic simulation model. Considering that the distribution of an electromagnetic field is closely related to boundary conditions, the existence of a human body can cause electric field distortion, the human body models which are positioned on the open ground (namely outdoor), the ground central point in a building model (namely indoor) and the bottom central point in an automobile model (namely in an automobile) are respectively arranged at different horizontal distances from a direct current transmission line with different voltage grades, wire ground clearance and wire loop numbers, the field intensity of the human body surface is calculated based on the simulation of an electric field calculation model near a direct current wire on the infinite ground, and corresponding models are established according to the maximum electric field fitting of the human body surface of the human body model at different positions.
The building model simplify to hollow parallelepiped, its top surface and four side thickness homogeneous phases are the same, the bottom surface is the ground, offers a plurality of rectangular holes on each side of building model for simulate open window.
The automobile model constructed after simplification is a hollow parallelepiped, the thicknesses of six faces of the automobile model are the same, and a plurality of rectangular holes or rectangular areas with different dielectric constants and conductivities are arranged on each side face of the automobile model and used for simulating an open or closed glass window (including a front windshield of an automobile).
The human body is simplified into a male human body model with the total height of 1.73m and composed of five parts, namely a head part, a neck part, an upper body part, arms and legs, and the method is specifically referred to the human body size of Chinese adult (GB/T10000) 1988. Wherein, the head is a sphere with the radius of 0.11 m; the neck is a cylinder with the radius of the bottom surface of 0.05m and the height of 0.06m (the lowest distance from the head sphere); the upper body is a cuboid with a cambered surface, the length of the cuboid is 0.32m, the height of the cuboid is 0.55m, and the thickness of the cuboid is 0.23 m; the arm is two cylinders spliced by a cylinder and a bent cylinder, the radius of the bottom surface is 0.045m, and the total height is 0.58 m; the legs are two cylinders with the radius of the bottom surface of 0.07m and the height of 0.9 m.
In step 4), the shorter the sampling interval Δ T is theoretically, the more accurate the measurement result of the exposure dose of the electric field in each interval is. Generally, the minimum moving speed of an adult is about 0.8m/s (in slow walking), the positioning precision provided by a mobile phone positioning system is about 5-10 m, and in order to ensure that the moving distance of a human body (not at rest) in a position sampling time interval Delta T is greater than the positioning precision, the Delta T is not too short, and preferably 10-15 s.
And 7) taking the ground clearance of the closest point on the wire to the human body according to the ground clearance of the wire of the ultrahigh voltage direct current circuit closest to the human body.
The method for determining the position state of the human body (outdoor, indoor or in-car) in step 8) refers to sampling time interval delta t (delta t-t) according to the mobile phone satellite positioning signal k+1 -t k K is 1,2,3, …, l, k is the satellite positioning signal sampling time sequence number, t k For the kth sampling time, the delta T is an integral fraction of the sampling time interval delta T of the human body position, preferably 1/10), and a satellite is sequentially acquired at different times T by a mobile phone positioning system k Signal intensity W of k (signal-to-noise ratio), human body position longitude and latitude, and human body moving speed v k And acceleration a k . Not provided with t k The moment is at the sampling moment t of the human body position i And t i+1 In between, will t i Longitude and latitude of human body position and moving speed v at any moment i And acceleration a i As t k Corresponding data of time of day, wherein v i =s i /ΔT,a i =(v i+1 -v i )/ΔT,s i Is t i To t i+1 The distance the human body moves between moments. Normally, the satellite signal intensity inside and outside the automobile is greater than the critical value W max (W max Signal-to-noise ratio for weak satellite signals, preferably 20dB), satellite signal strength less than or equal to W for indoor and outdoor partial areas (e.g. mountainous areas) max . Furthermore, if t k-1 To t k The interventricular body enters the outdoor from the indoor and has stronger signal t k-1 To t k Absolute value | Δ W of inter-signal strength variation k-1 |
(ΔW k-1 =W k -W k-1 ) Greater than or equal to a critical value | Δ W min |(|ΔW min | is the absolute value of the difference in satellite signal strength between where the indoor and outdoor signals are stronger, preferably 20dB), if t k-1 To t k The interventricular body enters the stronger signal place from the weaker outdoor signal place, | Δ W k-1 I is less than | Δ W min L. the method is used for the preparation of the medicament. If a certain satellite signal W k Always small sudden appearance of W k >W max Or W k Always large and suddenly appearing W k ≤W max The moving speed, acceleration and satellite signal intensity of human body are combined with t k Time (including t) k ) The previous human body position state is judged, and the specific judgment steps are shown in figure 2.
And 9), if the horizontal distance between the actual human body and the nearest direct current transmission line center exceeds 100m, taking the exposure field intensity (which is close to the background value) when the distance is 100 m.
Compared with the prior art, the invention has the advantages that:
the method for determining the public electric field exposure dose near the ultra-high voltage direct current transmission line by using the mobile phone positioning signal can be realized based on an electromagnetic simulation result and by means of the mobile phone positioning signal, can determine the ultra-high voltage direct current transmission line closest to a human body and the position state of the human body (outdoor, indoor or in an automobile) in real time, calculates the actual electric field exposure dose of the human body in an off-line manner, does not need to use a professional electric field dose instrument for actual measurement, is low in cost, is more convenient and faster, and is more suitable for the public near the ultra-high voltage direct current transmission line.
Drawings
FIG. 1 is a flow chart of a method of determining public electric field exposure dose in the vicinity of an extra-high voltage direct current transmission line in an embodiment of the present invention;
FIG. 2 is a schematic diagram of a process for determining the position status (outdoor, indoor, and in-car) of a human body according to an embodiment of the present invention, wherein (A) is a general flow diagram, and (B), (C) and (D) respectively represent a part of the flow diagram in FIG. A;
FIG. 3 is a schematic view of a building model according to an embodiment of the present invention;
FIG. 4 is a schematic representation of a mannequin according to an embodiment of the present invention;
FIG. 5 is a schematic view of a model of an automobile according to an embodiment of the invention;
FIG. 6 is a schematic illustration of a human body positioned within (a) a building and (b) a model automobile, respectively, in an embodiment of the invention;
fig. 7 is a schematic diagram of a certain day of public motion trajectory near the extra-high voltage linear transmission line in the embodiment of the present invention, wherein a right-side solid line trajectory represents the public motion trajectory identified by the mobile phone positioning system in this situation.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described below with reference to the following embodiments and the accompanying drawings. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the described embodiments without inventive step, are within the scope of protection of the invention.
Unless defined otherwise, technical or scientific terms used herein shall have the ordinary meaning as understood by one of ordinary skill in the art to which this invention belongs. The use of "including" or "comprising" and the like in the present invention is intended to mean that the elements or items listed before the word "comprise" or "comprising" and the like, include the elements or items listed after the word and their equivalents, but do not exclude other elements or items. The terms "connected" or "coupled" and the like are not restricted to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "upper", "lower", "left", "right", and the like are used merely to indicate relative positional relationships, and when the absolute position of the object being described is changed, the relative positional relationships may also be changed accordingly.
Examples
The embodiment carries out off-line calculation on the public actual electric field exposure dose in a certain day near the extra-high voltage direct current transmission line.
Referring to fig. 3, in this embodiment, the building model is simplified to a hollow rectangular parallelepiped with a length × width × height of 12m × 10m × 4m, and the thicknesses of the top surface and four side surfaces are 0.2m, and a rectangular hole with a length × width of 1m × 0.8m is provided at the center of the side surface (12m × 10m) of the building model facing the power transmission line, for simulating an open window.
Referring to fig. 4, the human body is simplified to a male human body model with a total height of 1.73m, which is composed of five parts, namely a head part, a neck part, an upper body part, arms and legs, in this embodiment, specifically referring to the human body size of adult people in china (GB/T10000-. Wherein, the head is a sphere with the radius of 0.11 m; the neck is a cylinder with the radius of the bottom surface of 0.05m and the height of 0.06m (the lowest distance from the head sphere); the upper body is a cuboid with a cambered surface, the length of the cuboid is 0.32m, the height of the cuboid is 0.55m, and the thickness of the cuboid is 0.23 m; the arm is two cylinders spliced by a cylinder and a bent cylinder, the radius of the bottom surface is 0.045m, and the total height is 0.58 m; the legs are two cylinders with the radius of the bottom surface of 0.07m and the height of 0.9 m.
Referring to fig. 5, in this embodiment, the automobile model is simplified to a hollow rectangular parallelepiped with a length × width × height of 4m × 1.7m × 1.5m and six surface thicknesses of 0.03m, two side surfaces of the automobile model are provided with 4 rectangular holes with a length × width of 0.6m × 0.4m for simulating an open glass window, and a front portion of the automobile model is provided with 1 rectangular region with a length × width of 1.6m × 0.6m for simulating a windshield.
Referring to fig. 6, the mannequin located at the center point of the bottom surface in the model building in this embodiment faces the window, and the mannequin located at the center point of the bottom surface in the model car faces the front windshield of the car.
Referring to fig. 1, the electric field exposure dose is determined according to a public positioning signal of a certain day near the extra-high voltage direct current transmission line, and the steps are as follows:
s100, after the ultrahigh voltage direct current transmission line, the earth, the building, the automobile and the human body are properly simplified by using Commol Multiphysics software, a geometric model is constructed, model parameters are set, and an ultrahigh voltage direct current transmission line-building (or automobile) -earth-human body electromagnetic simulation model is constructed. The operating parameters of the direct current transmission line, except the voltage level, the ground height of the wire and the number of loops, are shown in table 1. The geometric shape of the building model is shown in FIG. 3, the material is concrete, and the conductivity is 5 × 10 -2 And (5) S/m. The geometric shape of the phantom is shown in FIG. 4, and the medium parameters are shown in Table 2. The geometric shape of the automobile model is shown in figure 5, the automobile body is regarded as a metal conductor, and the front windshield is not conductive at normal temperature. The ground is regarded as an infinite ideal conductor, the ground potential is 0, the terrain is flat, and the ground does not contain media such as trees, water flow and the like.
Table 1 operating parameters of the embodiment of the extra-high voltage dc transmission line
Figure BDA0002861439010000091
Table 2 medium parameters of the phantom in the example
Figure BDA0002861439010000092
Figure BDA0002861439010000101
S200, the building model is positioned on the open ground (namely outdoor), at the center point of the ground in the building model (namely indoor) and at the center point of the bottom surface in the automobile model (namely in the automobile)The human body models are respectively arranged on the ground positions with different horizontal distances x (m) from the direct current transmission line with different voltage grades, wire ground clearance and wire loop numbers, when x is more than or equal to 0<At 20, the position interval d between adjacent manikins is 0.5 m; when x is more than or equal to 20<At 40, d is 1 m; when x is more than or equal to 40<At 60, d is 2 m; when x is more than or equal to 60<At 100, d is 4 m. Simulating and calculating the human body surface field intensity based on an electric field calculation model near an infinite ground direct current wire, respectively fitting and establishing relational models of outdoor, indoor and automobile human body exposed field intensity E, the voltage level U of an extra-high voltage direct current transmission line, the ground height H of the wire, the number g of loops and the horizontal distance x between the center of the wire and the human body according to the human body surface maximum electric fields of the human body model at different positions, wherein the relational models are respectively E out =f out (U,H,g,x)、E in-b =f in-b (U, H, g, x) and E in-c =f in-c (U, H, g, x); in the formula, E out Exposing field intensity for outdoor human body; e in-b Exposing the field intensity for the indoor human body; e in-c Exposing the field intensity for human body in the automobile.
S300, acquiring all domestic parameters such as voltage classes, loop numbers and the like of the ultra-high voltage direct current transmission lines, information such as positions of all iron towers on the lines and height of the overhead lines at the iron towers, and establishing an ultra-high voltage direct current transmission project information base.
S400, numbering the iron tower from south to north by using 1,2,3, n-1 and n in sequence according to the position of the iron tower of the extra-high voltage direct current transmission line.
S500, sequentially acquiring the positions of the human body through a mobile phone positioning system according to equal sampling time intervals delta T (taking 10S), wherein the delta T is T i+1 -t i Where i is 1,2,3, …, m, i is the sampling point number of the human body position, t i M is the time T for calculating the exposure dose of the electric field at the moment of acquiring the position of the human body for the ith time t The number of sampling points of the inner human body position. The smooth curves are sequentially connected with the human body positions at the sampling moments to obtain the human body motion trail (see fig. 7), and the result shows that after the public starts along the point A on a certain day in the embodiment, the public returns to the point A after passing through the point B, C, D, E in sequence.
S600, at each sampling time t i The number j of the iron tower closest to the position of the human body is searched by adopting a dichotomy (j belongs to [1,n]) And further finding two iron tower numbers (marked as j 'and j') adjacent to the iron tower j on the line on the basis of the line where the iron tower j is located.
S700, establishing a three-dimensional coordinate system, respectively determining coordinates of the human body and iron towers numbered as j, j 'and j', simulating a direct current line between two basic iron towers by using a straight line segment connecting height points of overhead lines on adjacent iron towers, and calculating the horizontal distance d between the human body and the line between the iron towers j and j j&j’ And the horizontal distance d from the line between the iron towers j and j ″ j&j” Setting the result to obtain d j&j’ <d j&j” And d is j&j’ <100m, so the direct current transmission lines between iron towers with the serial numbers of j and j' are the nearest direct current transmission lines within the range of 100m from the human body, and the nearest horizontal distance between the direct current transmission lines and the human body is d j&j’
S800, searching an extra-high voltage direct current transmission project information base, and setting t obtained in the embodiment i The ultra-high voltage direct current line (namely the direct current transmission line where the iron towers j and j' are) closest to the human body at the moment is the +/-1100 kV double-loop direct current transmission line, and the ground clearance of the lead (namely the ground clearance of the point on the lead closest to the human body) is 26 m.
S900, sampling time t of different human body positions in the embodiment based on mobile phone positioning signals i The human body position state (outdoor, indoor, or in the automobile) is judged, and the judgment method is shown in fig. 2. According to the sampling time interval delta t (1 s) of the positioning signal of the mobile phone satellite, the mobile phone positioning system sequentially obtains the satellite at different moments t k Signal strength W of (k 1,2,3, … l, l 10m) k (signal-to-noise ratio), human body position longitude and latitude, and human body moving speed v k And acceleration a k . Without setting t k The moment is at the sampling moment t of the human body position i And t i+1 In between, will t i Human body position longitude and latitude and moving speed v at any moment i And acceleration a i As t k Corresponding data of time of day, wherein v i =s i /ΔT,a i =(v i+1 -v i )/ΔT,s i Is t i To t i+1 The human body movement distance between moments). In general, the interior and the compartment of a carThe intensity of the external satellite signal is greater than the critical value W max (W max 20dB in the embodiment) for the signal-to-noise ratio when the satellite signal is weak, the satellite signal intensity in the indoor and outdoor partial areas (such as mountainous areas) is less than or equal to W max . Furthermore, if t k-1 To t k The interventricular body enters the outdoor from the indoor and has stronger signal t k-1 To t k Absolute value | Δ W of inter-signal strength variation k-1 |(ΔW k-1 =W k -W k-1 ) Greater than or equal to a critical value | Δ W min |(|ΔW min I is the absolute value of the difference in satellite signal strength between where the indoor and outdoor signals are stronger, 20dB in this example), if t k-1 To t k The interventricular body enters the stronger signal place from the weaker outdoor signal place, | Δ W k-1 I is less than | Δ W min |。
After the public starts from A, the satellite signal intensity is always smaller and is at t k1 (k 1 W suddenly appears at the moment of belonging to (1, l)) k1 >W max The public is now at point B. From t k1 Firstly, the human body position state at each moment is sequentially acquired according to the step length-delta t, and in the embodiment, the t is p (p<k 1 ) At time, the human body position state is known as indoor, at t p The human body position state at the next moment is obtained when the moment stops, the public is positioned at the point A, and the movement from the point A to the point B (namely t) can be determined p To t k1 M (not containing t) k1 ) When the public is located indoors.
Further determination of t k1 The human body position state at all times. From t k1 Firstly, the human body acceleration at each moment is sequentially obtained according to the step length + delta t until a certain moment t q (k 1 < q) the absolute value of the acceleration of the human body is greater than or equal to the absolute value of the acceleration | a when the automobile is started or braked for the first time max [ in this example, 2.7m/s is taken 2 ),t k1 To t q The moving speed of the human body is not 0, and the point B (namely t) can be determined k1 Time of day) the public is located outdoors.
After the public starts from the point B, the satellite signal strength is always larger and is at t k2 (k 2 ∈(k 1 L)) suddenly appear W k2 <W max The public is now at point E. t is t k2-1 To t k2 Absolute value | Δ W of satellite signal strength variation k2-1 |(ΔW k2-1 =W k2 -W k2-1 ) Greater than the absolute value of the difference in satellite signal strength between the indoor and outdoor signal strengths | Δ W min I (20 dB in this example), it can be determined at point E (i.e., t) k2 Time of day) the public is located indoors.
Further determination of t k2 Before time (without t) k2 ) The position state of the human body. From t k2 Firstly, the human body position state at each moment is acquired in sequence according to the step length-delta t, and at t k1 (k 1 <k 2 ) The state of the human body position is known as outdoor at time t k1 And stopping acquiring the human body position state at the next moment, wherein the public is positioned at the point B.
Further determination of t k1 To t k2 The position state of the human body. From t k1 To t k2 Sequentially acquiring the acceleration of the human body at each moment according to the step length + delta t at t q (k 1 <q<t k2 ) The absolute value of the acceleration of the human body at the moment is greater than or equal to | a max At t | q Stopping acquiring the acceleration of the human body at the next moment, wherein a q > 0, public at C, and t q-1 A certain time does not exist before, and the time is enabled to reach t q-1 The speed of the human body is continuously 0, and the movement from the point B to the point C (namely t) can be determined k1 To t q-1 Time) the public is located outdoors, t q The public is located in the vehicle at the moment.
Further determination of t q To t k2 The position state of the human body. From t q To t k2 Sequentially acquiring the acceleration of the human body at each moment according to the step length + delta t at t r (q<r<t k2 ) The absolute value of the acceleration of the human body at the moment is more than or equal to | a max At t | r Stopping acquiring the acceleration of the human body at the next moment, wherein a r < 0, the public is at D, and it can be determined to move from C to D (i.e. t) q To t r Time) the public is located in the car.
Further determination of t r To t k2 State of inter-body position. From t r To t k2 Sequentially acquiring the human body acceleration at each moment according to the step length + delta t, wherein the human body acceleration at each moment is less than | a | max L, and t r At a later time, let t r Until the time when the human body speed is continuously 0, the movement from the point D to the point E (namely t) can be determined r To t k2 (not containing t) k2 ) When the public is located outdoors.
The satellite signal intensity of the public moving from the point E to the point A is always small, and the satellite signal intensity is not suddenly larger than W max In this case, it is not necessary to determine the position state of the human body in the time zone. When t is k2 The satellite signal intensity is suddenly greater than W after the moment max Then, the electric field exposure dose of the public during this period is calculated again.
And S1000, calculating the electric field exposure dose of the public outdoors. When the public moves from B to C and D to E, the public is outdoors, and the sampling times of the human body position are M 1 Taking a value less than [ (q-1-k) 1 )+(k 2 -1-r)+2]Maximum integer of/10, sampling time t of public at human body position i1 (1≤i 1 ≤M 1 ) Time-distance closest extra-high voltage direct current transmission line (namely direct current line between iron towers j and j') center horizontal distance x i1 E is substituted by 1100kV line voltage class U, 26m wire ground clearance H and 2 loop number g out =f out (U, H, g, x), calculating to obtain t i1 The exposure field strength E of the moment i1 According to D i1 =E i1 Δ T calculation to obtain T i1 To t i1+1 Dose of electric field exposure D between moments i1 According to
Figure BDA0002861439010000141
Calculating to obtain the electric field exposure dose D of the public outdoors out
And S1100, calculating the electric field exposure dose of the public in the automobile. When the public moves from the point C to the point D, the public is in the vehicle, and the sampling times of the human body position are M 2 Taking the maximum integer less than (r-q +1)/10, and sampling the public at the human body position at the time t i2 (1≤i 2 ≤M 2 ) Time-distance minimum extra-high voltage direct current line (namely direct current line between iron towers j and j') center horizontal distance x i2 Substituting the line voltage class U of 1100kV, the ground height H of the lead of 26m and the loop number g of 2 into E in-c =f in-c (U, H, g, x) to calculate t i2 Exposed field strength E at time i2 According to D i2 =E i2 Δ T calculation to get T i2 To t i2+1 Dose of electric field exposure D between moments i2 According to
Figure BDA0002861439010000142
Calculating to obtain the electric field exposure dose D of the public in the automobile in-c
And S1200, calculating the electric field exposure dose of the public in the room. When the public moves from the point A to the point B (the embodiment does not consider the electric field exposure dose when the public moves from the point E to the point A), the public stays indoors for a time T in-b =T t -(t k2 -t k1 ) Wherein T is t In order to calculate the time of public electric field exposure dose in this embodiment, the central horizontal distance x between the central point in the public housing and the nearest extra-high voltage direct current line (i.e. direct current line between iron towers j and j') is set in-b Substituting the line voltage class U of 1100kV, the ground height H of the lead of 26m and the loop number g of 2 into E in-b =f in-b (U, H, g, x) calculating the exposure field intensity E of the public in the room in-b According to D in-b =E in-b ·T in-b Calculating to obtain the electric field exposure dose D of the public in the room in-b

Claims (8)

1. A method for determining public electric field exposure dose near an extra-high voltage direct current transmission line is characterized by comprising the following steps:
1) establishing an extra-high voltage direct current transmission line-building or vehicle-ground-human body electromagnetic simulation model, respectively establishing a relationship model between the human body exposed field intensity E of outdoor, indoor and automobile and the voltage grade U of the extra-high voltage direct current transmission line, the ground clearance H of a lead, the loop number g and the horizontal distance x between the center of the line and the human body through simulation calculation, wherein the relationship model is respectively E out =f out (U,H,g,x)、E in-b =f in-b (U, H, g, x) and E in-c =f in-c (U,H,g,x);E out Exposing field intensity for outdoor human body; e in-b Exposing the field intensity for the indoor human body; e in-c Exposing the field intensity for human bodies in the automobile;
2) acquiring voltage grades and loop numbers of all domestic ultrahigh voltage direct current transmission lines, positions of all iron towers on a line and stringing height information of the iron towers, and establishing an ultrahigh voltage direct current transmission engineering information base;
3) according to the position of an iron tower of the extra-high voltage direct current transmission line, numbering the iron tower from south to north by using 1,2,3, n-1 and n;
4) sequentially acquiring the positions of the human body through a mobile phone positioning system according to equal sampling time intervals delta T, wherein the delta T is T i+1 -t i I is 1,2,3, is, m, i is a human body position sampling point serial number, t i The moment when the position of the human body is acquired for the ith time;
5) at t i At the moment, a dichotomy is adopted to search the number j (j belongs to [1, n ]) of the iron tower closest to the position of the human body]) According to a line where an iron tower j is located, two iron towers adjacent to the iron tower j on the line are further found on the basis, and the two iron towers adjacent to the iron tower j are marked as j 'and j';
6) establishing a three-dimensional coordinate system, respectively determining the coordinates of a human body and iron towers with the numbers of j, j 'and j', simulating a direct current line between two base iron towers by using a straight line segment connecting height points of overhead lines on adjacent iron towers, and calculating and determining a direct current transmission line closest to the human body and the closest horizontal distance x between the direct current transmission line and the human body;
7) searching the extra-high voltage direct current transmission project information base to obtain t i The voltage grade, the number of loops and the ground clearance of the conducting wire of the ultra-high voltage direct current transmission line which is closest to the human body at any moment;
8) the signal intensity of a GPS or Beidou satellite and the positioning position of a human body are obtained through a mobile phone positioning system, the moving speed and the acceleration of the human body are determined according to the signal intensity and the positioning position of the human body, and t is judged i The human body position state at any moment, including outdoors, indoors or in a car, and storing t i Such information of the time of day;
9) according to the position state of the human body determined in the step 8), substituting the horizontal distance from the actual human body to the center of the nearest extra-high voltage direct current line, the ground clearance of the point on the line closest to the human body, the voltage grade of the line and the number of loops into the corresponding electric field exposure intensity calculation model established in the step 1), and calculating to obtain the sampling time t i Human body exposure field strength E i
10) According to D i =E i Δ T calculation to obtain T i To t i+1 Electric field exposure dose D between moments i (ii) a According to
Figure FDA0003727825570000021
Calculating to obtain the exposure dose D of the electric field of the human body within the time T T Wherein m is the number of sampling points in time T.
2. The method for determining the public electric field exposure dose near the extra-high voltage direct current transmission line according to claim 1, wherein in the step 1), a relation model between the human body exposure field intensity E and the extra-high voltage direct current transmission line is established by the following method:
1-1) simplifying actual extra-high voltage direct current transmission lines, the earth, buildings, automobiles and human bodies by using software, constructing a geometric model, setting model parameters, wherein the model parameters comprise the voltage grade and the loop number of the direct current transmission lines, and the relative dielectric constants and the conductivities of the earth, the buildings, the automobiles and the human bodies, and constructing an extra-high voltage direct current transmission line-building or automobile-earth-human body electromagnetic simulation model;
1-2) considering that the distribution of an electromagnetic field is closely related to boundary conditions and the existence of a human body can cause electric field distortion, respectively placing the human body models which are positioned on the open ground, at the central point of the ground in a building model and at the central point of the bottom in an automobile model at different horizontal distances from a direct current transmission line with different voltage grades, ground height of a wire and the number of loops of the wire, and carrying out simulation calculation on the surface field intensity of the human body based on an electric field calculation model near the direct current wire on the infinite ground;
1-3) establishing corresponding models according to the maximum electric field fitting of the human body surface of the human body model at different positions.
3. The method for determining the public electric field exposure dose near the extra-high voltage direct current transmission line according to claim 2, wherein in the step 1-1), the simplified constructed building model is a hollow parallelepiped, the top surface and four side surfaces of the building model are the same in thickness, the bottom surface is the ground, and a plurality of rectangular holes are formed in each side surface of the building model and used for simulating an open window.
4. The method for determining the public electric field exposure dose near the extra-high voltage direct current transmission line according to claim 2, wherein in the step 1-1), the simplified constructed automobile model is a hollow parallelepiped with six surfaces having the same thickness, and a plurality of rectangular holes or rectangular regions with different dielectric constants and conductivities are arranged on each side surface of the automobile model and used for simulating an open or closed glass window.
5. The method for determining the public electric field exposure dose in the vicinity of the extra-high voltage direct current transmission line according to claim 1, wherein in the step 4), the sampling time interval Δ T is 10-15 s.
6. The method for determining the public electric field exposure dose near the extra-high voltage direct current transmission line according to claim 1, wherein in the step 7), the ground clearance of the conductor of the extra-high voltage direct current transmission line closest to the human body is taken as the ground clearance of the point on the conductor closest to the human body.
7. The method for determining the public electric field exposure dose in the vicinity of an extra-high voltage direct current transmission line according to claim 1, wherein in the step 8), the human body position state determination method is based on a sampling time interval Δ t of a mobile phone satellite positioning signal, where Δ t is t k+1 -t k K is a satellite positioning signal sampling time sequence number, t k At the kth sampling moment, delta T is an integer fraction of the sampling time interval delta T of the human body position, and a mobile phone positioning system is usedThe system sequentially acquires the satellite at different moments t k Signal intensity W of k Longitude and latitude of human body position and human body moving speed v k And acceleration a k
Let t k The moment is at the sampling moment t of the human body position i And t i+1 In between, will t i Human body position longitude and latitude and moving speed v at any moment i And acceleration a i As t k Corresponding data of time of day, wherein v i =s i /ΔT,a i =(v i+1 -v i )/ΔT,s i Is t i To t i+1 The human body movement distance between moments;
if a certain value is always less than W max Of the satellite signal W k Suddenly appear W k >W max Or some is always greater than W max Of the satellite signal W k Suddenly appear W k ≤W max ,W max The signal-to-noise ratio when the satellite signal is weak is combined with the human body movement speed, the acceleration or the satellite signal intensity to t k And judging the position state of the human body before the moment.
8. The method for determining the public electric field exposure dose near the extra-high voltage direct current transmission line according to claim 1, wherein in the step 9), the human body exposure field strength is taken as the exposure field strength when the distance between the actual human body and the center of the nearest direct current transmission line exceeds 100 m.
CN202011564338.1A 2020-12-25 2020-12-25 Method for determining public electric field exposure dose near extra-high voltage direct current transmission line Active CN112765773B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011564338.1A CN112765773B (en) 2020-12-25 2020-12-25 Method for determining public electric field exposure dose near extra-high voltage direct current transmission line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011564338.1A CN112765773B (en) 2020-12-25 2020-12-25 Method for determining public electric field exposure dose near extra-high voltage direct current transmission line

Publications (2)

Publication Number Publication Date
CN112765773A CN112765773A (en) 2021-05-07
CN112765773B true CN112765773B (en) 2022-08-23

Family

ID=75694595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011564338.1A Active CN112765773B (en) 2020-12-25 2020-12-25 Method for determining public electric field exposure dose near extra-high voltage direct current transmission line

Country Status (1)

Country Link
CN (1) CN112765773B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435873A (en) * 2011-09-09 2012-05-02 国网电力科学研究院 Method for determining protection space between ultra-high-voltage DC (direct current) transmission line and wireless station

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001343411A (en) * 2000-03-29 2001-12-14 Sony Corp Forecasting apparatus for intensity of radiation electric field, method thereof and distributing medium
CN101221202A (en) * 2007-12-20 2008-07-16 国网武汉高压研究院 Method for appraising electromagnetic environment level of 1000KV ultrahigh voltage alternating current transmission lines
US9265331B1 (en) * 2010-01-21 2016-02-23 Hd Electric Company Wearable carrying case
CN103344810A (en) * 2013-06-25 2013-10-09 国家电网公司 Method for calculating electric shock effects of vehicle on human body under high-voltage alternating current transmission line
KR102078988B1 (en) * 2013-09-30 2020-02-19 한국전력공사 Apparatus for measuring dc electric field
CN104814805B (en) * 2015-04-15 2016-08-24 浙江大学 A kind of cage for rearing poultry of animal model under extra-high voltage direct-current transmission line
CN105260583B (en) * 2015-11-27 2020-07-24 国网重庆市电力公司电力科学研究院 Method and system for calculating biological effect of ultra-high voltage power frequency electromagnetic field on human body
FR3071652B1 (en) * 2017-09-26 2019-09-27 Ineo Hauts-De-France DEVICE FORMING INDIVIDUAL EQUIPMENT FOR PROTECTING PEOPLE EVOLVING IN AN AREA OF ELECTRICAL RISK
CN108345743B (en) * 2018-02-09 2021-08-03 武汉大学 Electromagnetic safety analysis method for wireless power supply device of high-voltage line online monitoring terminal
CN109374991A (en) * 2018-07-31 2019-02-22 国网江苏省电力有限公司徐州供电分公司 A kind of ultra-high-tension power transmission line voltage class confirmation device and its confirmation method based on the measurement of different location electric field change
CN109918813A (en) * 2019-03-14 2019-06-21 国网上海市电力公司 A kind of power frequency electric field acquisition methods based on around ultra-high-tension power transmission line
CN110516309B (en) * 2019-07-31 2022-10-11 国网福建省电力有限公司建设分公司 Method for predicting and controlling power frequency electric field of roof near same-tower double-circuit 500kV power transmission line

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435873A (en) * 2011-09-09 2012-05-02 国网电力科学研究院 Method for determining protection space between ultra-high-voltage DC (direct current) transmission line and wireless station

Also Published As

Publication number Publication date
CN112765773A (en) 2021-05-07

Similar Documents

Publication Publication Date Title
CN109413578B (en) Indoor positioning method based on fusion of WIFI and PDR
Zuo et al. Multi-phase fingerprint map based on interpolation for indoor localization using iBeacons
Li et al. How feasible is the use of magnetic field alone for indoor positioning?
CN105717483B (en) A kind of location determining method and device based on multi-source positioning method
CN106162555A (en) Indoor orientation method and system
CN113923596B (en) Indoor positioning method, device, equipment and medium
CN103167606B (en) Based on the WLAN indoor orientation method of rarefaction representation
CN109945865A (en) The indoor orientation method merged based on WiFi with earth magnetism
CN109218983A (en) Positioning method and positioning system
CN108801267B (en) Indoor anchor-free positioning method integrating multiple sensors
CN112954594A (en) Wireless sensor network node positioning algorithm based on artificial bee colony
CN109188360A (en) A kind of indoor visible light 3-D positioning method based on bat algorithm
CN108738132A (en) A kind of three base station movement communication positioning methods based on TDOA
Wattananavin et al. A comparative study of RSSI-based localization methods: RSSI variation caused by human presence and movement
CN106358233B (en) A kind of RSS data smoothing method based on Multidimensional Scaling algorithm
CN106197418B (en) A kind of indoor orientation method merged based on the fingerprint technique of sliding window with sensor
CN112765773B (en) Method for determining public electric field exposure dose near extra-high voltage direct current transmission line
CN105044659A (en) Indoor positioning device based on environment spectrum fingerprint
CN106658538A (en) Mobile phone base station signal coverage range simulation method based on Thiessen polygon
CN107105404B (en) Pedestrian indoor positioning method based on step length matching
Al-Sabbagh et al. A centralized multi-floor indoor navigation system for a large mall
CN109640251B (en) Indoor positioning method and device
CN114173281B (en) TOA-based positioning system beacon node optimization layout method in indoor NLOS environment
CN110187306A (en) A kind of TDOA-PDR-MAP fusion and positioning method applied to the complicated interior space
Jia et al. An indoor and outdoor seamless positioning system based on android platform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant