卫星地球站对微波接力站的干扰计算步骤和具体方法分析

由于卫星通信系统与微波接力系统共用频段，会产生同频干扰。为了避免新建卫星地球站建成后的上行频率对地球站周围微波接力系统造成有害干扰，依据中华人民共和国信息产业部第21号令《建立卫星通信网和设置使用地球站管理规定》【2002.6.21】和《建设卫星通信网和设置使用地球站的暂行规定》【国无管（1995.1）4号】，要求建站单位在卫星地球站确定站址之前，在申报有关技术资料的同时，提供数据要准确、论证要充分、分析要科学、结论要明确的卫星地球站对微波接力站的干扰分析报告。无线电管理部门根据建站单位提供的拟建站址电磁环境测试报告和干扰分析报告，作为审批站址时的技术参考。本文主要介绍卫星地球站周围（以卫星地球站天线为中心）半径100 km以内共用频段的卫星地球站对微波接力站的干扰计算步骤和具体方法，供大家参考。

1.相关技术标准

根据GB13616－92“微波接力站电磁环境保护要求”规定，卫星地球站上行信号对微波接力通信系统的干扰允许值为——

（1） 对于模拟微波接力通信系统2500 km假设参考电路的任何话路的相对零电平点，干扰功率应符合以下要求：

a.任何月份20％以上的时间内，噪声计加权1分钟平均功率不超过1000 pW。

b.任何月份0.01％以上的时间内，噪声计加权1分钟平均功率不超过50000 pW。

（2）对于数字微波接力通信系统2500 km假设参考数字通路道的64 kbit/s输出端，干扰功率应符合以下要求：

a.任何月份0.04％以上的时间内，任意1分钟射频干扰功率引起的平均误码率应不超过10-6。

b.任何月份0.0054％以上的时间内，任意1秒钟射频干扰功率引起的平均误码率应不超过10-3。

c.任何月份由于射频干扰功率引起的误码秒累积时间应不超过0.032％。

2.相关参数的计算

在实际的干扰分析计算中，预选卫星地球站站址的经纬度和静止卫星星下点的经纬度及微波接力站（以下简称：微波站）站址的经纬度是已知的，通过这些已知的参数可以计算出：卫星地球站的通信方位角和仰角、卫星地球站与微波站之间的距离、卫星地球站与微波站的空间损耗、卫星地球站鉴别角、微波站鉴别角、卫星地球站天线增益、微波站天线增益、实际干扰信号电平值、干扰信号允许值（干扰容限值）。

2.1 地球站通信方位与仰角的计算

设卫星地球站A的经纬度分别为φ1和θ1，静止卫星W的星下点S的经纬度为Ψ和0，φ0=Ψ-φ1为星下点S对卫星地球站A的经度差，微波站B点的经纬度分别为φ2和θ2。

当卫星地球站A天线对准卫星W时，其仰角EL和方位角AZ分别为：

对于对地静止卫星而言，

式中，RE为地球半径，6378公里;hE为卫星距地面距离大约35786.6公里。

注意：方位角有正有负，此公式是以正南为基准。我国基本上以正北为准，因此：

2.2 卫星地球站和微波站间的通信方位及路径的计算

已知卫星地球站和微波站两点A、B的经纬度，可进行以下计算。

2.2.1 计算A点到B点的真北方位角

A点到B点的真北方位角AZ用下式计算：

当tgAZ 》0，θ1《θ2，方位角为AZ;

若θ1》θ2，方位角为π+AZ;

当tgAZ《0，θ1》θ2，方位角为π-｜AZ｜;

若θ1《θ，方位角为2π-｜AZ｜。

式中： φ1为A点经度，弧度;φ2为B点经度，弧度;

θ1为A点纬度，弧度;θ2为B点纬度，弧度。

2.2.2 计算A点到B点的距离

① 大圆路径计算方法

d=cos-1（cosαcosβ+sinαsinβcosC）*α （7）。

式中，C=ψ2-ψ1; θ1，θ2，ψ2和ψ1的定义同上。

② 平坦地面计算方法（ d《10 km）

式中，R为地球半径，为6378公里; θ1，θ2，ψ2和ψ1的定义同上。

2.3 卫星地球站和微波站干扰鉴别角的计算

卫星地球站通信方向AC与微波站被干扰方向AB之间夹角为ψ，即卫星地球站通信方向对微波站的干扰鉴别角。φA为卫星通信方向在其通信方位角上的投影与该卫星地球站至微波站方向在水平面的夹角。CB垂直于AB，C点向地面投影交于O点（见图1）。

由图2可知：

φA′=｜AZ-AZ′｜ （ rad ） ，

φA =φA′ 0≤φA′≤π ，

或φA=2π－φA′ φA′＞π。

2.3.1 卫星地球站通信方向对微波站被干扰方向的干扰鉴别角的计算

由图1可知： ψ= cos-1 （cosEL* cosφA） （9）。

（以上计算假设微波站工作仰角ε=0时）若ε≠0，则（9）式应写成如下式：

ψ=cos-1（cosε*cosEL*cosφA＋sinε*sinEL）（10）。

2.3.2 微波站通信方向对被卫星地球站 干扰方向的鉴别角的计算

由图2可知：

φB′=|φD-AZ′+π| （ rad ） ，

φB =φB′ 0≤φB′≤π;

或φB= 2π－φB′ φB′＞π。

图2中，BD为有用信号传输路径，AB为干扰信号传输路径， φD为微波站通信方向。

2.4 天线增益的确定

天线增益是通过实际天线方向图来确定的，当缺少天线方向图资料时，可利用下述公式计算天线在某一方向的增益。

设D为天线直径（m）;λ为波长（m）;当D/λ》100时，

其中：Gmax—天线主瓣增益dB，Gmax=10lg［0.6（πD/λ）2］，0.6为天线效率;

G1—天线第一旁瓣增益dB，G1=2+15 lg （D/λ）;

ψ—偏离主波束中心轴的张角，度;

ψ1—主瓣半功率角宽度的一半，可用下式表示：

2.5 传输损耗的计算

2.5.1 视距路径传输损耗计算

视距路径传输损耗主要包括自由空间传输损耗及氧气和水汽的吸收损耗：

Lsd=Lo+（γ0+γw）d （13）。

Lo=92.5+20lgf+20lgd （14）。

式中：Ls—视距路径传输损耗，dB;

Lo—自由空间传输损耗，dB;

f —频率，GHz; d—路径长度，km;

γ0—氧气吸收衰减系数，dB/km;

γw—水汽吸收衰减系数，dB/km。

当f＜15 GHz时，γW＝０，则

式中：ρ—水汽浓度，取决于无线电气候区。

2.5.2 超视距路径传输损耗计算

超视距传播机制主要是绕射（包括障碍物绕射和光滑球面绕射）和对流层散射。对于距离稍超过视距的传输路径，在大多数情况下绕射是主要传播机制，散射可忽略不计;对于远远超过视距的路径，散射是主要传播机制，绕射可忽略不计;对于介于两者之间的中等长度的路径，两种传播机制都需要考虑。考虑参数有：HC—路径余隙，m;HO—自由空间余隙，m。主要传播损耗有四种：

（1）光滑球面绕射损耗。

（2）不规则地形障碍物上的绕射损耗。

在传播路径上往往会遇到一个或多个障碍物，为了估算这些障碍物的附加绕射损耗，通常是将障碍物的形状理想化。当障碍物的厚度相对较窄时，可视为刃形障碍;当障碍物的厚度相对较宽、顶部比较平缓时，可视为圆形障碍。因此又可细分为①刃形绕射损耗的计算;②单圆形障碍物绕射损耗的计算。

（3）多重障碍物绕射计算。

（4）对流层散射传播。

以上四种传播损耗具体计算比较复杂，由于受篇幅的限制在这里就不详细赘述，具体方法见《微波接力通信系统干扰计算方法》（国标GB/T 13619-92）。

2.5.3 不同类型路径传输损耗Lsd的计算

（1）当路径余隙HC＞HO，并且路径距离不超过100 km时，Lsd按自由空间传输损耗计算：

Lsd=Lo=92.5+20lgf+20lgd （17）。

（2）当路径余隙HC≤HO，并且路径距离不超过100 km时，Lsd按自由空间传输损耗计算并附加绕射传输损耗：

Lsd=Lo+Ld （18）。

3.微波接收机输入端实际干扰信号电平值和干扰信号允许值（干扰容限）的计算

3.1 实际干扰信号电平值的计算

实际干扰信号电平值计算公式为：

I=PI-LI+GI+GR-LR-Lsd-XPD+K （19）。

式中，I—微波收信机输入端干扰信号功率，dBm;

PI—干扰发信机输出端信号功率，dBm;

GI—卫星地球站天线在干扰站方向的天线增益，dBi;

GR—微波天线在被干扰源方向的天线增益，dBi;

LR—接收端馈线系统损耗，dB;

LI—干扰站发射端馈线系统损耗，dB;

Lsd—路径传输损耗，dB;

XPD—交叉极化去耦，dB;（由于影响不大可以省略）

信号带宽

K—带宽修正因子，K=10lg ———— 。

干扰带宽

3.2 干扰信号允许值（干扰容限值）的计算

3.2.1 模拟微波站干扰信号允许值（干扰容限值）的计算

模拟TV/FM系统和FDM/FM系统干扰信号允许值的计算比较复杂，需要计算有用信号和干扰信号分别为低、中、高调制指数时的载噪比。本文介绍两种比较简单的工程上常用的计算方法。

A. 采用FM制式的微波接力系统的微波接收机输入端对来自同一频段，采用FDM/FM/FDMA制式的卫星地球站干扰电平允许值采用下述公式估算（用于选址）。

1.任何月份20％以上时间，干扰信号电平允许值为：

I0=10lgTR-176.6 （任一4 kHz带宽） （20）。

2.任何月份0.01％以上时间，干扰信号电平允许值为：

I0=10lgTR-159.6 （任一4 kHz带宽） （21）。

式中：I0—微波接收机输入端干扰信号允许值，dBW;

TR—微波接收机的噪声温度，K。

对于噪声温度TR，可用下式求出：

TR=（F－1）×290 （22）。

式中：F—微波接收机的噪声系数。

B.通过门限恶化量来估算干扰信号允许值（干扰容限值）。门限电平是指恰好满足中断概率指标的接收信号电平。描述门限恶化的一种方法就是给出最大干扰信号允许值。式（23）可用来估算TV/FM系统和FDM/FM系统的最大干扰信号允许值。（23）式表示如下：

I0=Rth-10+Se （23）。

式中：I0—最大干扰信号允许值dBm;

Rth—被干扰的微波接收机门限电平值，dBm;

Se—被干扰的微波接收机有效选择性，dB。

3.2.2 数字微波站干扰信号允许值（干扰容限值）的计算

A.由干扰允许值的误码率计算公式可知：

①MPSK调制系统

Pe≈0.5exp（Eb/N0）=0.5expC/N （DBPSK）;

当M=2时 Pe≈0.5erfc（X） （BPSK）;

当M＞2时 Pe≈erfc［XSin（t）］ （MPSK）。

式中，erfc（Y）—补余误差函数，可用下式表示：

（C/N）0—平均载波功率与平均噪声功率之比，即对应于某一误码率的信噪比理论值;

PL—MQAM 调制系统两个正交分量之一的基带信号误码率;

M—调制电平数;

L＝√M。

B.实际门限信噪比的计算

根据3.2.2 A中公式可算出信噪比的理论值（C/N）0，另外还有设备不完善引起的恶化，系统内部干扰引起的恶化和系统外部干扰引起的恶化；因此，实际的门限信噪比可表达如下：

（C/N）th =（C/N）0 +δ1 +δ2 +δ3 （24）。

式中：（C/N）th—对应于某一误码率（如10-3或10-6）的门限信噪比值，dB;

（C/N）0 —对应于某一误码率（如10-3或10-6）的信噪比理论值，dB;

δ1—设备恶化量（设备厂家给出），dB;

δ2—系统内部干扰的恶化量，dB;

δ3—系统外部干扰的恶化量，dB。

C.信号干扰比允许值计算

信号干扰比允许值可按下式近似估算：

（C/I） 0=（C/N）th+Δ （dB） （25）。

式中：Δ—由于系统外部干扰所要求的信号干扰比增量，dB。

Δ与δ3有如下关系：

Δ=-10lg（100.1δ3-1），

通常取δ3=0.04～0.4 dB，则Δ=10～20 dB。

D.干扰信号允许值（干扰容限值）的计算

I0=C－（C/I）0 （26）。

式中： I0—干扰信号允许值;

C—门限接收电平;

（C/I）0—对应于某一误码率（如10-3或10-6）的信噪比理论值，dB。

４。微波接收机输入端实际干扰信号电平和干扰信号允许值（干扰容限）的比较

经计算可以得出微波接收机输入端实际干扰信号功率I的数值（dBm）和干扰信号允许值I0（dBm）（干扰容限值）。

当卫星地球站周边只要有一个微波站的I≥I0 时，不可以在拟选的站址上架设卫星地球站，卫星地球站会对工作在同频段100 km以内的微波接力站产生干扰。

当I＜I0 时可以在拟选的站址上架设卫星地球站，卫星地球站不会对工作在同频段100 km以内的微波接力站产生干扰。

５。结束语

卫星地球站的电磁环境测试和干扰分析是卫星地球站建设前的一项非常重要的工作。认真做好这项工作，可提高无线电频谱资源的利用率，有效地避免卫星地球站建成后上行信号可能对周边同频段内微波站造成的干扰，降低建设卫星地球站的成本。本文按照“暂行规定”要求，只给出卫星地球站周围（以卫星地球站天线为中心）半径100 km以内的计算方法。在实际的工程设计中，要进行发射卫星地球站协调区计算，画协调区图，协调区内所有同频段的微波站都要进行干扰计算，最后确定协调区内哪几个微波站和卫星地球站间存在有害干扰，尽量避免新建的卫星地球站可能对周边同频段内微波站造成的干扰。