code: def calc_geo_lines_in_buffer(network, node, grid, radius, radius_inc):
      lines = []
      while not lines:
          node_shp = transform(proj2equidistant(network), node.geom)
          buffer_zone_shp = node_shp.buffer(radius)
          for line in grid.graph.lines():
              nodes = line['adj_nodes']
              branch_shp = transform(proj2equidistant(network), LineString([nodes[0].geom, nodes[1].geom]))
              if buffer_zone_shp.intersects(branch_shp):
                  lines.append(line)
          radius += radius_inc
      return sorted(lines, key=lambda _: repr(_))