CN111355545A - Simplified directional function-based base station electromagnetic radiation visual simulation prediction method - Google Patents

Simplified directional function-based base station electromagnetic radiation visual simulation prediction method Download PDF

Info

Publication number
CN111355545A
CN111355545A CN202010165360.2A CN202010165360A CN111355545A CN 111355545 A CN111355545 A CN 111355545A CN 202010165360 A CN202010165360 A CN 202010165360A CN 111355545 A CN111355545 A CN 111355545A
Authority
CN
China
Prior art keywords
antenna
electromagnetic radiation
base station
power
angle
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.)
Granted
Application number
CN202010165360.2A
Other languages
Chinese (zh)
Other versions
CN111355545B (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.)
Xiamen Environmental Monitoring Center Station Of Fujian Province Jiulongjiang River Basin Ecological Environment Monitoring Center
Original Assignee
Xiamen Environmental Monitoring Center Station Of Fujian Province Jiulongjiang River Basin Ecological Environment Monitoring Center
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 Xiamen Environmental Monitoring Center Station Of Fujian Province Jiulongjiang River Basin Ecological Environment Monitoring Center filed Critical Xiamen Environmental Monitoring Center Station Of Fujian Province Jiulongjiang River Basin Ecological Environment Monitoring Center
Priority to CN202010165360.2A priority Critical patent/CN111355545B/en
Publication of CN111355545A publication Critical patent/CN111355545A/en
Application granted granted Critical
Publication of CN111355545B publication Critical patent/CN111355545B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3913—Predictive models, e.g. based on neural network models
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/373—Predicting channel quality or other radio frequency [RF] parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses a base station electromagnetic radiation visual simulation prediction method based on a simplified directional function, which comprises the following steps of; the practical simplified calculation of any point in space is realized by mathematically simplifying the directional attenuation function of the antenna; predicting each equivalent position curve of the antenna space limit value to realize the prediction of the single-antenna three-dimensional visual equivalent surface; meanwhile, the prediction of a spatial electromagnetic radiation power density visualization curve in the vertical direction and the prediction of a ground 1.7m high electromagnetic radiation power density visualization curve are realized. According to the invention, through open design and conventional multi-parameter adjustment, the established prediction model can accurately calculate and visualize the maximum electromagnetic radiation intensity at a certain position of a single base station, so that the method has great reference value for quick site selection of the base station, visualization evaluation of electromagnetic radiation environment influence and visualization of environmental protection, and has certain social benefit.

Description

Simplified directional function-based base station electromagnetic radiation visual simulation prediction method
Technical Field
The invention relates to the field of base station wire radiation, in particular to a base station electromagnetic radiation visual simulation prediction method based on a simplified directivity function.
Background
At present, aiming at the electromagnetic radiation prediction around a communication base station, the calculation of the directional attenuation function of an antenna is complicated, has no practicability, and has no method for simply realizing visualization and effectiveness. In fact, practical simplified calculation of any point in space is realized by simplifying the directional attenuation function of the antenna mathematically, and the method is the basis of space calculation, three-dimensional display, simplicity and easiness in use of the communication base station electromagnetic radiation visual simulation prediction system. If the maximum transmitting power is adopted to predict the electromagnetic radiation around the base station, the deviation between the predicted value and the actual measured value can be corrected through the actual evaluation of the network load. And three-dimensional simulation visualization can embody the range of the electromagnetic radiation limit value of a single-antenna space, and show the superimposed electromagnetic radiation power density curve of the space and the ground multi-antenna, so that the influence evaluation of the electromagnetic radiation environment around the base station can be quickly made only by visual inspection. The invention simplifies the directional attenuation function of the antenna through mathematics, combines with network load correction, realizes simpler and more accurate electromagnetic radiation prediction and realizes the visualization of the superposition effect.
Disclosure of Invention
The invention mainly aims to overcome the defects in the prior art and provides a practical base station electromagnetic radiation visual simulation prediction method based on a simplified directional function.
The invention adopts the following technical scheme:
s1, simplifying a directional attenuation function according to the characteristic that the angle corresponding to half of the attenuation of the main lobe intensity of the antenna radiation unit wave beam is a half-power angle and the using frequency of the antenna, the level of an antenna directional diagram and the half-power angle;
s2, calculating the transmitting power P of the antenna by superposing network load according to the antenna parameters through a loss model;
s3, calculating equivalent position distances r corresponding to space power density limit values of the corresponding vertical horizontal direction main lobe center position, half power angle bilateral symmetry position and double half power angle bilateral symmetry position in groups according to the antenna radiation main lobe center position, 5 vertical sections of the horizontal direction position, the half power angle bilateral symmetry position and the double half power angle bilateral symmetry position according to the antenna radiation main lobe center position, the half power angle bilateral symmetry position and the double half power angle bilateral symmetry position by combining the directional attenuation function;
s4: repeating S3, forming a position curve on three-dimensional software according to each equivalent position distance r obtained by 5 vertical sections of the main lobe horizontal direction position, the half-power angle bilateral symmetry position and the double half-power angle bilateral symmetry position, and realizing three-dimensional visual simulation prediction of the single-antenna electromagnetic radiation power density isosurface model space;
s5: through the said dayLinear transmitting power P, and calculating electromagnetic radiation power density P of a certain point in space in the vertical direction of each antenna in the same direction by combining the directional attenuation function, the horizontal distance between the mobile terminal and the base station, the downward inclination angle between the mobile terminal and the opposite antenna of the base station and the height differenced1；
S6: and repeating S5, forming a vertical electromagnetic radiation intensity curve on the three-dimensional software, and realizing visual simulation prediction of the spatial electromagnetic radiation power density curve in the vertical direction of the multi-antenna in the same direction.
S7: calculating the power density P of the high-level equal-distance electromagnetic radiation of 1.7m on the ground by the energy addition through the antenna transmitting power P and combining the corrected directional attenuation function, superposing according to the ground mirror reflection power density, not considering the phase difference of the direct wave and the ground reflection waved2；
S8: and repeating S7, forming a horizontal electromagnetic radiation intensity curve on the three-dimensional software, and realizing the visual simulation prediction of the high-level equal-distance electromagnetic radiation power density curve of 1.7m on the ground in the same-direction multi-antenna horizontal direction.
Specifically, the simplified mathematical expression of the directional attenuation function is:
Figure BDA0002407256140000021
wherein:
Figure BDA0002407256140000022
represents the horizontal tilt angle of the main lobe horizontal plane; x represents a directional attenuation factor corresponding to a horizontal half-power angle; θ represents a declination angle of the main lobe vertical plane; y represents a directional attenuation factor corresponding to the vertical half-power angle.
Specifically, the simplified expression of the antenna's inherent vertical directivity attenuation function is:
cosyθ1/2＝0.5；y＝lg(0.5)/lg(cosθ1/2)
a simplified expression of the horizontal directivity attenuation function inherent to the antenna is:
Figure BDA0002407256140000023
wherein: 2 theta1/2Representing the inherent vertical half-power angle of the main lobe of the antenna radiation, rounded by y,
Figure BDA0002407256140000024
representing the inherent horizontal half-power angle of the antenna radiating main lobe, x is rounded.
E.g. y 4, when the antenna vertical half power angle 2 theta1/2＝65°；
y is 90, and the antenna vertical half power angle 2 theta1/2＝14°；
y 370, and antenna vertical half power angle 2 θ1/2＝7°；
Similarly, x is 4, and the antenna horizontal half power angle
Figure BDA0002407256140000031
x is 90, and antenna horizontal half power angle
Figure BDA0002407256140000032
370, antenna horizontal half power angle
Figure BDA0002407256140000033
Specifically, in step S2, the expression of the antenna transmission power is:
P＝K×Pt×G
in the formula: p represents the antenna transmitting power and has the unit of W; ptRepresents the transmit power fed into the antenna port in units of W; g represents antenna gain in multiples; k represents the superposition of the various attenuations.
Specifically, in step S3, the expression of each equivalent position distance r corresponding to the spatial power density limit value is:
Figure BDA0002407256140000034
wherein: p represents the antenna transmission power in W, Pd1Representing a given power density limit in W/m2，
Figure BDA0002407256140000035
Is a directional attenuation function.
Specifically, in step S5, the electromagnetic radiation power density expression at a certain point in the space in the vertical direction of a certain antenna in the same direction is as follows:
Figure BDA0002407256140000036
Figure BDA0002407256140000041
wherein: pd1The electromagnetic radiation power density of a certain point of a communication base station is represented, and P represents the unit of the antenna transmission power as W/m2，siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
Specifically, in step S7, the expression for calculating the maximum electromagnetic radiation power density at a high level and equal distance of 1.7m on the ground is as follows:
Figure BDA0002407256140000042
the modified directional attenuation function
Figure BDA0002407256140000043
Figure BDA0002407256140000044
Figure BDA0002407256140000045
The normalized directivity basis function is:
Figure BDA0002407256140000046
Figure BDA0002407256140000047
wherein: r is2The unit of the propagation distance from the position of the radiation lobe to the center point to the prediction point is m; theta2To predict the ground projection angle of the point reflection, P represents the antenna transmission power in W/m2；siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
As can be seen from the above description of the present invention, compared with the prior art, the present invention has the following advantages:
according to the invention, through open design, conventional multi-parameter adjustment, mathematical simplification of an antenna directional attenuation function and network load correction, the established prediction model can accurately calculate and visualize the maximum electromagnetic radiation intensity at a certain position of a single base station, and has great reference value for quick site selection of the base station, visualization evaluation of electromagnetic radiation environment influence and visualization of environment protection, and certain social benefit is achieved.
Drawings
FIG. 1 is a flow chart of the present invention
FIG. 2 is a schematic diagram of a visualization result of evaluation of electromagnetic radiation environment influence of a base station;
FIG. 3 is a diagram of the predicted maximum value of high power density of 1.7m on the ground of a certain base station, and the average value and instantaneous maximum value of the measured 6 minutes in multiple days:
FIG. 3(a) is the power density calculation and measured value of a 1.7m antenna on the ground for a layer 1 GSM-900MHz antenna of a base station;
FIG. 3(b) is the power density calculation and measured value of a 1.7m antenna on the ground for a base station layer 2 TDD-2600MHz antenna;
FIG. 3(c) is the calculation and measured values of the ground power density of 1.7m for a base station layer 3 FDD-1870MHz antenna;
FIG. 3(d) is the calculated and measured values of the total power density of the 1.7m antenna on the ground for the layer 3 FDD-875MHz antenna of a certain base station;
FIG. 3(e) is the calculated and measured values of the total power density of the antennas of 1.7m on the ground of a certain base station;
FIG. 4 is a comparison graph of the mathematical simplified antenna directivity attenuation function ground calculation of the present invention:
fig. 4(a) is a directional attenuation diagram of a certain fixed parameter antenna, height difference 20m, downtilt angle 6 °, and vertical half-power angle 65 °;
fig. 4(b) is a directional attenuation diagram of a certain fixed parameter antenna, height difference 20m, downtilt angle 6 °, and vertical half-power angle 30 °;
fig. 4(c) is a directional attenuation diagram of a certain fixed parameter antenna, height difference 20m, downtilt angle 6 °, and vertical half-power angle 14 °;
FIG. 4(d) is a directional attenuation plot for a fixed parameter antenna, height difference 20m, downtilt angle 6 °, vertical half-power angle 7 °;
Detailed Description
The invention is further described below by means of specific embodiments.
The ground is actually measured and verified only by adopting multiple groups of ground for multiple days, the experimental ground implemented by the group is near a No. 1 pole-holding base station, the main radiation direction is in an open and flat area along the road, the obstacles on the road are few, the measurement is started from a position 5 meters away from the base station, and the horizontal movement is started at a measurement interval of 5 meters along the maximum radiation direction of the base station to perform fixed-point measurement to 200 m. The test method comprises the following steps: according to standard HJ/T10.2-1996 electromagnetic radiation monitoring instrument and method for radiation environment protection management guide and mobile communication base station electromagnetic radiation environment monitoring method (HJ 972-2018), the power density value and the maximum value of each measuring point are averagely tested according to 6 minutes, wherein the distance between each measuring point and the ground is 1.7 m. Testing equipment: the portable electromagnetic radiation spectrum analyzer SRM3006, the probe K-1054(27MHz-3GHz) and the probe E-0141(420MHz-6 GHz). And 3 layers of base stations have different frequency standards, the single-antenna double-frequency of the 3 rd layer is calculated according to double antennas, and each layer of 3 horizontal half-power angle 65-degree directional antennas sector covers a circular area with the radius of 500 meters. The parameters of the 1-3 layers of antennas are respectively antenna types GSM, TD-LTE and FDD-LTE, typical integrated loss attenuation is-4.7 dBi, -3dBi, -4.7dBi and-4.7 dBi, transmitting power is 40W, 10W, 20W and 20W, carrier number is 3, 1, 2 and 2, channel number is 1, 8 and 2, network load is preset for 100%, 100% and 100%, channel resource conversion rate is 100%, 50% and 100%, power control is 100%, 100% and 100%, antenna energy efficiency is 80%, 80% and 80%, service beam antenna gain is 10dBi, 19.5dBi and 24dBi &17.5dBi, antenna vertical half-power angle is 65 °, 14 °, 7 ° and 14 °, total downward inclination angle is 6 °, 11 °, 8 ° and 10 ° and test height difference is 18m, 16m, 15m and 15m &15 m.
Fig. 1 is a flow chart of the present invention, and the specific steps are as follows:
s1, simplifying a directional attenuation function according to the characteristic that the angle corresponding to half of the attenuation of the main lobe intensity of the antenna radiation unit wave beam is a half-power angle and the using frequency of the antenna, the level of an antenna directional diagram and the half-power angle;
specifically, the simplified mathematical expression of the directional attenuation function is:
Figure BDA0002407256140000071
wherein:
Figure BDA0002407256140000072
represents the horizontal tilt angle of the main lobe horizontal plane; x represents a directional attenuation factor corresponding to a horizontal half-power angle; θ represents a declination angle of the main lobe vertical plane; y represents a directional attenuation factor corresponding to the vertical half-power angle.
Specifically, the simplified expression of the antenna's inherent vertical directivity attenuation function is:
cosyθ1/2＝0.5；y＝lg(0.5)/lg(cosθ1/2)
a simplified expression of the horizontal directivity attenuation function inherent to the antenna is:
Figure BDA0002407256140000073
wherein: 2 theta1/2Representing the inherent vertical half-power angle of the main lobe of the antenna radiation, rounded by y,
Figure BDA0002407256140000074
watch (A)
The inherent horizontal half-power angle of the main lobe of the antenna radiation is shown, and x is rounded.
E.g. y 4, when the antenna vertical half power angle 2 theta1/2＝65°；
y is 90, and the antenna vertical half power angle 2 theta1/2＝14°；
y 370, and antenna vertical half power angle 2 θ1/2＝7°；
Similarly, x is 4, and the antenna horizontal half power angle
Figure BDA0002407256140000075
x is 90, and antenna horizontal half power angle
Figure BDA0002407256140000076
370, antenna horizontal half power angle
Figure BDA0002407256140000077
In this embodiment, the y values of the 1-3 layers of antennas are 4, 90, 370, and 90, respectively; in this embodiment, according to the requirement of environmental protection, only the maximum value in the vertical direction is obtained, and the simplified mathematical expression is expressed as cos4θ，cos90θ，cos370θ，cos90θ；
S2: calculating the transmitting power P of the antenna by a loss model according to the antenna parameters and superposed network loads;
specifically, in step S2, the expression of the antenna transmission power is:
P＝K×Pt×G
in the formula: p represents the antenna transmission powerRate, in units of W; ptRepresents the transmit power fed into the antenna port in units of W; g denotes the antenna gain in multiples.
And according to PtConventional calculation formulas can be found:
calculating the layer 1 antenna P value as: p10(-4.7/10)×(40×3×1×1×1×1×0.8)×10(10/10)＝325.3
Calculating the layer 2 antenna P value as: p10(-3/10)×(10×1×8×1×0.5×1×0.8)×10(19.5/10)＝1429.4
Calculating the layer 3 antenna frequency-P value as: p10(-4.7/10)×(20×2×2×1×1×1×0.8)×10(24/10)＝5447.3
Calculating the second layer 3 antenna frequency as: p10(-4.7/10)×(20×2×2×1×1×1×0.8)×10(17.5/10)＝1219.5
S3: calculating equivalent position distances r corresponding to space power density limits of corresponding vertical horizontal direction positions, half-power angle bilateral symmetry positions and double half-power angle bilateral symmetry positions in groups according to 5 vertical sections of the antenna radiation main lobe center position, horizontal direction positions, half-power angle bilateral symmetry positions and double half-power angle bilateral symmetry positions through the antenna transmission power P and the directional attenuation function in combination with the national standard electromagnetic radiation power density limit;
specifically, the expression of each equivalent position distance r corresponding to the spatial power density limit value is as follows:
Figure BDA0002407256140000081
wherein: p represents the antenna transmission power in W, Pd1Representing a given power density limit in W/m2，
Figure BDA0002407256140000082
Is a directional attenuation function.
S4: repeating S3, forming a position curve on three-dimensional software according to each equivalent position distance r obtained by 5 vertical sections of the main lobe horizontal direction position, the half-power angle bilateral symmetry position and the double half-power angle bilateral symmetry position, and realizing three-dimensional visual simulation prediction of the single-antenna electromagnetic radiation power density isosurface model space;
and during evaluation, the isosurface is displayed by a spatial three-dimensional limit value isosurface, and the isosurface is visually observed not to touch the standing position of a person, so that the evaluation of the electromagnetic radiation environment influence reaches the standard.
S5: calculating the electromagnetic radiation power density P of a certain point in the space in the vertical direction of each antenna in the same direction by the antenna transmission power P and combining the directional attenuation function, the horizontal distance between the mobile terminal and the base station, the downward inclination angle between the mobile terminal and the opposite antenna of the base station and the height differenced1；
Specifically, the electromagnetic radiation power density expression at a certain point in the space in the vertical direction of a certain antenna in the same direction is as follows:
Figure BDA0002407256140000091
Figure BDA0002407256140000092
wherein: pd1The electromagnetic radiation power density of a certain point of a communication base station is represented, and P represents the unit of the antenna transmission power as W/m2，siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
S6: and repeating S5, forming a vertical electromagnetic radiation intensity curve on the three-dimensional software, and realizing visual simulation prediction of the spatial electromagnetic radiation power density curve in the vertical direction of the multi-antenna in the same direction.
Similarly, after prediction, evaluation is carried out, the intensity maximum curve of the base station from the horizontal to the vertical plane is displayed, and the position of a standard line of 40 muW-cm & lt-2 & gt on the upper part of the curve is not touched by visual observation, so that the electromagnetic radiation environment influence evaluation reaches the standard;
s7: transmitting power P through the antennaAnd calculating the power density P of the electromagnetic radiation with the high level and the equal distance of 1.7m on the ground by combining the corrected directional attenuation functiond2；
Specifically, the expression for calculating the maximum electromagnetic radiation power density at a high level equal distance of 1.7m on the ground is as follows:
Figure BDA0002407256140000093
the modified directional attenuation function
Figure BDA0002407256140000094
Figure BDA0002407256140000095
The normalized directivity basis function is:
Figure BDA0002407256140000101
Figure BDA0002407256140000102
wherein: r is2The unit of the propagation distance from the position of the radiation lobe to the center point to the prediction point is m; theta2To predict the ground projection angle of the point reflection, P represents the antenna transmission power in W/m2；siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
In the correction of the directional function, in view of the environmental requirement characteristics, the influence of the mirror ground on the electromagnetic radiation sight distance propagation is only considered, ① total reflection (the reflection coefficient is 1), ② energy addition is adopted without considering the phase difference between the direct wave and the ground reflected wave, and only the maximum value in the vertical direction is obtained.
In the embodiment of the invention, a point and a baseHorizontal distance of station antenna is siHorizontal distance of hjAnd the main radiation direction is along an open and flat area of the road, the obstacles on the road are fewer, the electromagnetic radiation magnetic radiation power density at the height of 1.7m on the ground is measured from the position 5m away from the base station, and the main radiation direction starts to horizontally move along the maximum radiation direction of the base station at the measurement interval of 5m to perform fixed-point measurement to 200 m.
S8: and repeating S7, forming a horizontal electromagnetic radiation intensity curve on the three-dimensional software, and realizing the visual simulation prediction of the high-level equal-distance electromagnetic radiation power density curve of 1.7m on the ground in the same-direction multi-antenna horizontal direction.
The evaluation method comprises the following steps: the intensity maximum curves at different horizontal distances of 1.7m above the base station floor are shown so as to be visually observed without touching the upper 40 μ W cm of the curve-2And the standard line position is the standard of the electromagnetic radiation environment influence evaluation.
Fig. 2 shows a three-dimensional visualization model for prediction calculation and implementation of the present invention. The middle three-dimensional model in FIG. 2 is calculated by mathematically simplifying the directional attenuation function of the antenna, so as to realize the three-dimensional visual simulation prediction of the isosurface model space of the electromagnetic radiation power density of a plurality of single antennas. The evaluation method comprises the following steps: the equivalent surface is not touched with the standing position of a person by visual observation, and the evaluation of the electromagnetic radiation environment reaches the standard. And (3) adopting a mathematical simplified antenna directivity attenuation function calculation to realize the same-direction multi-antenna space electromagnetic radiation superposition power density intensity curve visual simulation prediction of the left vertical curve model of the figure 2. The evaluation method comprises the following steps: the standard line position of 40 muW cm-2 at the upper part of the curve is not touched by visual inspection, so that the evaluation of the electromagnetic radiation environment influence reaches the standard. The right horizontal curve model in fig. 2 considers the influence of antenna radiation side lobes on the ground, and adopts the directional attenuation function calculation of a uniform linear array to realize the visualized simulation prediction of the high-level equal-distance electromagnetic radiation superposition power density curve of 1.7m on the ground in the same-direction multi-antenna horizontal direction. The evaluation method is that the standard line position of 40 muW cm-2 at the upper part of the curve is not touched by visual inspection, and the evaluation of the electromagnetic radiation environment effect reaches the standard. In contrast to the simplified calculation, fig. 3, 4 can be seen.
FIG. 3 is a graph of two predicted calculated maximum values of ground 1.7m high power density, and a comparison of the measured 6-minute average value and the measured instantaneous maximum value over multiple days, where the measured instantaneous maximum value is close to the predicted maximum value, in accordance with an embodiment of the present invention.
Fig. 4 is a comparison of the mathematically simplified antenna directivity attenuation function calculation of the present invention and the directivity function calculation using a uniform linear array. The deviation except the influence of the side lobe at the near part of the base station is less than 1 mu W cm under the antenna height difference of 20m and different vertical half power angles-2The two predictive calculations for a ground 1.7m high power density almost coincide.
The above description is only an embodiment of the present invention, but the design concept of the present invention is not limited thereto, and any insubstantial modifications made by using the design concept should fall within the scope of infringing the present invention.

Claims (7)

1. The base station electromagnetic radiation visual simulation prediction method based on the simplified directional function comprises the following steps:
s1, simplifying a directional attenuation function according to the characteristic that the angle corresponding to half of the attenuation of the main lobe intensity of the antenna radiation unit wave beam is a half-power angle and the using frequency of the antenna, the level of an antenna directional diagram and the half-power angle;
s2, calculating the transmitting power P of the antenna by superposing network load according to the antenna parameters through a loss model;
s3, calculating equivalent position distances r corresponding to space power density limit values of the corresponding vertical direction main lobe center position, half power angle bilateral symmetry position and double half power angle bilateral symmetry position in groups according to the antenna radiation main lobe center position, 5 vertical sections of the horizontal direction position, the half power angle bilateral symmetry position and the double half power angle bilateral symmetry position according to the national standard electromagnetic radiation power density limit value through the antenna transmission power P and the combination of the directional attenuation function;
s4, repeating S3, and forming a position curve on three-dimensional software according to distances r of each equivalent position obtained by 5 vertical sections of the main lobe horizontal direction position, the half-power angle bilateral symmetry position and the double half-power angle bilateral symmetry position, so as to realize three-dimensional visual simulation prediction of the single-antenna electromagnetic radiation power density isosurface model space;
s5, calculating the electromagnetic radiation power density P of a certain point in the space in the vertical direction of each antenna in the same direction by the transmission power P of the antenna and combining the directional attenuation function, the horizontal distance between the mobile terminal and the base station, the downward inclination angle between the mobile terminal and the opposite antenna of the base station and the height differenced1；
S6, repeating S5, forming a vertical electromagnetic radiation intensity curve on the three-dimensional software, and realizing visual simulation prediction of the spatial electromagnetic radiation power density curve in the same-direction multi-antenna vertical direction;
s7, transmitting power P through the antenna, combining with the corrected directional attenuation function, superposing according to the ground mirror reflection power density, not considering the phase difference of the direct wave and the ground reflection wave, and calculating the high-level equal-distance electromagnetic radiation power density P of 1.7m on the ground by energy additiond2；
And S8, repeating S7, and forming a horizontal electromagnetic radiation intensity curve on the three-dimensional software to realize the visual simulation prediction of the high-level equal-distance electromagnetic radiation power density curve of 1.7m on the ground in the same-direction multi-antenna horizontal direction.
2. The method for base station electromagnetic radiation visualization simulation prediction based on the simplified directivity function as claimed in claim 1, wherein in the step S1, the simplified mathematical expression of the directivity attenuation function is:
Figure FDA0002407256130000021
wherein:
Figure FDA0002407256130000022
represents the horizontal tilt angle of the main lobe horizontal plane; x represents a directional attenuation factor corresponding to a horizontal half-power angle; θ represents a declination angle of the main lobe vertical plane; y represents a directional attenuation factor corresponding to the vertical half-power angle.
3. The method for visualization simulation prediction of electromagnetic radiation of base station based on simplified directivity function as claimed in claim 2, wherein the simplified expression of the antenna inherent vertical directivity attenuation function is:
cosyθ1/2＝0.5；y＝lg(0.5)/lg(cosθ1/2)
a simplified expression of the horizontal directivity attenuation function inherent to the antenna is:
Figure FDA0002407256130000023
wherein 2 theta is1/2The inherent vertical half-power angle of the antenna radiation main lobe is represented, and y is rounded;
Figure FDA0002407256130000024
representing the inherent horizontal half-power angle of the antenna radiating main lobe, x is rounded.
4. The method for visualization, simulation and prediction of electromagnetic radiation of base station based on simplified directivity function as claimed in claim 3, wherein in said step S2, the expression of said antenna transmission power is:
P＝K×Pt×G
p represents the transmitting power of the antenna, and the unit is W; ptRepresents the transmit power fed into the antenna port in units of W; g represents antenna gain in multiples; k represents the superposition of the various attenuations.
5. The method for base station electromagnetic radiation visualization simulation prediction based on the simplified directional function as claimed in claim 4, wherein in step S3, the expression of each equivalent position distance r corresponding to the spatial power density limit is:
Figure FDA0002407256130000025
wherein, the unit is W, P representing the transmitting power of the antennad0Representing a given power density limit in W/m2，
Figure FDA0002407256130000031
Is a directional attenuation function.
6. The method for base station electromagnetic radiation visualization simulation prediction based on the simplified directivity function as claimed in claim 5, wherein in step S5, the electromagnetic radiation power density expression at a certain point in the space in the vertical direction of a certain antenna in the same direction is:
Figure FDA0002407256130000032
Figure FDA0002407256130000033
wherein P isd1The electromagnetic radiation power density of a certain point of a communication base station is represented, and P represents the unit of the antenna transmission power as W/m2，siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
7. The method for base station electromagnetic radiation visualization simulation prediction based on the simplified directivity function as claimed in claim 4, wherein in the step S7, the expression for calculating the maximum electromagnetic radiation power density at a high level equal distance of 1.7m on the ground is:
Figure FDA0002407256130000034
the modified directional attenuation function
Figure FDA0002407256130000035
Figure FDA0002407256130000036
The normalized directivity basis function is:
Figure FDA0002407256130000037
Figure FDA0002407256130000038
wherein r is2The unit of the propagation distance from the position of the radiation lobe to the center point to the prediction point is m; theta2To predict the ground projection angle of the point reflection, P represents the antenna transmission power in W/m2；siRepresenting the horizontal distance from a certain point in space to the base station, and the unit is m; h isjThe height difference of a certain point in space from the base station antenna is shown in m, and r is the distance of the certain point in space from the base station antenna and is shown in m.
CN202010165360.2A 2020-03-11 2020-03-11 Simplified directional function-based base station electromagnetic radiation visual simulation prediction method Active CN111355545B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010165360.2A CN111355545B (en) 2020-03-11 2020-03-11 Simplified directional function-based base station electromagnetic radiation visual simulation prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010165360.2A CN111355545B (en) 2020-03-11 2020-03-11 Simplified directional function-based base station electromagnetic radiation visual simulation prediction method

Publications (2)

Publication Number Publication Date
CN111355545A true CN111355545A (en) 2020-06-30
CN111355545B CN111355545B (en) 2021-12-17

Family

ID=71197312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010165360.2A Active CN111355545B (en) 2020-03-11 2020-03-11 Simplified directional function-based base station electromagnetic radiation visual simulation prediction method

Country Status (1)

Country Link
CN (1) CN111355545B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113613282A (en) * 2021-08-12 2021-11-05 广州广电计量检测股份有限公司 Method and device for simulating radio propagation electromagnetic environment of urban scene base station
CN115078850A (en) * 2022-07-19 2022-09-20 沸蓝建设咨询有限公司 Communication base station radiation detection system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2006132681A (en) * 2006-09-13 2008-03-20 ООО "Лаборатори Информационных Систем" (RU) CORRELATION METHOD FOR RECOGNIZING ADVERSE ELECTROMAGNETIC RADIATION AND CALCULATION OF A COMPUTER ENGINEERING
KR20120101873A (en) * 2011-03-07 2012-09-17 삼성전자주식회사 Prediction apparatus and method for electromagnetic radiated emission of cable
CN103076503A (en) * 2012-12-27 2013-05-01 广东省辐射防护协会 Environmental electromagnetic radiation three-dimensional prediction method of GSM (Global System for Mobile Communication) mobile communication base station
CN104237651A (en) * 2014-09-18 2014-12-24 国家电网公司 Method for calculating radiation power density of antenna of mobile communication base station
CN105676006A (en) * 2016-01-22 2016-06-15 电子信息系统复杂电磁环境效应国家重点实验室 Electronic equipment harmful radiation field prediction method and system based on actual measurement data
CN106454918A (en) * 2016-09-18 2017-02-22 广东智环创新环境技术研究有限公司 Method for predicting electromagnetic radiation level of surrounding environment of TD-LTE mobile communication base station
CN106771668A (en) * 2017-01-05 2017-05-31 西南交通大学 A kind of electromagnetic radiation parameter test system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2006132681A (en) * 2006-09-13 2008-03-20 ООО "Лаборатори Информационных Систем" (RU) CORRELATION METHOD FOR RECOGNIZING ADVERSE ELECTROMAGNETIC RADIATION AND CALCULATION OF A COMPUTER ENGINEERING
KR20120101873A (en) * 2011-03-07 2012-09-17 삼성전자주식회사 Prediction apparatus and method for electromagnetic radiated emission of cable
CN103076503A (en) * 2012-12-27 2013-05-01 广东省辐射防护协会 Environmental electromagnetic radiation three-dimensional prediction method of GSM (Global System for Mobile Communication) mobile communication base station
CN104237651A (en) * 2014-09-18 2014-12-24 国家电网公司 Method for calculating radiation power density of antenna of mobile communication base station
CN105676006A (en) * 2016-01-22 2016-06-15 电子信息系统复杂电磁环境效应国家重点实验室 Electronic equipment harmful radiation field prediction method and system based on actual measurement data
CN106454918A (en) * 2016-09-18 2017-02-22 广东智环创新环境技术研究有限公司 Method for predicting electromagnetic radiation level of surrounding environment of TD-LTE mobile communication base station
CN106771668A (en) * 2017-01-05 2017-05-31 西南交通大学 A kind of electromagnetic radiation parameter test system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHIGANG WANG 等: "Test and research on the electromagnetic environment of electromagnetic radiation from base station", 《 2015 7TH ASIA-PACIFIC CONFERENCE ON ENVIRONMENTAL ELECTROMAGNETICS (CEEM)》 *
周红华 等: "通讯基站50m内轴向发射电磁辐射预", 《海峡科学》 *
许尧 等: "移动通信基站电磁辐射强度的分析与评价", 《通信管理与技术》 *
赵晓峰: "基于MATLAB的电磁辐射数据分析与仿真", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113613282A (en) * 2021-08-12 2021-11-05 广州广电计量检测股份有限公司 Method and device for simulating radio propagation electromagnetic environment of urban scene base station
CN113613282B (en) * 2021-08-12 2023-11-17 广电计量检测集团股份有限公司 Urban scene base station radio propagation electromagnetic environment simulation method and device
CN115078850A (en) * 2022-07-19 2022-09-20 沸蓝建设咨询有限公司 Communication base station radiation detection system

Also Published As

Publication number Publication date
CN111355545B (en) 2021-12-17

Similar Documents

Publication Publication Date Title
Medeđović et al. Wireless insite software verification via analysis and comparison of simulation and measurement results
CN103076503B (en) Environmental electromagnetic radiation three-dimensional prediction method of GSM (Global System for Mobile Communication) mobile communication base station
EP1191804B1 (en) Path loss data normalization for growth management of a cellular system
CN104749447B (en) A kind of evaluation method and device of the electromagnetic radiation from environment of base station
CN111355545B (en) Simplified directional function-based base station electromagnetic radiation visual simulation prediction method
CN103076505B (en) Three-dimensional space prediction method for electromagnetic radiation of TD-SCDMA (Time Division-Synchronization Code Division Multiple Access) mobile communication base station environment
CN103336182B (en) A kind of antenna phase center calibration system based on place insert loss
CN104237651B (en) Antenna of mobile communication base station radiosity computational methods
CN106937300B (en) A kind of method of dynamic and intelligent adjustment Downtilt
CN103412227B (en) Frequency selection radar antenna housing transmission performance testing system and method based on diffraction suppression
JP2010096706A (en) Multi-antenna measurement method and multi-antenna measurement system
CN111257884B (en) Artificial rainfall experiment method for link length adjustable microwave rain attenuation rule
CN107087286A (en) A kind of indoor electromagnetic radiation Forecasting Methodology of window towards communication base station
CN111615138B (en) Strong, medium and weak field region dividing method for mobile base station co-located high-voltage power transmission iron tower
Linhares et al. Determination of measurement points in urban environments for assessment of maximum exposure to EMF associated with a base station
CN105388368A (en) Impedance load point selection method for electromagnetic scattering of high tension overhead transmission lines
Colombi et al. Implications of ICNIRP 2020 exposure guidelines on the RF EMF compliance boundary of base stations
US10750373B1 (en) Determining whether a site supports a supercell base station
Pinem et al. Characterization of Path Loss for Radio Wave Propagation Over the Sea in 4G Network
Sallomi et al. A novel theoretical model for cellular base station radiation prediction
Olivier et al. A practical method for compliance testing of base stations for mobile communications with exposure limits
US7180456B2 (en) Antennas supporting high density of wireless users in specific directions
JP5342219B2 (en) Antenna measurement system and method
JPWO2019181667A1 (en) Radio wave measurement method
CN207164169U (en) The measuring system of extensive MIMO radiation of equipment handling capacities

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