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(_))