网上给出的计算公式是这样的

Lng1表示A点纬度和经度，Lat2 Lng2 表示B点纬度和经度
a = Lat1 – Lat2为两点纬度之差 b = Lng1 -Lng2 为两点经度之差
6378.137为地球半径，单位为公里
计算出来的结果单位为公里

由大牛改装过的代码是这样的：
Mysql排序公式一：

select * from sb_user_geo order by ACOS(SIN(('.$city_lat.' * 3.1415) / 180 ) *SIN((lat * 3.1415) / 180 ) +COS(('.$city_lat.' * 3.1415) / 180 ) * COS((lat * 3.1415) / 180 ) *COS(('.$city_lng.' * 3.1415) / 180 - (lng * 3.1415) / 180 ) ) * 6380  asc  limit 30

city_lat,city_lng 分别变量经纬度，lat，lng为数据库资源

Mysql排序公式二：

select *,(2 * 6378.137* ASIN(SQRT(POW(SIN(PI()*(111.86141967773438-lng)/360),2)+COS(PI()*33.07078170776367/180)* COS(lat * PI()/180)*POW(SIN(PI()*(33.07078170776367-lat)/360),2)))) as juli from area  order by juli asc limit 0,20

Mysql 语句，根据经纬度，计算5公里范围内的数据：

select ROUND(6378.138*2*ASIN(SQRT(POW(SIN(($latitude*PI()/180-latitude*PI()/180)/2),2)+COS($latitude*PI()/180)*COS(latitude*PI()/180)*POW(SIN(($longitude*PI()/180-longitude*PI()/180)/2),2)))*1000) AS distance FROM shop having distance <= 5000 order by distance asc

PHP计算两点之间的距离：

/**
 *求两个已知经纬度之间的距离,单位为千米
 *@param lng1,lng2 经度
 *@param lat1,lat2 纬度
 *@return float 距离，单位千米
 **/
private function _distance($lng1,$lat1,$lng2,$lat2)//根据经纬度计算距离
{
    //将角度转为弧度
    $radLat1=deg2rad($lat1);
    $radLat2=deg2rad($lat2);
    $radLng1=deg2rad($lng1);
    $radLng2=deg2rad($lng2);
    $a=$radLat1-$radLat2;//两纬度之差,纬度<90
    $b=$radLng1-$radLng2;//两经度之差纬度<180
    $s=2*asin(sqrt(pow(sin($a/2),2)+cos($radLat1)*cos($radLat2)*pow(sin($b/2),2)))*6378.137;
    return $s;
}

https://www.cnblogs.com/wrld/p/9377680.html
https://www.cnblogs.com/atjs/p/5053005.html