In data-centric wireless sensor networks (WSNs), sensing data have a high time–space correlation. Most queries are spatial and used to obtain data in a defined region. Geographic routing (GR) protocols are the optimal choice for routing spatial queries. However, several drawbacks still exist in GRs, and these the include premature death of nodes and communication latency, which result in reduced network life and query efficiency. A new clustering GR protocol called quadtree grid (QTGrid) was proposed in this study to save energy and improve spatial query efficiency. First, the monitoring area was logically divided into clusters by a quadtree structure, and each grid’s location was encoded to reduce the memory overhead. Second, cluster head (CH) nodes were selected based on several metrics, such as distance from the candidate node to the grid center and adjacent CHs and residual energy. Third, the next-hop routing node was selected depending on the residual energy of the candidate node and its distance to the sink node. Lastly, a lossless data aggregation algorithm and a flexible spatial query algorithm were adopted to reduce the transmission of redundant data and meet the application requirements, respectively. Simulation results showed that compared with three related protocols, QTGrid has lower energy consumption and higher spatial query efficiency and is more suitable for large-scale WSN spatial query application scenarios.