☰
雷达回波为什么是发射波*调制?从物理到仿真的深度解析
2026/10/18 11:29:21 网站建设 项目流程

雷达回波为什么是发射波调制?我第一次在雷达原理的总结材料里看到这个表达式时,心里是先愣了一下:“”是什么意思?乘号?卷积?还是排版错误?等到后来自己动手写脉冲压缩代码、处理真实回波数据时,才明白这句话其实是在用最简练的方式概括整个雷达的观测本质——回波不是凭空产生的另一段电磁波,而是发射信号被目标这个“调制器”改造后的产物,目标的距离、速度、尺度、姿态、甚至细微振动,全都以调制项的形式写进了回波里。

这篇文章想把这件事彻底讲透。我会从物理过程一步步推到信号表达式,再把幅度、时延、多普勒这三类最核心的调制逐项拆开,接着聊波形设计为什么直接影响“解调”能力,最后给出一段可以上手的仿真代码和几个真实工程里最容易踩的坑。适合正在啃雷达原理的电子、通信方向学生,也适合刚接触雷达信号处理的工程师。

1. 先把问题拆开:雷达回波到底是什么

1.1 从发射到回波:电磁波在“途中”经历了什么

要理解“调制”两个字,先得把发射波和回波之间的物理过程走一遍。雷达发射机产生一段射频波形,经过功率放大后由天线定向辐射出去。电磁波离开天线后,会以球面波的形式向空间扩展,照射到目标上。目标本质上是一个由金属或其他导电材料构成的散射体,入射波会在目标表面激励起感应电流,感应电流再向外产生散射波。接收天线真正收到的,是目标朝雷达方向这个特定角度的“后向散射”分量。

这里有一个很多人容易忽略的点:雷达里的所谓“反射”,大多不是光学意义上的镜面反射。镜面反射只在目标表面法线恰好对准雷达时才能形成强回波,而实际目标外形复杂,大量回波来自多个散射中心的叠加散射。雷达更依赖的是分布式散射机制,目标的雷达散射截面(RCS)就是衡量这种散射能力的等效参数。你可以把它理解成:在目标位置放一个各向同性的理想辐射体,它需要多大的等效接收面积、才能产生同样强度的后向散射功率。目标越能“抓住”电磁波再重新辐射出来,RCS就越大,回波自然就越强。

1.2 从物理过程到信号表达式:回波长什么样

把物理过程写成信号,就清晰了。假设发射信号是 s(t),目标在距离 R 处,电磁波往返路径为 2R,所以回波相对发射波有一个双程时延:

τ = 2R / c

其中 c 是光速。如果目标相对雷达还有径向速度 v,回波还会产生一个多普勒频移。一个典型的点目标回波模型可以写为:

s_r(t) = A · s(t − τ) · exp(j2πf_d t) + n(t)

逐项来看:A 是幅度因子,它综合了发射功率、天线增益、波长、目标RCS以及距离四次方衰减;τ 是目标距离带来的时间平移;f_d 是运动目标带来的多普勒频率;n(t) 是接收机噪声和外界干扰。

这个表达式本身就说明了一切:回波仍然是“发射波”,只是被缩放、平移、并且额外乘了一个复指数。目标没有产生新的独立信号,它只是对发射波做了改造。当存在多个目标时,回波就是多个这样的改造项叠加。

更严格地说,对于扩展目标或者宽带信号,回波其实是发射波形与目标冲激响应函数的卷积。卷积在频域中表现为相乘,因此“回波 = 发射波 × 目标传递函数”这句话在频域反而是绝对成立的。而本文讨论的“发射波×调制”,是点目标窄带假设下更直观的时域说法——本质并没有跑偏。

1.3 为什么用“调制”这个词来概括

“调制”原本的含义,是让低频信息去控制高频载波的某个参数,比如调幅广播里声音控制载波幅度。雷达回波恰好也遵循同一套逻辑:发射波就是“载波”,目标的距离、速度、尺寸等物理量就是“信息”,只不过这些信息分别被装进了时间位置、频率和幅度这些参数里。

距离信息调制到时间轴上,表现为整段波形的平移;速度信息调制到频率上,表现为回波载频的偏移;姿态和尺度信息调制到幅度上,表现为回波强弱的起伏。接收机要做的,就是从被调制的回波中把这些参数“解调”出来。

之所以不能把回波简单理解为“发射波加上目标产生的新信号”,是因为目标自身不会发光、不会独立辐射电磁波,它只对入射电磁波产生散射。线性系统里,响应从来都是激励与系统冲激响应的乘积关系,而不是简单叠加。这个区分直接决定了雷达信号处理的框架——你已知发射波形,未知的是目标加在波形上的参数,测距、测速、成像、识别,本质上全是从一个“已知激励”的反向调制问题。

2. 三种调制的底细:目标信息怎么“写”在回波上

2.1 幅度调制:目标的“个头”和姿态

先看最常见也最直觉的幅度项。雷达方程把回波功率写成:

P_r = P_t · G² · λ² · σ / [(4π)³ · R⁴]

这里 P_t 是峰值发射功率,G 是天线增益,λ 是波长,R 是距离,σ 是目标雷达散射截面。在脉冲体制里,平均功率还与脉冲宽度和重复频率有关。

σ 这个参数非常有意思。一个半径为 a 的金属球,RCS 等于 πa²,是一个只跟投影面积有关的几何量。一块平板在法向照射时,RCS 可以达到 4πA²/λ²,这往往比它的几何面积大几个数量级。而某些经过外形设计的隐身目标,RCS 可以被压到 0.001 平方米甚至更低。普通汽车大约有十几平方米的 RCS,小型无人机可能只有 0.01 平方米量级。所以“回波强不强”这个看似简单的问题,背后是目标外形、材料、电尺寸和观测角度的综合结果。

幅度调制还有一个动态维度:目标只要在机动、旋转或者姿态角发生变化,多个散射中心的干涉叠加就会让回波幅度剧烈起伏。目标跟踪里常说的“目标闪烁”,就是这个原因。同样一个飞机,侧面对着雷达时回波可能很强,转个角度就急剧变弱。因此雷达不能只看单次回波幅度就下结论,必须结合多帧积累判断目标是否存在。

要注意的是,幅度 A 里同时混合了距离四次方衰减。距离近的小目标,回波可能比距离远的大目标更强。所以雷达做检测前通常做距离归一化或灵敏度时间控制(STC),把距离带来的大动态范围先压掉一部分,否则中近程目标会把后端的动态范围直接顶满。

2.2 延迟调制:距离信息写在时间轴上

雷达测距最基础的手段就是测延迟。电磁波往返一个距离 R 的目标,总时延是 2R/c。粗略算一下:1 公里往返约 6.7 微秒,10 公里往返约 67 微秒。发射一个脉冲后,只要数一数回波在多久后回来,距离就出来了。

不过“能测延迟”和“能测得多准”是两回事。距离分辨率由信号带宽决定:

δR = c / (2B)

B 是信号带宽。这个结论初看有点反直觉,很多人会以为是脉冲越窄测距越准。没错,普通单频脉冲的带宽约等于脉宽的倒数,窄脉宽自然带宽大。但如果用了调频或者相位编码,就可以在脉宽不变的前提下把带宽做大,所以真正决定距离分辨能力的是带宽,不是脉宽本身。

下面这组数字能更直观地感受带宽的影响:

信号带宽理论距离分辨率
1 MHz150 m
5 MHz30 m
10 MHz15 m
100 MHz1.5 m
1 GHz0.15 m

从“调制”的视角看,目标距离就是把发射波形在时间轴上平移了一个 τ,属于位置调制。接收端要做的事,就是在接收到的长时间序列里寻找“发射波被平移到了哪里”。工程上最常用的做法是匹配滤波,本质是把回波信号与自己复制的发射波形共轭反转后做相关。相关峰出现的位置,就是延迟的估计值。

2.3 多普勒调制:运动目标改写了回波频率

目标运动时,回波会产生多普勒频移。经典计算公式:

f_d = 2v / λ = 2v·f_c / c

v 是径向速度,f_c 是载频。多普勒效应在声学里很好理解,火车鸣笛靠近时音调变高、离开时音调变低。电磁波也一样,只不过频移的大小与载频直接挂钩,载频越高,同样速度对应的多普勒频率越大。

以 X 波段 10 GHz 为例,不同径向速度对应的多普勒频移大概是:

径向速度多普勒频移
10 m/s666.7 Hz
50 m/s3333 Hz
100 m/s6667 Hz
300 m/s20 kHz

相对于 10 GHz 的载频,这些都是很小的相对偏移,几十千赫兹而已。所以雷达接收机必须做到非常低的相位噪声和很高的频率稳定度,否则多普勒信息会被系统自身的相位抖动淹没。

多普勒测速在工程上不是直接“测量回波频率比发射频率高了多少”,而是通过多个脉冲之间的相位变化来计算。相邻两个脉冲的时间间隔是 PRI,如果一个目标以多普勒频率 f_d 运动,那么相邻脉冲间回波相位会累积一个增量:

Δφ = 2π · f_d · PRI

在数字域,对慢时间采样序列做 FFT,得到的频谱峰值位置就是 f_d。所以多普勒测量本质上是对“相位随时间变化率”的估计。这也是为什么雷达要有相参性——非相参雷达只能看包络幅度,看不了相位,就完全没有多普勒能力。

反向运动的目标产生负多普勒,所以接收机必须采用 I/Q 双通道或者复基带处理,保留复数回波的符号信息。只用实数包络,正负速度会混在一起,方向就丢了。

此外,脉冲重复频率 PRF 会带来多普勒模糊。数字采样定律在慢时间维同样成立,可观测的 f_d 范围只有 −PRF/2 到 +PRF/2。超过这个范围,目标多普勒就会折叠到可见区里。想要更高的不模糊速度,就得提高 PRF;但 PRF 提高又会压缩最大不模糊距离。距离模糊和速度模糊的矛盾,是所有脉冲雷达设计者都要面对的第一道关口。

2.4 三种调制同时存在时,怎么拆开看

真实目标通常是三种调制同时作用:距离平移波形,速度偏移频率,RCS 改变幅度,而且可能同时有多个目标,回波是多个分量的叠加。

对付多目标,要靠距离维和多普勒维的二维分辨。快时间维上,匹配滤波能把不同时延的目标分开;慢时间维上,FFT 能把不同多普勒的目标分开。一个目标只要在两个维度上至少有一个分得开,就能被单独检测出来。这也是为什么现代雷达几乎都做“距离—多普勒”二维处理而不是只看一维距离谱。

这里要提一下线性调频信号特有的距离—多普勒耦合。LFM 信号做脉冲压缩后,目标的时延和多普勒都会影响相关峰的位置,表现为峰点在距离—多普勒平面上沿一条斜线移动。目标速度越快,距离测量偏差越大。这个问题等到第三章细说。

所以“解调”从来不是简单看一个幅度峰值,而是要在二维、甚至多维空间里把每个分量的延时和频移联合估计出来。理解这一点,再看雷达数据处理中的各种跟踪滤波算法,就顺理成章了。

3. 波形设计:如何“解锁”调制信息

3.1 单频连续波和简单脉冲的局限

既然回波是被目标调制过的发射波,那“解调”能力的天花板,其实是发射波形自己给的。不同波形的调制可分辨性差别极大。

单频连续波(CW)只发射一个固定频率,频谱是一条极窄的线。这种波形对多普勒极为敏感,测速非常准。但它没有时间上的标记,接收机无法判断回波对应的是哪个时刻发射的,自然也就测不了距离。想用连续波测距,就得加上三角波调频或者其他调制方式,比如调频连续波(FMCW)雷达。

单频脉冲的情况好一些:脉冲有明确的起止时间,可以通过回波延迟测距。但在一个脉冲内部没有任何频率或相位变化,频谱宽度约等于脉宽的倒数,所以速度分辨力取决于脉冲持续时间,距离分辨力又取决于脉冲持续时间——这两者是硬冲突。窄脉冲有好的距离分辨率,但频谱很宽,一个脉冲里能观测到的频率细节少,测速能力弱;宽脉冲频谱窄,测速准,但距离分辨率又差得离谱。

简单脉冲的所有麻烦,都源于它没有“内部调制结构”。想要摆脱束缚,就必须主动给脉冲加上波形调制。

3.2 线性调频脉冲与脉冲压缩

最经典的解决方案是线性调频(LFM)。脉内频率随时间线性变化:

s(t) = rect(t/T) · exp(j2π(f_c t + ½Kt²))

瞬时频率从 f_c 线性增加(或降低)到 f_c + K·T,K 是调频斜率,带宽 B = K·T。这个波形一方面有足够长的脉宽 T,能装下足够的能量;另一方面带宽 B 很大,距离分辨率可以做得很高。时宽带宽积 D = T·B,是这个波形最核心的参数。

匹配滤波(脉冲压缩)后,信号输出是一个主瓣宽度约为 1/B 的 sinc 型窄脉冲,灵敏度峰值比压缩前提高约 T·B 倍。一句话总结:LFM 用脉内频率“调制”换取了大带宽,再用匹配滤波把长脉冲的能量压缩成短脉冲的等效宽度。

举个例子:脉宽 10 微秒、带宽 10 MHz 的 LFM 信号,时宽带宽积是 100。压缩后主瓣宽度等效于 0.1 微秒的脉冲,距离分辨率 15 米,灵敏度增益 20 dB。如果只用 0.1 微秒的方波脉冲去直接获得这个分辨率,峰值功率必须大 100 倍才能有相同的能量。这就是 LFM 能成为现代脉冲雷达“默认选项”的原因。

从这个角度看,雷达波形设计就是设计“调制方式”:让目标携带的信息尽可能以正交、无歧义的形式写进回波。除了 LFM,巴克码相位编码、步进频、频率编码等也都在干同一件事——只不过它们调制的是相位或者频率跳变模式。

3.3 模糊函数:波形的“分辨率身份证”

一个波形能多干净地把目标的延时和频移分辨出来,需要用模糊函数来衡量。模糊函数本质上是发射波形对不同的时延和多普勒组合做匹配相关后输出的响应幅值分布。理想情况下,它应该是一个“图钉”——只有在真实(τ, f_d)处有一个尖锐峰,其余地方都是零,这意味着任何两个不同状态的目标都能被完全分开。

但现实波形总有妥协。单一矩形脉冲的模糊函数在频率维有很宽的响应,说明它的多普勒分辨能力弱。LFM 的模糊函数是一条倾斜的“刀刃”,主峰沿距离—多普勒平面斜向延伸,这正是前面提到的距离—多普勒耦合:只要目标有速度,匹配滤波峰就会沿着刀刃方向偏移,测距结果就会带有一个与速度成比例的系统误差。

所以设计波形时,不能只看理论分辨率参数,最好直接画出模糊函数图。你很快会发现,抑制旁瓣、保持主瓣宽度、减轻距离—多普勒耦合这些需求之间存在冲突,只能针对任务取舍。比如某些成像雷达更在乎旁瓣水平,就会对匹配滤波器加窗加权,牺牲一些峰值分辨能力来换取低旁瓣。

4. 实操:从回波里把调制参数“解”出来

4.1 仿真一把:构造目标回波再测距测速

纸上谈兵不如动手写代码。这里给出一段思路清晰的 Python 仿真:生成一个 LFM 波形,构造一个带延时和多普勒的点目标回波,加上噪声后先做匹配滤波测距离,再做慢时间 FFT 测速度。

import numpy as np from scipy.signal import fftconvolve fs = 100e6 # 快时间采样率 100 MHz T = 10e-6 # 脉宽 10 us B = 10e6 # 带宽 10 MHz K = B / T # 调频斜率 pri = 200e-6 # 脉冲重复间隔 200 us n_pulse = 128 # 慢时间脉冲数 fc = 10e9 # 载频 10 GHz c = 3e8 # 目标:距离 3000 m,径向速度 -80 m/s(负号代表接近雷达) R0 = 3000.0 v = -80.0 fd = 2 * v / c * fc N_fast = int(pri * fs) t_fast = np.arange(N_fast) / fs # 发射LFM(只在前 T 秒内有信号) s_tx = np.zeros(N_fast, dtype=complex) idx = np.arange(int(T * fs)) s_tx[:len(idx)] = np.exp(1j * np.pi * K * (idx / fs) ** 2) # 时延对应的采样点数 delay_samp = int(round(2 * R0 / c * fs)) # 构造128个脉冲的回波 data = np.zeros((n_pulse, N_fast), dtype=complex) for m in range(n_pulse): # 窄带近似下,脉内多普勒相位变化很小,取该脉冲开始时刻的相位即可 data[m, delay_samp:delay_samp + len(idx)] = ( s_tx[:len(idx)] * np.exp(1j * 2 * np.pi * fd * m * pri) ) # 加复高斯白噪声 noise_amp = 0.05 data += noise_amp * (np.random.randn(n_pulse, N_fast) + 1j * np.random.randn(n_pulse, N_fast)) # 匹配滤波(脉冲压缩) h = np.conj(s_tx[::-1]) pc = np.array([fftconvolve(data[m], h, mode='same') for m in range(n_pulse)]) # 距离维峰值:对慢时间平均后取最大距离门 rng_profile = np.mean(np.abs(pc), axis=0) peak_bin = np.argmax(rng_profile) # 该距离门上的慢时间序列做FFT测多普勒 slow_sig = pc[:, peak_bin] spec = np.fft.fftshift(np.fft.fft(slow_sig)) freq_axis = np.fft.fftshift(np.fft.fftfreq(n_pulse, pri)) # 谱峰换算速度 peak_fd = freq_axis[np.argmax(np.abs(spec))] v_est = peak_fd * c / (2 * fc) print(f"峰值距离门: {peak_bin}, 对应距离约 {peak_bin / fs * c / 2:.1f} m") print(f"多普勒谱峰频率: {peak_fd:.1f} Hz") print(f"估计径向速度: {v_est:.2f} m/s")

运行起来,你会看到距离维匹配滤波出现一个明显的峰值,所在距离门就是 3000 米附近;慢时间 FFT 的谱峰位置对应多普勒频率约 −5333 Hz(10 GHz 下 −80 m/s 的目标),换算速度后非常接近设定值。简单说,距离解出来了靠“脉冲在时间上的平移”,速度解出来了靠“相位随时间的变化”,两次“解调”分别完成。

我把噪声幅度调到 0.05 时,单脉冲距离剖面还能比较清晰地看到峰。如果把噪声加到 0.5 再跑,单次脉冲就完全找不到峰值了——这时就体现出了慢时间积累的价值:128 个脉冲做 FFT 之后,又多了大约 21 dB 的信噪比增益,目标又从噪声里“浮”出来了。这也是相参积累的本质含义。

4.2 真实雷达数据怎么读:A显、R显和距离—多普勒图

仿真看完,回到真实界面。雷达显示方式本质上都是回波调制参数的投影。

A 显最简单:横轴是时间(或距离),纵轴是回波幅度。一个目标就是一段距离上的一个包络凸起。A 显只能看到幅度和距离,看不到速度,适合观察单个脉冲的原始回波。

R 显(距离—方位显示)把距离和方位角画成极坐标,亮点的径向位置表示距离,角度位置表示方位,亮度或颜色表示回波强度。这是监视雷达、航管雷达最常见的画面。你看到的每个飞机回波,都对应一个或多个距离分辨单元里的幅度调制峰值,只是屏幕又把幅度映射成了亮度。

距离—多普勒图是最直观体现“调制”的二维图。横轴是快时间距离,纵轴是慢时间多普勒频率,每个像素的亮度代表该距离、该速度上有多少能量。静止的地物杂波会压在多普勒零频这条线上;运动目标会出现在偏离零频的位置;如果目标存在微动部件,还会在主体峰周围看到周期性的边带分量。换句话说,你在距离—多普勒图上看到的每一个亮点,背后正是“发射波被目标调制后的两个核心参数在二维空间的落点”。

4.3 必须避开的坑

仿真和真实数据之间,有一批典型的坑,几乎每个人都会踩:

距离模糊。PRF 太高,远处目标的回波会延迟到下一个甚至下下个脉冲周期里回来,于是它被显示在比真实距离近得多的地方。消除办法是使用参差 PRI 或者多 PRF 解模糊,但代价是处理复杂度上升。

速度模糊。PRF 太低,高速目标的多普勒超过 ±PRF/2,频谱折返后看起来像个低速目标,甚至方向都反了。多普勒模糊在雷达里有个直观名字叫“多普勒混叠”。

盲速。当目标多普勒频率恰好是 PRF 的整数倍时,慢时间相位每脉冲都旋转整周,看起来完全静止。这个现象处理不好,高速目标会被当成静止杂波滤掉。

LFM 距离—多普勒耦合。快速目标测距离会出现系统性偏差,必须先用粗测速度修正距离。在成像雷达里,这个校正通常是运动补偿的一部分。

还有匹配滤波本身的问题:发射信号经过功放和天线之后会有一点失真,如果复制信号和实测发射波形不完全一致,脉压会出现失配损失和旁瓣恶化。最好用定向耦合器实时采集一版发射波样本,作为模板做数字脉压,而不是拿着理想波形函数硬凑。

5. 常见误解与避坑速查

5.1 误解一:回波波形只是发射波延迟后的“完全副本”

点目标窄带模型确实把这个近似用到了极限,很多教材开头就这么讲,但它带偏了不少初学者。真实目标是分布式的,由多个散射中心组成,每个散射中心本身又可能具有距离上的分散和不同的姿态变化。宽带雷达照射一架飞机时,回波里叠加了来自机头、机翼、发动机进气道等位置的多个散射分量,它们的相对相位随观测角变化,于是合成回波的幅度、相位甚至在极短时间内剧烈变化。再加上极化响应,雷达在目标识别上能看到的信息其实非常丰富。

一个典型例子是微多普勒。直升机旋翼、喷气发动机叶片、车轮旋转都会给回波叠加周期性的频率调制,频谱上表现为主峰附近出现成簇的边带。这些特征常被用来做机型识别、步态识别、甚至区分飞鸟和无人机。如果你只盯着“延迟副本”的模型,根本想不到去看这些边带。

5.2 误解二:多普勒只有高速目标才明显

微多普勒恰恰来自很慢的运动。人类呼吸时胸壁的起伏速率只有每秒几次,幅度也不过毫米级,但在生物雷达的频域观测量里,呼吸频率、心跳频率都清晰可见。穿墙雷达就是靠提取这类微多普勒信号,判断墙后是否有人以及人在做什么动作。地震救援、老年人跌倒检测、驾驶疲劳监测都用到了同一套原理:回波属于发射波的调制项,不管目标动得多慢、动得多微小,只要相位能被相参系统捕捉到,就能暴露出来。

我最初调试非相参模拟器的数据时,一直读不出微动特征,后来意识到根本不是信噪比不够,而是硬件链路不保持相参,回波相位在两次采样之间随机跳变,微多普勒信息早就被系统自身的不确定性抹掉了。想观察弱调制,第一步是保证雷达系统自己是相参的、相位噪声足够低,这个前提比功率更重要。

5.3 误解三:回波弱就加发射功率

雷达方程告诉我们,回波功率与发射功率成正比,但现实中“看得清”需要的往往不是单一功率,而是整个链路的能力。提高峰值功率会遇到器件成本、电磁兼容、辐射安全等多重限制;更容易被忽视的是积累增益:把 128 个脉冲相参积累起来,等效信噪比提升约 21 dB,这比把峰值功率提升 100 倍还划算。很多雷达设计里,真正决定探测距离上限的是系统相位噪声、杂波对消能力、接收机噪声系数和匹配滤波增益,而不是发射机有多“大力”。

结合常见误区,整理一个速查表供日常排查使用:

现象直觉判断通常正解
回波弱发射功率不够先查相参积累、系统噪声系数、相位噪声
目标距离测不准采样率不够先看信号带宽和脉压失配、距离—多普勒耦合
速度方向不对信号极性接反检查多普勒模糊范围和 I/Q 通道一致性
脉冲压缩后旁瓣高波形问题先排查模板失配、加窗函数、非线性放大失真
目标时有时无目标真的不存在先确认RCS闪烁、遮挡、姿态角变化,再谈灵敏度

最后分享一点我自己的体会。刚接触雷达信号处理那会儿,我总觉得回波就是“把发射波延迟一下、放大或缩小一下”,测距仿真做得顺风顺水,一到测速就怎么都对不上。后来才发现问题出在没用复基带信号、丢了相位信息,多普勒调制根本不可能从纯实数包络里解出来。从那时起,我习惯把所有雷达问题都拆成三个维度来想:幅度调制对应“目标有多大”,时延调制对应“目标有多远”,频率调制对应“目标动多快”。一旦用这个框架去审视,阅读雷达论文、理解新算法都顺畅多了。雷达回波是发射波乘以目标调制,这句话不是什么玄学公式,它只是把整个雷达工作的本质,压进了一行字里。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询